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A review of CP violation from the Standard Model to strings is given which includes a broad 
landscape of particle physics models, encompassing the non-supersymmetric 4D extensions of the 
standard model, and models based on supersymmetry, on extra dimensions, on strings and on 
branes. The supersymmetric models discussed include complex mSUGRA and its extensions, 
while the models based on extra dimensions include 5D models including models based on warped 
geometry. CP violation beyond the standard model is central to achieving the desired amount of 
baryon asymmetry in the universe via baryogenesis and leptogenesis. They also affect a variety of 
particle physics phenomena: electric dipole moments, g — 2, relic density and detection rates for 
neutralino dark matter in supersymmetric theories, Yukawa unification in grand unified and string 
based models, and sparticle production cross sections, and their decays patterns and signatures at 
hadron colliders. Additionally CP violations can generate CP even-CP odd Higgs mixings, afi'ect 
the neutral Higgs spectrum and lead to phenomena detectable at colliders. Prominent among 
these are the CP violation effects in decays of neutral and charged Higgs bosons. Neutrino masses 
introduce new sources of CP violation which may be explored in neutrino factories in the future. 
Such phases can also enter in proton stability in unified models of particle interactions. The 
current experimental status of CP violation is discussed and possibilities for the future outlined. 
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I. INTRODUCTION 

Wc begin with a brief history of the considerations 
that led to question the vahdity of CP symmetry as an 
exact symmetry for elementary particles. The history is 
tied to the issue of electric dipole moments and we need 
to retrace the steps back to 1950 when it was generally 
accepted that the particle electric dipole moments 
vanished due to parity symmetry. However, in 1950 
it was first observed by Purcell and Ramsey (Purcell 
and Ramsey, 1950), that there was no experimental 
evidence for the parity symmetry for nuclear forces and 
for elementary particles, and thus the possible existence 
of an electric dipole moment for these needed to be 
tested experimentally. They and their graduate student 
James Smith then carried out the first such test by 
showing experimentally in 1951 that the magnitude of 
the electric dipole moment of the neutron was less than 
3 X 10~^° e.cm where e is the charge of the proton ^. 
After the violation of parity symmetry proposed by T.D. 
Lee and C.N. Yang (Lee and Yang, 1956) was confirmed 
(Wu et al., 1957), it was argued by many that the 
elementary electric dipole moments would vanish due to 
the combined charge conjugation and parity symmetry, 
i.e., CP symmetry (or cquivalcntly under a time reversal 
symmetry under the assumption of CPT invariance). 
However, it was then pointed out by Ramsey(Ramsey, 
1958) and independently by .Jackson and coUabora- 
tors(Jackson et aZ., 1957) that T invariance was also an 
assumption and needed to be checked experimentally (A 
brief review of early history can be found in (Ramsey, 
1998)). Since then the search for CP violations has been 



The experimental results of Purcell, Ramsey and Smith while 
completed in 1951 were not published till much later(Smith et al., 
1957). However, they were quoted in other publications(Lee and 
Yang, 1956; Ramsey, 1956; Smith, 1951). 



vigorously pursued. The CP violation was eventually 
discovered in the Kaon system by Val Fitch, James 
Cronin and collaborators in 1964(Christenson et al., 
1964). Shortly thereafter it was pointed out by Andre 
Sakharov(Sakharov, 1967) that CP violations play an 
important role in generating the baryon asymmetry in 
the universe. However, it has recently been realized 
that sources of CP violation beyond what exist in the 
Standard Model arc needed for this purpose. In this 
context over the past decade a very significant body 
of work on CP violation beyond the Standard Model 
has appeared. It encompasses non-supersymmetric 
models, supersymmetric models, models based on extra 
dimensions and warped dimensions, and string models. 
There is currently no review which encompasses these 
developments. The purpose of this review is to bridge 
this gap. Thus in this review we present a broad 
overview of CP violation starting from the Standard 
Model and ending with strings. CP violation is central 
to understanding the phenomena in particle physics as 
well as in cosmology. Thus CP violation enters in K and 
B physics, and as mentioned above CP violation beyond 
the Standard Model is deemed necessary to explain the 
desired baryon asymmetry in the universe. Further, new 
sources of CP violation beyond the Standard Model 
could also show up in sparticle production at the LHC, 
and in the new generation of experiments underway 
on neutrino physics. In view of the importance of CP 
violation in particle physics and in cosmology it is 
also important to explore the possible origins of such 
violations. These topics arc the focus of this review. 
We give now a brief outline of the contents of this review. 

In Sec. (II) we give a discussion of CP violation in 

the Standard Model and of the strong CP problem. 
The electroweak sector of the Standard Model contains 
one phase which appears in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. The CKM matrix satisfies 
unitarity constraints including the well known unitarity 
triangle constraint where the three angles a, (3, 7 defined 
in terms of ratios involving the products of CKM matrix 
elements and their complex conjugates sum to tt. In 
addition the quantum chromo dynamic (QCD) sector 
of the Standard Model brings in another source of CP 
violation - the strong CP phase Oqcd- The natural 
size of this phase is 0(1) which would produce a huge 
contribution to the electric dipole moment (EDM) of 
the neutron in violation of the existing experimental 
bounds. A brief discussion of these issues is given in 
Sec. (II). A review of the experimental evidence for 
CP violation and of the searches for evidence of other 
CP violation such as in the electric dipole moment of 
elementary particles and of atoms is given in Sec. (III). 
Here we discuss the current experimental situation in the 
K and B system. In the Kaon system two parameters 
e (indirect CP violation) and e' (direct CP violation) 
have played an important role in the discussion of CP 
violation in this system. Specifically the measurement 
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of e'/e rules out the so called superweak theory of CP 
violation while the measurement is consistent with the 
Standard Model prediction. In this section we also give 
an analysis of experimental constraints on the angles 
a, /3, 7 of the unitarity triangle discussed in Sec. (II). The 
current experimental limits of the EDMs of the electron, 
of the neutron and of ^^^Hg are also discussed. 

In Sec. (IV) we give a discussion of the CP viola- 
tion in some non-supersymmetric extensions of the 
Standard Model. These include the Left-Right (LR) 
extensions, the two Higgs doublet model and extensions 
with more than two Higgs doublets. It is shown that 
such extensions contain more sources of CP violation. 
For example, the LR extensions with the gauge group 
SU{2) L X SU{2)iiX U{1)y and three generations contains 
seven CP phases instead of one phase that the Standard 
Model has. Similarly it is shown that the number of 
allowed CP phases increases with the number of Higgs 
doublets. Further, new sources of CP violation arise as 
one increases the number of allowed generations. CP 
violation in the context of supersymmetric extensions 
of the Standard Model are discussed in Scc.(V). Here 
one finds that the minimal supersymmetric standard 
model (MSSM) has a large number (i.e., 46) of phases 
which, however, is reduced to two phases in the minimal 
supergravity unified model (mSUGRA). However, more 
phases are allowed if one considers supergravity unified 
models with non-univcsal soft breaking at the grand 
unified (GUT) scale consistent with flavor changing 
neutral current (FCNC) constraints. A discussion of CP 
violation in extra dimension models is given in Sec. (VI). 
In this section we give an exhibition of the phenomena 
of spontaneous vs explicit CP violation. In this section 
we also discuss CP violation in the context of warped 
extra dimensions. 

A discussion of CP violation in strings is given in 
Sec. (VII). It is shown that soft breaking in string models 
is parametrized by vacuum expectation values (VEVs) of 
the dilaton (5) and of the moduli fields (Tj) which carry 
CP violating phases. Additionally CP phases can oc- 
cur in the Yukawa couplings. Thus CP violation is quite 
generic in string models. We give specific illustration of 
this in a Calabi-Yau compactification of an Es x het- 
erotic string and in orbifold compactifications. Here we 
also discuss CP violation in D brane models. Finally in 
this section we discuss the possible connection of SUSY 
CP phases with the CKM phase. 

A discussion of the computation of the EDM of an 
elementary Dirac fermion is given in Sec. (VIII) while 
that of a charged Icpton in supersymmetric models is 
given in Sec. (IX). In Scc.(X) we give an analysis of the 
EDM of quarks in supersymmctry. The supersymmetric 
contributions to the EDM of a quark involve three differ- 
ent pieces which include the electric dipole, the chromo 
electric dipole and the purely gluonic dimension six oper- 
ators. The contributions of each of these are discussed in 



Sec.(X). Typically for low lying sparticle masses the su- 
persymmetrtic contribution to the EDM of the electron 
and of the neutron is generally in excess of the current 
experimental bounds. This poses a serious difficulty for 
supersymmetric models. Some ways to overcome these 
are also discussed in Sec.(X). Two prominent ways to 
accomplish this include either a heavy sparticle spectrum 
with sparticle masses lying in the TeV region, or the 
cancellation mechanism where contributions arising 
from the electric dipole, the chromo electric dipole and 
the purely gluonic dimension six operators largely cancel. 



If the large SUSY CP phases can be made consistent 
with the EDM constraints, then such large phases can 
affect a variety of supersymmetric phenomena. We 
discuss several such phenomena in Sec. (XI). These 
include analyses of the effect of CP phases on — 2, 
on CP even-CP odd Higgs mixing in the neutral Higgs 
sector, and on the b quark mass. Further, CP phases 
can affect significantly the neutral Higgs decays into 
bb and rf and the decays of the charged Higgs into 
tbjDrT and the decays X^X°- These phenomena 

are also discussed in Sec. (XI). Some of the other 
phenomena affected by CP phases include the relic 
density of ncutralino dark matter, proton decay via 
dimension six operators, the decay A*^A*~, decays 

of the sfermions, and the decay B — > (f)K. These are 
all discussed in some detail in Sec. (XI). Finally in this 
section we discuss the T and CP odd operators and their 
observability at colliders. An analysis of the interplay 
between CP violation and fiavor is given Sec. (XII). 
Here we first discuss the mechanisms which may allow 
the muon EDM to be much larger than the electron 
EDM, and accessible to a new proposed experiment on 
the muon EDM which may extend the sensitivity of 
this measurement by several orders of magnitude and 
thus make it potentially observable. In this section 
an analysis of the effect of CP phases on B — > Xgj 
is also given. This FCNC process is of importance as 
it constrains the parameter space of MSSM and also 
constrains the analyses of dark matter. Sec. (XIII) is 
devoted to a study of CP violation in neutrino physics. 
Here a discussion of CP violation and leptogenesis is 
given, as well as a discussion on the observability of 
Majorana phases. 



Future prospects for improved measurement of CP 
violation in experiments are discussed in Sec. (XIV). 
These include improved experiments for the measure- 
ments of the EDMs, B physics experiments at the LHCb 
which is dedicated to the study of B physics, Super Belle 
proposal, as well as superbeams which include the study 
of possible CP violation in neutrino physics. Conclu- 
sions are given in Sec. (XV). Some further mathematical 
details are given in the Appendices in Sec. (XVI). 
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II. CP VIOLATION IN THE STANDARD MODEL AND 
THE STRONG CP PROBLEM 

The electroweak sector of the Standard Model with 
three generations of quarks and leptons has one CP vi- 
olating phase which enters via the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix V. Thus the electroweak inter- 
actions contain the CKM matrix in the charged current 
sector 

92ua^,Vij{^ - ^5)djW^ + H.c. (1) 

where Ui = u,c,t and dj = d,s,b quarks. The CKM 
matrix obeys the unitarity constraint {VV^)ij = Sij and 
can be parameterized in terms of three mixing angles and 
one CP violating phase. For the case i j the unitarity 
constraint can be displayed as a unitariy triangle, and 
there are six such unitarity triangles. Thus the unitarity 
of the CKM matrix for the first and the third column 
gives 

VudV:, + VcdV:, + VtdV,l = 0. (2) 

One can display this constraint as a unitarity triangle by 
defining the angles a, /3, 7 so that 

a = arg{~VtdV;jVudV:,), (3 = arg{-VedV:jVtdV;,), 

7 = arg{-VudV:jV,dV:,) (3) 

which satisfy the constraint a + (3 + j = n. One can 
parameterize CP violation in a way which is independent 
of the phase conventions. This is the so called Jarlskog 
invariant (Jarlskog, 1985) J which can be defined in nine 
different ways, and one of which is given by 

J = /m(KsKrbK5K*s)- (4) 

An interesting observation is that the CKM is hierarchi- 
cal and allows for expansion in A ~ 0.226 so one may 
write y as a perturbative expansion in A which up O(A^) 

is given by 

/ 1-f A AX^p-t7^)\ 

-A i_2f AA2 (5) 

\AX^{l-p-iri) -A)? 1 / 

In this representation the Jarlskog invariant is given by 

J ~ A'^X^T], and the CP violation enters via 77. 

The Standard Model has another source of CP vio- 
lation in addition to the one that appears in the CKM 
matrix. This source of CP violation arises in the strong 
interaction sector of the theory from the term O^GG, 
which is of topological origin. It gives a large contribu- 
tion to the EDM of the neutron and consistency with 
current experiment requires 6 = 6 + ArgDet{MuMd) to 
be small 9 < O(10^^°). One solution to the strong CP 
problem is the vanishing of the up quark mass. However, 
analyses based on chiral perturbation theory and on 
lattice gauge theory appear to indicate a non-vanishing 
mass for the up quark. Thus a resolution to the strong 



CP problem appears to require beyond the Standard 
Model physics. For example, one proposed solution 
is the Peccei-Quinn mechanism (Pcccei and Quinn, 
1977) and its refinements(Dine et al., 1981; Kim, 1979; 
Zhitnitskii, 1980) which leads to axions. But currently 
severe limits exist on the corridor in which axions can 
exist. There is much work in the literature regarding 
how one may suppress the strong CP violation effects 
(for a review see (Dine, 2000)). In addition to the use of 
axions or a massless up quark one also has the possibility 
of using a symmetry to suppress the strong CP effects 
(Barr, 1984; Nelson, 1984). 



The solution to the strong CP in the framework of 

Left-Right symmetric models is discussed in (Babu et al., 
2002; Mohapatra et al., 1997). Specifically in the analy- 
sis of (Babu et al., 2002) the strong CP parameter 6 is 
zero at the tree level, due to parity (P), but is induced 
due to P -violating effects below the unification scale. 
In the analysis of (Hiller and Schmaltz, 2001) a solution 
to the strong CP problem using supersymmetry is pro- 
posed. Here one envisions a solution to the strong CP 
problem based on supersymmctric non-renormalization 
theorem. In this scenario CP is broken spontaneously 
and its breaking is communicated to the MSSM by radia- 
tive corrections. The strong CP phase is protected by a 
SUSY non-renormalization theorem and remains exactly 
zero while the loops can generate a large CKM phase 
from wave function renormalization. Another idea ad- 
vocates promoting the U{1) CP violating phases of the 
supersymmctric standard model to dynamical variables, 
and then allowing the vacuum to relax near a CP con- 
serving point (Dimopoulos and Thomas, 1996). In the 
analysis of (Demir and Ma, 2000) an axionic solution of 
the strong CP problem with a Peccei-Quinn mechanism 
using the gluino rather than the quarks is given and the 
spontaneous breaking of the new U(l) global symmetry is 
connected to the supersymmetry breaking with a solution 
to the n problem (Demir and Ma, 2000). Finally, in the 
analysis of (Aldazabal et al., 2004) a solution based on 
gaiiging away the strong CP problem is proposed. Thus 
the work of (Aldazabal et al., 2004) proposes a solution 
that involves the existence of an unbroken gauged U {l)x 
symmetry whose gauge boson gets a Stueckelberg mass 
term by combining with a pseudoscalar field 77(2;) which 
has an axion like coupling to GG. Thus the 6 parameter 
can be gauged away by a C/(l)x transformation. The ad- 
ditional U{l)x generates mixed gauge anomalies which 
are canceled by the addition of an appropriate Wess- 
Zumino term. We will assume from here on that the 
strong CP problem is solved by one or the other of the 
techniques outlined above. 
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III. REVIEW OF EXPERIMENTAL EVIDENCE ON CP 
VIOLATION AND SEARCHES FOR OTHER EVIDENCE 

There axe currently four pieces of experimental evi- 
dence for CP violation. These consist of (i) the obser- 
vation of indirect CP violation (c), and (ii) of direct CP 
violation (e'/e) in the Kaon system, (iii) the observation 
of CP violation in B physics, and (iv) an indirect evi- 
dence for CP violation due to the existence of baryon 
asymmetry in the universe. Thus far the experimental 
evidence indicates that the CP violation in the K and B 
physics can be understood within the framework of the 
standard model. However, an understanding of baryon 
asymmetry in the universe requires a new source of CP 
violation We briefly review these below. 
A. CP violations in the Kaon system^ 
Historically the first indication for CP violation came 
from the observation of the decay tt+tt". In order 

to understand this phenomenon we begin with the states 
K'^ (with strangeness S = +1) and K'^ (with strangeness 
S = —1). From the above one can construct CP even 
and CP odd eigenstates. 



(6) 



One can arrange to be the CP conjugate of K'^, 
i.e., CP\K" >= \K" >, and in that case Ki is the CP 
even and K2 is the CP odd state. The decay of neu- 
tral K's come in two varieties: Ks{Kl) with lifetimes 
Ts = 0.89 X 10^°s(tl = 5.2 x 10"**) with dominant de- 
cays Ks ,Tr'^n^{KL 3tt,ttIv). If these were 
the only decays one would identify Ks with Ki and Kl 
with K2- However, the decay of the Kl tt+tt" pro- 
vided the first experimental evidence for the existence of 
CP violation (Cliristenson et ai, 1964). This experiment 
indicates that the Ks{Kl) are mixtures of CP even and 
CP odd states and one may write 



Ks 



Ki + eK2 K2 + eKi _ 



Experimentally one attempts to measure two indepen- 
dent CP violating parameters e and e' which are defined 
by 

^ < {Tni)i=Q\Cw\KL > . . 

^~ <{ni^)i=o\Cw\Ks>' ^' 

where Cvm is the Lagrangian for the weak AS' = 1 inter- 
actions, and 

, ^ < {■mT)i=2\Cw\KL > _ < in'K)i=2\Cw\Ks > . . 
^ < {7nr)i^o\Cw\KL > < {mr)i=o\Cw\Ks >' ^ ' 



^ For a review of this topic see (Bertolini et al., 2000; Winstein 
and Wolfenstein, 1993). 



The parameter e' is often referred to as a measure of 
direct CP violation while e is referred to as a measure of 
indirect CP violation in the Kaon system. An accurate 
determination of e has existed for many years so that 



|e| = (2.266 ±0.017) x 10 



-3 



(10) 



The determination of direct CP violation is more recent 
and here one has (Alavi-Harati et al, 1999; Burkhardt 
et al, 1988; Fanti et al, 1999) 



e7e= (1.72 ±0.018) x 10" 



(11) 



The above result rules out the so called superweak theory 
of CP violation (Wolfenstein, 1964) but is consistent with 
the predictions of the Standard Model. A detailed dis- 
cussion of direct CP violation can be found in (Bertolini 
et al., 2000). There are other Kaon processes where CP 
violation effects can, in principle, be discerned. The most 
prominent among these is the decay Kl tt^vP. This 
process is fairly clean in that it provides a direct deter- 
mination of the quantity VtdV*g. The Standard Model 
prediction for the branching ratio is (Buras et al., 2004) 
BR{Kl ^ TT^iyP) = (3.0 ± 0.6) X 10-11 while the cur- 
rent experimental limit is (Anisimovsky et al., 2004) 
BR{Kl 'K^vu) < 1.7 X 10"^. Thus an improvement 
in experiment by a factor of around 10^ is needed to test 
the Standard Model prediction. On the other hand sig- 
nificantly larger contribution to this branching ratio can 
arise in beyond the Standard Model physics (Buras et al., 
2005, 2004; Colangelo and Isidori, 1998; Grossman and 
Nir, 1997). A new experiment, 391a, is underway at KEK 
which would have a significantly improved sensitivity for 
the measurement of this branching ratio and its results 
could provide a window to testing new physics in this 
channel. 

We turn now to B physics. There is considerable lit- 
erature in this area to which the reader is directed for 

details ((Bigi and Sanda, 1981, 1984; Carter and Sanda, 
1980; Dunietz and Rosner, 1986). For reviews see (Bar- 
berio, 1998; Harrison and Quinn, 1996; Hitlin and Stone, 
1991; Nakada, 1994; NarduUi, 1993; Peruzzi, 2004; Quinn, 
1998; Sanda, 2004; Stone, 2006)). CP violations can oc- 
cur in charged B or neutral B decays such as B^ = bd 
and Bg =bs. In the B° — B° system the mass eigenstates 
can be labeled as Bh and Bl with 



\Bl>=p\B° >+q\B° >, 
\Bh >=p|B°>-g|SO >, 

where p{q) may be parameterized by 

l±eB 



(12) 



P = 



V2(l + |eBP)' 
V2(l + |eBP)' 



(13) 



A quantity of interest is the mass difference between these 
states, i.e., Am^ = uibh ~ ''^Bl- Next let us consider a 
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state f which is accessible to both i?" and B". A quan- 
tity sensitive to CP violation is the asymmetry which is 
defined by 

a (t) = r(i3°(t)^/)-r(^°(f)^/) 

r(i?"(i)^/) + r(BO(t)^/) ^ ' 

where {B^{t)) denote the states which were ini- 

tially B^{B^). The analysis of the asymmetry becomes 
specially simple if the final state is an eigen state of CP. 
Af{t) may be written in the form 



Aj{t) = Ay cos(Ami) + Ay sin(Amt) 



where 



A 



-2/toA 



/ 



1 



lAI 



(15) 



(16) 



Here A ee qAf/pAf, where Af =< f\H\B° >, and 
Af ~< f\H\B'^ >. An interesting aspect of a/ is that 
it is free of hadronic uncertainties and for the Stan- 
dard Model case it is determined fully in terms of the 
CKM parameters. This would be the case if only one 
amplitude contributes to the decay B''\B'^) /. More 
generally one has more than one diagram contributing 
with different CKM phase dependence which make the 
extraction of CKM phases less transparent. Specifically 
B^{B^) decays may in general involve penguin diagrams 
which tend to contaminate the simple analysis outlined 
above. Gronau and London have proposed an isospin 
analysis which can disentangle the effect of the tree and 
penguin contributions when the final states in B^{B^) 
are tt+tt^ and tt'^tt'^ which is useful in the analysis of 
all the CKM angles (Gronau and London, 1990; Gronau 
and London., 1991). The decay final states J/'^Ks is 
interesting in that it is a CP eigen state and it has 
a large branching ratio and to leading order is domi- 
nated by a single CKM phase. Specifically, the relation 
Aj/^!KslAj/<ijKs — 1 holds to within a percent (Boos 
et ai, 2004), A}^^^^ = sin(2/3) and A}/^^^ = 0. Thus 

B^{B^) decay into this mode gives a rather clean mea- 
surement of sin 2/3. BaBar and Belle have both measured 
CP asymmetries utilizing the charm decays. Using the 
decays B^{B^) -> J/*i^s and B"(B") ^ J/^i^L BaBar 
and Belle have obtained a determination of the CP asym- 
metry sin(2/3) and the world average for this is (Barberio 
et ai, 2006) 



sin(2/3) = 0.685 ± 0.032 



(17) 



While the analysis of CP asymmetries in the J/'i'Ks 
system is the cleanest way to determine sin(2/3) there 
are additional constraints on (3 that are indirect such as 
from Ariid, and Anis- These lead to a constraint on /3 
with P lying in the range (13°, 31°) at 95% C.L. (Charles 
et ai, 2005; Long, 2005). 

The determination of a comes from the measure- 
ment of processes of type i?° — > tt+tt^, p~^p~ since 




FIG. 1 The penquin diagram that contributes to B decays. 



the combinations of phases that enter here are via 
sin(2(/3 -I- 7)) = — sin(2Q;). One problem arises due 
to the contribution of the penguin diagram Fig.(l) 
which does not contain any weak phase. The penguin 
diagram can thus contaminate the otherwise neat weak 
phase dependence of this process. A possible cure 
come from the fact that one can use the analysis of 
(Grossman and Quinn, 1998) to put an upper limit 
on the branching ratio for i3° p'^p'^. The current 
determination of a gives a = (96 ± 13 ± 11)° (Stone, 
2006). The determination of 7 comes from the charged 
decays B^ — > D'^K^. The current experimental values 
from BaBar and Belle are 7 = (67 ± 28 ± 13 ± 11)°, and 
7 = (67t}^ ± 13 ± 11)° (Asner and Sun, 2006; Stone, 
2006). A detailed analysis of global fits to the CKM ma- 
trix can be found in (Charles, 2006; Charles et ai, 2005). 

We discuss now Z?° — 13° system. In analogy with the 
neutral B system we introduce the two neutral mass eigen 
states Di, D2 defined by 



\Di >= p\D" > +q\D" >, 
\D2 >=p\D° > -q\W> . 



(18) 



The D mesons are produced as flavor eigen states but 
they evolve as admixtures of the mass eigen states which 
govern their decays. The analysis of Z)° and Z)° decays 
by BaBar(Aubert et ai, 2007) and by Belle(Staric et ai, 
2007) finds no evidence of CP violation. For further 
details the reader is directed to (Nir, 2007b). 

The fourth piece of experimental evidence for CP vio- 
lation in nature is indirect. It arises from the existence 
of a baryon asymmetry in the universe which is generally 
expressed by the ratio 



ne/n^ = (6-l°:^.2) x 10 



-10 



(19) 



An attractive picture for the understanding of the 
baryon asymmetry is that the asymmetry was generated 
in the very early history of the universe within the 
context of an inflationary universe starting with no 
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initial baryon asymmetry (for a recent review on matter- 
antimatter asymmetry see (Dine and Kusenko, 2004)). 
The basic mechanism how this can come about was 
already enunciated a long time ago by (Sakharov, 1967). 
According to Sakharov there are three basic ingredi- 
ents that govern the generation of baryon asymmetry, 
(i) One needs a source of baryon number violating 
interactions if one starts out with a universe which 
initially has no net baryon number. Such interactions 
arise quite naturally in grand unified models and in 
string models, (ii) One needs CP violating interactions 
since otherwise would be a balance between processes 
producing particles vs processes producing anti-particle 
leading to a vanishing net baryon asymmetry. (iii) 
Finally, even with baryon number and CP violating 
interactions the production of a net baryon asymmetry 
would require a departure from thermal equilibrium. 
Thus one finds that one of the essential ingredients for 
the generation of the baryon asymmetry in the early 
universe is the existence of CP violation. However, the 
CP violation in the Standard Model is not sufiicient to 
generate the desire amount of baryon asymmetry and 
one needs a source of CP violation above and beyond 
what is present in the Standard Model. Such sources 
of CP violation are abundant in supersymmetric theories. 

In addition to the baryon asymmetry in the universe 
there are other avenues which may reveal the existence 
of new sources of CP violation beyond what exists 
in the Standard Model. The EDMs of elementary 
particles and of atoms arc prime candidates for these. 
The largest values of EDMs in the framework of the 
Standard Model SM are very small. SM predicts for 
the case of the electron the value of ~ 10~^^ecm and 
for the case of the neutron the value that ranges from 
10~^^ to 10~^^ecm (Bernreuther and Suzuki, 1991; Bigi 
and Urahsev, 1991; Booth, 1993; Gavela et al, 1982; 
Khriplovich and Zhitnitsky, 1982; Shabalin, 1983). 

So far no electric dipole moment for the electron or for 
the neutron has been detected, and thus strong bounds 

on these quantities exist. For the electron the current 
experimental limit is (Regan et at, 2002), 

\de\ < 1.6 X IQ-^^ecm (90% CL). (20) 

For the neutron the Standard Model gives dn ^ lO^''^*^ 
ecm while the current experimental limit is (Baker et al., 
2006) 

\dn\ < 2.9 X IQ-^^ecm (90% CL). (21) 

In each case one finds that the Standard Model predic- 
tion for the EDM is several orders of magnitude smaller 

than the current experimental limit and thus far beyond 
the reach of experiment even with improvement in sen- 
sitivity by one to two orders of magnitude. On the 
other hand many models of new physics beyond the Stan- 
dard Model generate much larger EDMs and such models 



are already being constrained by the EDM experiment. 
Indeed improved sensitivities in future experiment may 
lead to a detection of such effects or put even more strin- 
gent constraints on the new physics models. The EDM of 
the atoms also provides a sensitive test of CP violation. 
An example is Hg-199 for which the current limits are 
(Romalis et al., 2001), 

IdjjJ < 2 X 10-^^ecm. (22) 



IV. CP VIOLATION IN SOME NON-SUSY EXTENSIONS 
OF THE STANDARD MODEL 

While the Standard Model contains just one CP phase 

more phases can appear in extensions of the Standard 
Model. In general the violations of CP can be either ex- 
plicit or spontaneous. The CP violation is called explicit 
if redefinitions of fields cannot make all the couplings real 
in the interaction structure of the theory. The remaining 
phases provide an explicit source of CP violation. CP 
violation is called spontaneous if the model starts 
out with all the couplings being real but spontaneous 
breaking in the Higgs sector generates a non-removable 
phase in one of the vacuum expectation values in the 
Higgs fields at the minimum of the potential. Returning 
to CP violation in the extension of the Standard Model, 
such extensions could be based on an extended gauge 
group, on an extended Higgs sector, or on an extended 
fermionic content (see, for example, (Accomando et al., 
2006)). An example of a model with an extended gauge 
sector is the left-right (LR) symmetric model based 
on the gauge group SU{2)l x SU{2)r x [/(I) (Moha- 
patra and Pati, 1975). For Ug number of generations 
the number of phases is given by Nl + Nfi where 
= {ng — l)(?ig — 2)/2 is exactly what one has in 
SU{2)l X U{1)y model and Nr = ng{ng + l)/2 are 
additional set of phases that arise in the LR model. For 
the case of three generations this leads to 7 CP phases 
instead of just one CP phase that one has in the Stan- 
dard Model. An analysis of EDM in LR models for the 
electron and for the neutron is given in (Prank, 1999a,b). 

The simplest extension of the Standard Model with 
an extended Higgs sector is the so called two Higgs dou- 
blet model (Lee, 1973, 1974) (2HDM) which contains two 
SU{2) doublets which have exactly the same quantum 
numbers $i = (0^,0°), i=l,2. One problem with the 
model is that it leads to flavor changing neutral cur- 
rents (FCNC) at the tree level if one allows couplings 
of both $j to the up and down quarks. The FCNC 
can be suppressed by imposing a discrete Z2 symme- 
try (Glashow and Weinberg, 1977) such that under Z2 
one has $2 —^2 and UiR — > —UiR and the remain- 
ing fields are unaffected. Under the above symmetry 
the most general renormalizable scalar potential one can 
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write is 

Vq = - ^il^\^2 + Al($I$l)2 + A2($|$2)' + 

+A3($l$l)($^$2)^ + A4|$l$2|' + 

+{\n{^\^2f +H.C.) (23) 

However, with an exact Z2 discrete symmetry CP cannot 
be broken either exphcitly or spontaneously in a 2HDM 
model (Branco, 1980a, b; Mendez and Pomarol, 1991). 
Thus to have CP in the 2HDM model one must allow 
for violations of the discrete symmetry, but arrange for 
suppression of FCNC. If the couplings allow for FCNC 
at the tree level, then they must be suppressed cither by 
heavy Higgs masses (Branco et al., 1985; Lahanas and 
Vayonakis, 1979) or by adjustment of couplings or fine 
tunings so that FCNC are suppressed but CP violation 
is allowed (Liu and Wolfenstein, 1987). 

However, the hard breaking of the Z2 discrete sym- 
metry is generally considered not acceptable. A more 
desirable possibility is violation of the discrete symmetry 
only via soft terms (Branco and Rebelo, 1985). Here the 
FCNC are not allowed at the tree level but the inclusion 
of the soft terms allows for CP violation. Such a term is 
of the form 

Vsoft = -l4^\^2 + H.c. (24) 

Soft breaking of the Z2 symmetry can allow both explicit 
and spontaneous CP violation. Thus explicit CP viola- 
tion can occur inV — Vq + Vsofi if one has(Grzadkowski 
et al, 1999) Im{iJLl'^X^) ^ 0. For the case when 
Im{iJ,l'^ X5) = a spontaneous violation of CP can arise. 
Specifically, in this case one can choose phases so that 
< ^-^ >= V1/V2 {vi > 0) and < $2 >= e*^t;2/v^ 
{v2 > 0) with the normalization 

^Jvl+vl = 2mwl92 = 246GeV. (25) 

The conditions for CP violation in a 2HDM model, 
both explicit and spontaneous, have more recently been 
studied using basis independent potentially complex 
invariants which are combinations of mass and coupling 
parameters. These invariants also are helpful in distin- 
guishing between explicit and spontaneous CP violation 
in the Higgs sector. For further discussion, the reader is 
refereed to the works of (Botella and Silva, 1995; Branco 
et al., 2005; Davidson and Haber, 2005; Ginzburg and 
Krawczyk, 2005; Gunion and Haber, 2005; Lavoura 
and Silva, 1994). While the spontaneous breaking of 
CP discussed above involves SU{2) Higgs doublets 
which may enter in the spontaneous breaking of the 
electro-weak symmetry, similar spontaneous violations 
of CP can occur in sectors not related to electro-weak 
symmetry breaking. 

In the absence of CP violation, the Higgs sector of the 
theory after spontaneous breaking of the SU (2)^ x [/(1)y 



symmetry gives two CP even, and one CP odd Higgs in 
the neutral sector. In the presence of CP violation, either 
explicit or spontaneous, the CP eigcnstates mix and the 
mass eigenstates are admixtures of CP even and CP odd 
states. The above leads to interesting phenomenology 
which is discussed in detail in (Grzadkowski et al., 1999; 
Mendez and Pomarol, 1991). The number of indepen- 
dent CP phases increases very rapidly with increasing 
number of Higgs doublets. Thus, suppose we consider 
an n/3 number of Higgs doublets. In this case the num- 
ber of independent CP phases that can appear in the 
unconstrained Higgs potential is (Branco et al., 2005) 
Np — n|)(n|) — l)/4 — (nu — 1). For no — 1, 2, 3 one gets 
Np — 0, 2, 16, and thus the number of independent CP 
phases rises rather rapidly as the number of Higgs dou- 
blets increases. An analysis of the EDMs in the two Higgs 
model is given in (Barger et al, 1997; Hayashi et al., 
1994). Finally, one may consider extending the fermionic 
sector of theory with inclusion of additional generations. 
Such an extension brings in more possible sources of CP 
violation. Thus, for example, with four generation of 
quarks the extended CKM matrix will be 4 x 4. Such a 
matrix can be parameterized in terms of six angles and 
three phases (Barger et al, 1981; Oakes, 1982). Thus 
generically extensions of the Standard Model will in gen- 
eral have more sources of CP violation than the Stan- 
dard Model . We discuss CP violation in supersymmet- 
ric theories next. While the spontaneous breaking of 
CP discussed above involves SU (2) Higgs doublets which 
may enter in the spontaneous breaking of the electro- 
weak symmetry, similar spontaneous violations of CP can 
occur in sectors not related to electro-weak symmetry 
breaking. 



V. CP VIOLATION IN SUPERSYMMETRIC THEORIES 

Supersymmetric models are one of the leading can- 
didates for new physics (for review see (Haber and 
Kane, 1985; Martin, 1997; Nath et al; Nilles, 1984)) 
since they allow for a technically natural solution to 
the gauge hierarchy problem. However, supersymmtetry 
is not an exact symmetry of nature. Thus one must 
allow for breaking of supersymmetry in a way that does 
not violate the ultraviolet behavior of the theory and 
destabilize the hierarchy. This can be accomplished 
by the introduction of soft breaking. However, the 
soft breaking sector in the minimal supersymmetric 
standard model (MSSM) allows for a large number of 
arbitrary parameters (Dimopoulos and Georgi, 1981; 
Girardello and Grisaru, 1982). Indeed in softly broken 
supersymmetry with the particle content of MSSM 
additionally 21 masses, 36 mixing angles and 40 phases 
(Dimopoulos and Sutter, 1995). which makes the model 
rather unpredictive. 

The number of parameters is significantly reduced 
in the minimal supergravity unified models under the 
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assumptions of a flat Kahlcr metric as explained below. 
The minimal supcrgravity model and supergravity model 
in general are constructed using techniques of applied 
N=l supergravity, where one couples chiral matter mul- 
tiplets and a vector multiplet belonging to the adjoint 
representation of a gauge group to each other and to 
supergravity. The supergravity couplings can then be 
described in terms of three arbitrary functions: the 
superpotential W{zi) which is a holomorphic function 
of the chiral fields Zi, the Kahler potential K{zi,zl) 
and the gauge kinetic energy function fa/di^iizj) which 
transforms like the symmetric product of two adjoint 
representations. In supergravity models supersymmetry 
is broken in a so called hidden sector and is commu- 
nicated to the physical sector where quarks and Icpton 
live via gravitational interactions. The size of the soft 
breaking mass, typically the gravitino mass ma, is 

^ K^l < Wh > [, where Wh is the superpotential in 
the hidden sector where supersymmetry breaks and 
K = 1/Mpi, where Mpi is the Planck mass. The simplest 
model where supersymmetry breaks in the hidden sector 
via a super Higgs effect is given by Wh = m^z where 
z is the Standard Model singlet super Higgs field. The 
breaking of supersymmetry by supergravity interactions 
in the hidden sector gives z a VEV of size ~ k~^, and 
thus with TO ^ 10^'^"^^ GeV, the soft breaking mass is 
of size ~ 103 GeV. 

In the minimal supergravity model one assumes that 
the Kahler potential has no generational dependence and 
is flat and further that the gauge kinetic energy function 

is diagonal and has no field dependence, i.e., one has 
effectively fa^ ~ Sa/s- In this case one finds that the 
low energy theory obtained after integrating the GUT 
scale masses has the following soft breaking potential 
(Chamseddine et at, 1982; Hall et al., 1983; Nath et al., 
1983) 

VsB = mi A«A« + mlzazl + {AqW^^^ + BqW^^^ + H.c.) 

(26) 

where W^"^^ is the quadratic and W^^^ is cubic in the 
fields. 

The physical sector of supergravity models consist of 

the MSSM fields, which include the three generations of 
quarks and leptons and their superpartners, and a pair 
of SU{2)l Higgs doublets Hi and H2 and their super- 
partners which are the corresponding Higgsino fields Hi 
and if 2- For the case of MSSM one has 

W^^^ = fj.aHiH2, 
W^^^ = QYuH2U^ + QYoHid" + LYEH2e^ (27) 

Here Hi is Higgs doublet that gives mass to the bottom 
quark and the lepton, and H2 gives mass to the up 

quark. As is evident from Eqs(26) and (27) the minimal 
supergravity theory is characterized by the parameters 



: 1710,111,1 , Aq, Bq and /iq- An interesting aspect of 
supergravity models is that they allow for sponta- 
neous breaking of the SU{2)l x U{1)y electroweak 
symmetry (Chamseddine et ai, 1982). This can be ac- 
complished in an efficient manner by radiative breaking 
using renormalization group effects (Alvarez-Gaume 
et al., 1983; EUis et al, 1983; Ibanez and Lopez, 1984; 
Ibanez et al, 1985; Ibanez and Ross, 1982, 2007; Inoue 
et al, 1982). 

To exhibit spontaneous breaking one considers the 

scalar potential of the Higgs fields by evolving the po- 
tential to low energies by renormalization group effects 
such that 

V = Vo + AV (28) 

where Vq is the tree level potential (Haber and Kane, 
1985; Nath et al; Nilles, 1984) 

Vo = ml\Hif + ml\H2f + {mlHi.H2 + H.c.) 

+ ^l±^\Hi\^ + - §\Hi.H2\' 

+ ^l^\Hi\^\H2\^ (29) 
8 

and Ay is the one loop correction to the effective poten- 
tial and is given by (Arnowitt and Nath, 1992; Carena 
et al., 2000; Coleman and Weinberg, 1973; Weinberg, 
1973) 

= ^.StriM\Hi,H2)(Xo,'J^^^- |)). 

(30) 

Here Str = T,iCi{2,Ji + 1)(— 1)^'^% where the sum runs 
over all particles with spin Ji and Ci{2Ji + 1) counts 
the degrees of freedom of the particle i and Q is the 
running scale which is to be in the electroweak region. 
The gauge coupling constants and the soft parameters 
are subject to the supcrgravity boundary conditions: 
0:2(0) = aa = §ar(0); TOf(O) = ml + i = 1,2; 
and toKO) = Bq^iq- As one evolves the potential down- 
wards from the GUT scale using renormalization group 
equations(Jack et al., 1994; Machacek and Vaughn, 1983, 
1984, 1985; Martin and Vaughn, 1994), a breaking of 
the electro-weak symmetry occurs when the determi- 
nant of the Higgs mass^ matrix turns negative so that 
(i) m\m\ — 2to| < 0, and further for a stable minimum 
to exist one requires that the potential be bounded from 
below so that (m) m\ + m\ — 2|to§| > 0. Additionally 
one must impose the constraint that there be color and 
charge conservation. Defining Vi =< Hi > as the VEV 
of the neutral component of the Higgs Hi, the neces- 
sary conditions for the minimization of the potential, i.e., 
dV/dvi = , gives two constraints. One of these can be 
used to determine the magnitude |/Lto| and the other can 
be used to replace Bq by tan/3 =< H2 > / < Hi >. In 
this case the low energy supergravity model or mSUGRA 
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can be parameterized by toq ,mi,AQ, tan f3 and sign(/io ) • 
It should be noted that fixing the value using radia- 
tive breaking does entail fine tuning but a measure of this 
is model dependent (see, for example, (Chan:1997bi) and 
the references therein). The above discussion is for the 
case when there arc no CP violating phases in the theory. 
In the presence of CP phases mi , Aq, /io become complex 
and one may parameterize them so that 

TOi=|mi|e'*i, Ao = \Ao\e"'°,no = \fio\e''-o. (31) 

Now not all the phases are independent. Indeed, in this 
case only two phase combinations are independent, and 
in the analysis of the EDMs one finds these to be + 
6fj,Q and ao + . Often one rotates away the phase of 
the gauginos which is equivalent to setting = 0, and 
thus one typical choice of parameters for the complex 
mSUGRA '(cmSUGRA) case is 

mo, |mi|, tan/3, \Ao\; ao, 0^^^ (cmSUGRA).(32) 

However, other choices are equally valid: thus, for 
example, the independent soft breaking parameters can 
be chosen to be mo, |mi|, tan/3, |^o|) oiOi 
mSUGRA model was derived using a super Higgs effect 
which breaks supersymmctry in the hidden sector by 
VEV formation of a scalar super Higgs field. Alternately 
one can view breaking of supersymmctry as arising 
from gaugino condensation where in analogy with QCD 
where one forms the condensate qq one has that the 
strong dynamics of an asymptotically free gauge theory 
in the hidden sector produces a gaugino condensate 
with < A7''A >= A^. The above can lead typically to 
supersymmctry breaking and a gaugino mass of size 
ma ~ K^k^. With |A > I ~ (1012-13) GeV one will have 
an ma again in the electro- weak region (Dine et al., 
1985; Ferrara et al, 1983; Nilles, 1982; Taylor, 1990). 

The assumption of a flat Kahler potential and of a fiat 
kinetic energy function in supergravity unified models is 
essentially a simplification, and in general the nature of 
the physics at the Planck scale is largely unknown. For 
this reason one must also consider more general Kahler 
potentials (Kaplunovsky and Louis, 1993; Soni and Wel- 
don, 1983) and also allow for the non-universality of the 
gauge kinetic energy function. In this case the number 
of soft parameters grows, as also do the number of CP 
phases. Thus, for example, the gaugino masses will be 
complex and non-universal, and the trilinear parameter 
Aq, which is in general a matrix in the generation space, 
will also be in general non-diagonal and complex. A 
simplicity assumption to maintain the appropriate con- 
straints on flavor changing neutral currents is to assume a 
diagonal form for Aq at the GUT scale. Additionally, the 
Higgs masses for Hi and H2 at the GUT scale could also 
be non-universal. Thus in general for the non-universal 
supergravity unification a canonical set of soft parame- 
ters at the GUT scale will consist of (Matalliotakis and 



Nilles, 1995; Nath and Arnowitt, 1997; Olechowski and 
Pokorski, 1995; Polonsky and Pomarol, 1995) 

rriHi = mo(l + Si), i = 1,2 
nia = Imale^^" , a = 1,2, 3 
Aa = |Aa|e^"% a = l,2,3 (33) 

which contain several additional CP phases beyond 
the two phases in complex mSUGRA. However, not 
all the phases are independent, as some phases can be 
eliminated by field redefinitions. Indeed in physical 
computations only a certain set of phases appear, as 
discussed in detail in (Ibrahim and Nath, 2000c) (also 
see Appendix XVI. E). It should be kept in mind that 
for the case of non-universalities the renormalization 
group evolution gives an additional correction term at 
low energies (Martin and Vaughn, 1994). 

As is apparent from the preceding discussion radiative 
breaking of the electroweak symmetry plays a central role 
in the supergravity unified models. An interesting phe- 
nomena here is the existence of two branches of radiative 
breaking: one is the conventional branch known since the 
early eighties (we call this the ellipsoidal branch (EB)) 
and the other was more recently discovered, i..e, it is 
the so called hyperbolic branch (HB). The two branches 
can be understood simply by examining the condition of 
radiative breaking which is a constraint on the soft pa- 
rameters mo,m'j^2>^o of the form (Chan et al, 1998) 

Cim^ + C3m[% + C'^Al + A/xLp = ^ + m'- (34) 

Here A/i^^^^ is the loop correction (Arnowitt and Nath, 
1992; Carena et al, 2000), and m[^^ = mi/a + iAoQ/Cs, 
where Ci are determined purely in terms of the gauge 
and the Yukawa couplings but depend on the renor- 
malization group scale Q. The behavior of radiative 
breaking is controlled in a significant way by the loop 
correction Afif^^^ especially for moderate to large values 
of tan/3. For small values of tan/3 the loop correction 
A/i^ is small around Q ~ Mz, and the C, are positive 
and thus Eq.(34) is an ellipsoidal constraint on the soft 
parameters. For a given value of /U, Eq.(34) then puts 
an upper limit on the sparticle masses. However, for 
moderate to large values of tan /3, A/i^ becomes sizable. 
Additionally Ci develop a significant Q dependence. It 
is then possible to choose a point Q = Qq where A/i^ 
vanishes and quite interestingly here one finds that one 
of the Ci (specifically Ci) turns negative, drastically 
changing the nature of the symmetry breaking con- 
straint Eq.(34) on the soft parameters. Thus in this 
case the soft parameters in Eq.(34) lie on the surface of 
a hyperboloid and thus for a fixed value of /i the soft 
parameters can get very large with mo getting as large 
as 10 TeV or larger. The direct observation of squarks 
and sleptons may be difficult on this branch, although 
charginos, neutralinos and even gluino may be accessible. 
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However, the HB does have other desirable features such 
as suppression of flavor changing neutral currents, and 
suppression of the SUSY EDM contributions. Further, 
HB still allows for satisfaction of relic density constraints 
with R parity conservation if the lightest neutralino is 
the lightest supersymmetric particle (LSP). We note in 
passing that the so called focus point region (Feng et al., 
2000) is included in the hyperbolic branch (Baer et al., 
2004; Chan et al., 1998; Lahanas et al, 2003). 

There is a potential danger in supergravity theories 
in that the hierarchy could be destabilized by non- 
renormalizable couplings in supergravity models since 
they can lead to power law divergences. This issue 
has been investigated by several authors: at one loop 
by (Bagger and Poppitz, 1993; Gaillard, 1995) and at two 
loop by (Bagger et al., 1995). The analysis shows that at 
the one loop level the minimal supersymmetric standard 
model appears to be safe from divergences (Bagger and 
Poppitz, 1993). In addition to the breaking of supersym- 
metry by gravitational interactions, there are a variety 
of other scenarios for supersymmetry breaking. These 
include gauge mediated and anomaly mediated breaking 
for which reviews can be found in (Giudicc and Rattazzi, 
1999; Luty, 2005). Finally as is clear from the preceding 
discussion in supergravity models and in MSSM there is 
no CP violation at the tree level in the Higgs sector of 
the theory. However, this situation changes when one in- 
cludes the loop correction to the Higgs potential. This 
leads to the generation of CP violating phase for one the 
Higgs VEVs and leads to mixings between the CP even 
and the CP odd Higgs fields. This phenomenon is very 
interesting from the experimental view point and will be 
discussed in greater detail later. 

While the Standard Model contribution to the EDMs 
of the electron and of the neutron is small and beyond 
the pale of observation of the current or the future experi- 
ment, the situation in supersymmetric models is quite the 
opposite. Here the new sources of CP violation can gen- 
erate large contributions to the EDMs even significantly 
above the current experimental limits. Here one needs 
special mechanisms to suppress the EDMs such as mass 
suppression (Kizukuri and Oshimo, 1992; Nath, 1991) or 
the cancelation mechanism to control the effect of large 
CP phases on the EDMs. (Chattopadhyay et al, 2001; 
Ibrahim, 2001b; Ibrahim and Nath, 1998a,b,c, 2000d). 
Specifically for the cancelation mechanism the phases can 
be large and thus affect a variety of CP phenomena which 
can be observed in low energy experiments and at accel- 
erators. The literature on this topic is quite large. A 
sample of these analyses can be found in (Akeroyd and 
Arhrib, 2001; Alan et at, 2007; Bartl et al, 2006, 2004c; 
Boz, 2002; Chattopadhyay et al, 1999; Demir, 1999; 
Falk and Olive, 1998; Gomez et al, 2004a,b, 2005, 2006; 
Huang and Liao, 2000a,b, 2002; Ibrahim et al, 2001; 
Ibrahim and Nath, 2000a,b,c, 2001a,b, 2002, 2003a,b,c, 
2004, 2005; Ibrahim et al, 2004). 



VI. CP VIOLATION IN EXTRA DIMENSION MODELS 

Recently there has been significant activity in the 
physics of extra dimensions (Antoniadis, 1990; Anto- 
niadis et al, 1998; Arkani-Hamed et al, 1998; Gog- 
berashvili, 2002; Randall and Sundrum, 1999a, b). One 
might speculate on the possibility of generating CP vio- 
lation in a natural way from models derived from extra 
dimensions (For an early work see (Thirring, 1972)). It 
turns out that it is indeed possible to do so (Branco et al, 
2001; Burdman, 2004; Chaichian and Kobakhidze, 2001; 
Chang et al, 2001; Chang and Mohapatra, 2001; Dienes 
et al; Grzadkowski and Wudka, 2004; Huang et al, 2002; 
Klilebnikov and Shaposhnikov, 1988; Sakamura, 1999). 
The idea is to utilize properties of the hidden compact 
dimensions in extra dimension models. Thus in extra di- 
mension models after compactification the physical four 
dimensional space is a slice of the higher dimensiional 
space and such a slice can be placed in different loca- 
tions in extra dimensions. In the discussion below we 
will label such a slice as a brane. We consider now a 
simple argument which illustrates how CP violation in 
extra dimension models can arise (Chang and Mohapa- 
tra, 2001). Thus consider a.U{l) gauge theory with left- 
handed fermions ^'^ (i=l-4), where i = 1, 2 have charges 
-hi and i = 3,4 have charges —1, and also consider a 
real scalar field $ which is neutral. We assume that 
the fermion fields are in the bulk and the scalar field 
is confined to the y = brane. The fields ^'i,\I/2 and 
i> are assumed to be even and ^zl , are assumed 
to be odd under y —y transformation. Further, un- 
der CP symmetry define the fields to transform so that 
*iL ^ (*3^,)^*2L ^ (*4L)^ and $ ^ where 
{^lY has the meaning of a 4D charge conjugate of ^. 
One constructs a 5D Lagrangian invariant under y — * —y 
transformation of the form 

M5-lA5^(t/)$[*r^C-l*2L - (*3l)^^C-1(*4l)1 

+M[*rL<^"'*2L - (*3l)^^C-^(*4l)T + H.c. (35) 

On integration over the y co-ordinate the interaction 
terms in 4D arise from the couplings on the y=0 brane 
and thus the zero modes of the fields odd in y are ab- 
sent, which means that the effective interaction at low 
energy in (A^ + m)^il^2L which violates CP provided 
Im{X* ji) ^ 0. Next we discuss a more detailed illustra- 
tion of this CP violation arising from extra dimensions. 
This illustration is an explicit exhibition of how viola- 
tions of CP invariance can occur in the compactification 
of a 5D QED (Grzadkowski and Wudka, 2004). Thus 
consider the Lagrangian in 5D of the form 

A = -Jv^Mjv + *(n'^-DM-rni)* + £gh. (36) 

Here Vm is the vector potential in 5d space with co- 
ordinates z^, where M = 0,1,2,3,5 so that z'^ = 
{x'^,y), where /i = 0,1,2.3, and where Dm = 9m + 
igzoYM is the gauge covariant derivative, with 55 the 
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U{1) gauge coupling constant, and q the charge of 
fermion field. The theory is invariant under the following 
gauge transformations 

VMiz)^VMiz)+dMX{z), (37) 
and additionally under the CP transformations in 5D 

^m^^M^M^ yM^^MyM^ ^ ^ (38) 

where r]^'"^'^ = —I = —rp'^ and P = 1. We compactify 
the theory in the fifth dimension on a circle with radius 
R assuming periodic boundary conditions for the gauge 
fields but assuming the twisted boundary condition for 
the fermion field 

i){x,y + R) = e"'i^{x,y). (39) 

One can now carry out a mode expansion in 4D and re- 
covers a masslcss zero mode V^{x) for the vector field 
(the photon). One also finds in addition a massless field 
(f){x) which is the zero mode of the V5{x,y) expansion. 
This is so because while V^,n ^ modes can be elimi- 
nated by an appropriate gauge choice, while is a gauge 
singlet and remains in the spectrum. We note in pass- 
ing that the presence of the zero mode is a consequence 
of the specific compactification chosen. Thus compacti- 
fication, on 5^/Z2, rather than on the circle will remove 
the field 0. Now while </> is massless at the tree level, 
it can develop a mass when loops contributions are in- 
cluded. Thus an analysis of one loop effective potential 
gives (Grzadkowski and Wudka, 2004). 

Veff = Yl^p^^M + WiLuili) + ^LMm 

i 

where j3i = mR, = exp{iujiR — (3i), and where 
= {ai + gsQiRipo), and (f)o —< (j) >) ^ind L,^ is the 
polylogarithm function. 

Now it turns out that for the case when one has a 
single fermion, there is no CP violation, but CP vio- 
lation is possible when there are two fermions and one 
can assume the boundary conditions in this case so that 
V'i(a;,y + R) = ipi{x,y) and ip2{x,y + R) = e^'4}2{x,y). 
In this situation the Yukawa couplings for the fermions 
violate CP. An interesting phenomenon here is that the 
above mechanism exhibits examples of both spontaneous 
CP violation as well as explicit CP violation. Thus for the 
case a = 0, TT one finds that the effective potential is sym- 
metric in and one has two degenerate minima away 
from 00 = and thus here one has spontaneous breaking 
of CP. For other choices of a, the effective potential is 
not symmetric in and one has explicit violation of CP. 
The fact that CP is indeed violated in this example can 
be tested by an explicit computation of the EDM of the 
fermions which is non-vanishing and suppressed by the 
inverse size of the extra dimension. 



V„(l<r'TeV') V^(l(r'TeV') 




FIG. 2 An exhibition of the phenomena of spontaneous vs ex- 
plicit breaking in a 5D compactification model (Grzadkowski 
and Wudka, 2004). The figure gives the efi^ectivc potential 
Veff for four cases of twist angles with a = 0, 7r/2, tt, Sn/2. 
The cases a = 0, tt correspond to spontaneous breaking and 
a = 7r/2,37r/2 correspond to explicit breaking. 

We turn now to another mechanism for the generation 
of CP violation in extra dimensional theories. This sce- 
nario is that of split fermions where the hierarchies of 
fermion masses and couplings are proposed to arise from 
a fermion location mechanism under a kink background 
wherein the quark and leptons of different generations 
being confined to different points in a fat brane (Arkani- 
Hamed and Schmaltz, 2000; Kaplan and Tait, 2000, 2001; 
Mirabelli and Schmaltz, 2000). To illustrate the fat brane 
paradigm consider the 4-|-l dimensional action of two 
fermions 

55 = J d^xdylQl^Md^ + ^Q{y)]Q + 

+U[i^Md^ + ^u{y)]U + kHQ'^U]. (41) 

The quantities ^q.u are potentials which confine the 
quarks at different points in the extra dimension. As 
a model one may consider these as Gaussian functions 
centered around points Iq ( i.e., functions of the form 
exp{—fi^{y — Iq)'^)) and Z„ where 1/2,^/7 is the width of 
the Gaussian. After expanding the fields in their nor- 
mal modes and integrating over the extra dimension the 
Yukawa interaction in 4D including the generation index 
will take the form 

jCy = X-jQiUjH + \%QiDjH*, (42) 

where is defined by 

A^, = Kye-5'''('''*-'"i\ (43) 

and \fj is similarly defined. The above structure 
indicates that the Yukawa textures are governed by the 

location of the quarks in the extra dimension. Detailed 
analyses, however, indicate that this scenario leads to an 
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insufficient amount of CP violation to explain the value 
of eK in Kaon decay. Thus the scenario above gives a 
value of the Jarskog invariant J < 5 x 10^^ while one 
needs J ~ 10~^ to get the proper value of ex- The above 
shortcoming can be corrected by extending the analysis 
to two extra dimensions (Branco et al., 2001). In this 
case one finds the Jarlskog invariant J ~ 2.2 x 10~^ 
which is of desired strength to explain CP violation in 
the Kaon decay. An extension to include masses for the 
charged leptons and neutrinos has been carried out in 
(Barenboim et al, 2001). 

An analysis using the fermion localization mechanism 
for generating quark-lepton textures within a supersym- 
metric SU(5) GUT theory is carried out in the analy- 
sis of (Kakizaki and Yamaguchi, 2004) where the differ- 
ent SU{5) chiral multiplets are localized along different 
points in the extra dimension. The analysis allows one 
to generate a realistic pattern of quark masses and mix- 
ings and lepton masses. The CP violation is of sufficient 
strength here since J ~ O(10~^). An additional feature 
of this model is that dimension 5 proton decay operators 
are also naturally suppressed due to the fact that these 
operators contain an overlap of wavefunctions of different 
chiral multiplets and are thus exponentially suppressed. 

Similar analyses can be carried out in the framework 
of a non-factorizable geometry (Abe et al, 2001; Chang 
et at, 2000; Grossman and Neubert, 2000; Huber and 
Shafi, 2001) based on the metric 

ds"^ = e-^''^y\dxf -dy'^, (44) 

where cr(y) = k\y\. Under the Z2 orbifold symmetry the 
5D fermion transform as *(— 2/)± = ±75\E'(j/)±. The 
have the mode expansion 

^{x,y)± = -^f;Vn±(2;)/i"^(y). (45) 

The zero modes of 'i>± are the left-handed and the right 
handed Weyl spinors. Masses for these are generated by 
the 5D Higgs couplings which are of the form 

J d^xdy^XijH^i+^j-. (46) 

For the zero mode they give rise to a Dirac mass term of 
the form (Huber and Shafi, 2001) 

dyX,,H{y)fil\y)f^°\y) (47) 

where 

/(o) = ( ; ' - V ^e(^-)- (48) 

where c is a parameter that characterizes the location 
of the fermion in the extra dimension. For c < 1/2 the 



fermion is localized near the y = brane while for r = ttTc 
it is localized near y = -kTc brane. With the appropri- 
ate choice of the c's one may generate a realistic pattern 
of quark masses and mixings and a realistic CKM ma- 
trix. However, an explicit determination of the Jarlskog 
invariant appears not to have been carried out. The tex- 
ture models using extra dimensions do generally require 
a high level of fine tuning in the selection of locations 
where the fermions are placed. Thus models of this type 
do not appear very natural. For related works on CP 
violation and extra dimensions see (Dooling et al., 2002; 
Huang et al, 2002; Ichinose, 2002; Sakamura, 1999). 



VII. CP VIOLATION IN STRINGS 

We discuss now the possible origins of CP violation in 
SUSY, string and brane models (for review of string the- 
ory see (Green et al, 1987a, b; Polchinski, 1998a, b)). One 
possible origin is string compactification(Bailin et al, 
1998a,b, 2000; Dent, 2001, 2002; Faraggi and Vives, 2002; 
Kobayashi and Lim, 1995; Witten, 1985; Wu et al, 1991). 
One may call this hard CP violation since this type of 
CP violations can exist even without soft terms. Now 
Yukawa couplings which are formed via string compacti- 
fication will carry this type of CP violation and the CKM 
phase 5cKM which arises from the Yukawas is therefore 
a probe of CP violation arising from string compactifi- 
cation (assuming there is no CP violation arising from 
the Higgs sector) . A second source of CP violation is via 
soft breaking. If SUSY contributions to K and B physics 
turn out to be small, then one has a plausible bifurca- 
tion, i.e., the CP violations in K and B physics are probe 
of string compactification, and baryogenesis and other 
CP phenomena that may be seen in sparticle decays etc 
become a probe of soft breaking. 

Regarding soft breaking in string theory, such an anal- 
ysis would entail specifying the Kahler potential, the su- 
perpotential, and the gauge kinetic energy function on 
the one hand and the mechanism of breaking on the 
other. Each of these are model dependent. However, 
it is possible to parameterize the breaking as in gravity 
mediated breaking in supergravity. Thus one can write 
the general form of the soft terms in the form 

+ ^{Ba0tta0CaC0 + H.C. )+•••, (49) 

where the general expressions for the scalar masses rUa, 
trilinear couplings A^fj-y and the bilinear term B can be 
given. For the case when K^j^ = d^^K^, one has (Brig- 
nole et al, 1994; Kaplunovsky and Louis, 1993) 

ml = ml/2 + - F'F'^didj\n{K^) , 
Ac,(i-, = cF' (diK + dj MY^py) - di MK^KpK^)) , 
B^p = cF' {diK + di ln(/xa^) - di HKo^Kf})) + ■■■ , (50) 
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while the gaugino masses are given by 

An efficient way to parameterize is given by (Brignole 
et al, 1994) 

= V3m|(5 + 5*)sin6'e-^^^ 
F' = V^m, 3 (T + T*) cos 9&^e-''^% (52) 

where 9, &i parameterize the Goldstino direction in the 
S, Ti field space and 7s and 7j are the F^ and F^ phases, 

and el + el + el = i. 

A. Complex Yukawa couplings in string compactifications 

The Yiikawa coiiplings arise at the point of string com- 
pactification, and it is interesting to ask how the Yukawa 
couplings develop CP phases. It is also interesting to 
determine if such phases are small or large. Consider, 
for example, the compactification of the Eg x Es het- 
erotic string on a six dimensional Calabi-Yau (CY) man- 
ifold. In this case the massless families are either (1,1) or 
(2,1) harmonic forms. For the case when hodge number 
hn > h2i, the massless mirror families are (1,1) forms 
while if /i2i > hn the massless families are (2,1) forms. 
For the case when the families are (1,1) the cubic cou- 
plings among the families have been discussed in (Stro- 
minger, 1985). The analysis for the case when /121 > hn 
is more involved. One specific model of interest that can 
lead to complex Yukawas corresponds to compactification 
on the manifold ifo(Gepner, 1988; Schimmrigk, 1987) 

3 

= ^zf + ao{ziZ2Z3) =0 

i=0 

P2^^^ix3 = (53) 

i=0 

which is deformed from the manifold Kq (corresponding 

to the case ao = 0) in the ambient space CP^ xCP^ by a 
single (2,1) form (21^2-23) ■ The has 35 /i2i forms and 
8 hn forms, giving an Euler characteric x = 2(/i2i — ^11) 
and the number of net mass less families is |x|/2 (Sotkov 
and Stanishkov, 1988). 

By modding out by two discrete groups and Z'-^ one 
gets a three generation model. The discrete symmetries 
are Z^ and Z'^ where 

Zz: g : {zq, zi, Z2, Z3 : xi,X2,X3) ^ 

{Z0,Z2,Z3,Zi;X2,X3,X3,xl), 

Z3 : h: {zo,Zi,Z2,Z3 : Xi,X2,X3) ^ 

{zo,zi,Z2,Z3;xi,ax2,a'^X3). (54) 

where = 1 , a 7^ 1 . The group is not freely acting 
and leaves three tori invariant. These invariant tori have 



to be blown up in order to obtain a smooth CY manifold. 
Such a blowing up procedure produces six additional 
(2,1) and (1,1) forms which, however, leave the net num- 
ber of generations unchanged. One considers now the flux 
breaking of E^ on this manifold. If one embeds a single 
factor, Z3 or in the Eq, then Eq can break to SU{3)^ 
or SU{6) X U (1) each of which leave the Standard Model 
gauge group unbroken. However, the case SU{6) x U{1) 
cannot be easily broken further since an adjoint repre- 
sentation does not arise in the massless spectrum. Thus 
typically one considers the SU{3)^ possibility. In this 
case there arc two possibilities : Case(A), where Z3 is 
embedded trivially and Z3 is embedded non-trivially, and 
case (B) where Z'^ is embedded trivially and Z3 is em- 
bedded non-trivially. Now for case (A) one may choose 
Ug = {id)c X {id)L X {id)R, Uh = {id)c x a{id)L x a{id)R, 
where Ug is defined so that g ^ Ug \s & homomorphism 
of Z3 into Eq 9 /7g(Witten, 1985), and similarly for Uh, 
where (id) stands for an identity matrix, and Cl, R stand 
for color, left and right -handed subgroups of SU{?>Y. 
The analysis of Yukawa couplings in this case has been 
carried out and the couplings can be made all real. Thus 
in this case there is no CP violation arising in the Yukawa 
sector at the compactification scale. 

We consider next case (B) where essentially one has an 
interchange in the definitions of Ug and Uf so that 

Ug = {id)c X a{id)L x a{id)ji, 

Uh = {id)c X {id)L X {id)R (55) 

In this case the massless states that survive flux breaking 
of Eq transform under Z3 as follows 

Z3L = L, Z3Q = aQ, ZsQ^ = OL^Q" 
Z3L = L, Z3Q = o?Q, = olQ" (56) 

where the leptons transform as L(l,3,3), quarks 
as (3(3,3,1), and conjugate quarks as (5'^(3, 1,3). The 
barred quantities represent the mirrors, so that L(l, 3, 3), 
(5(3,3,1), and (5'^(3, 1,3). In this model the number of 
generations and mirror generation are identical to that 
of the Tian-Yau model(Greene et al, 1986, 1987) so that 
there are 9 lepton generations and 6 mirror generations, 
7 quark generations and 4 mirror quark generations, 
7 conjugate quark generations and 4 mirror conjugate 
quark generations, providing us with three net families 
of quarks and leptons. The analysis of Yukawa couplings 
has been carried out on the manifold Kq by many author. 

Our focus here is the (27)^ couplings which are unaf- 
fected by the instantons (Distler and Greene, 1988) and 
here one can use the techniques of (Candelas, 1988) to 
determine the couplings. An analysis for case (B) was 
carried out in (Wu et ai, 1991). The Yukawa couplings 
determined in this fashion have unknown normalizations 
for the kinetic energy. However, symmetries can be used 
to obtain constraints on the normalizations. Including 
these normalization constraints into account it is found 
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that Yukawas depend on a in a non-trivial manner, and 
thus CP is violated in an intrinsie manner. Further, the 
CP phase entering in the coupling is large. The CP viola- 
tion on the Kq manifold persists even when the modulus 
ao is real, so in this sense CP violation is intrinsic. 

B. CP violation in orbifold models 

Next we discuss the possibility of spontaneous CP 
violation in some heterotic string models. What we 
consider are field point limits of such models so we 
arc essentially discussing supergravity models with the 
added constraint of modular invariance (T duality) . The 
duality constraints have been utilized quite extensively 
in the analysis of gaugino condensation and SUSY 
breaking (Binetruy and Gaillard, 1991; Cvetic et ai, 
1991; Ferrara et al, 1990; Font et ai, 1990; Gaillard and 
Nelson, 2007; Nilles and Olechowski, 1990) and have 
also been utilized recently in the analysis of spontaneous 
breaking of CP (Acharya et ai, 1995; Bailin et al., 1997; 
Dent, 2001, 2002; Giedt, 2002). 

The scalar potential in supergravity and string theory 
is given by (Chamseddine et al., 1982; Cremmer et al., 
1982) 

V = e^[{K-yjDiWDtw^ - 3WW^ + Vd, (57) 

where K is the Kahlcr potential, W is superpotential and 
DiW = Wi + KiW , with the subscripts denoting deriva- 
tives with respect to the corresponding fields. As noted 
above we now use the added constraint of T-duality sym- 
metry. Specifically we assume that the scalar potential in 
the effective four dimensional theory depends on the dila- 
ton field S and on the (Kahler) moduli fields Ti (i=l,2,3), 
and it is invariant under the modular transformations (to 
keep matters simple, we do not include here the depen- 
dence on the so called complex structure ?7-moduli) 

Ti^Tl^^^^^—^,{aidi-biC,) = l, (58) 
tCiii + di 

where ai.bi, Ci, di G Z. Under the modular transforma- 
tions, K and W undergo a Kahler transformation while 
the scalar potential V is invariant. For the Kahler poten- 
tial we assume essentially a no scale form (Lahanas and 
Nanopoulos, 1987) 

K = D{z) - ^ log{Ti + Ti) + KijQ\Qj + HjjQiQj, 

i 

where D{z) = —log{z), and for z one may consider 
1 ^ 

z={S + S+^Y.^f^^og{Ti + fi)), (59) 

i 

where is the one loop correction to the Kahler poten- 
tial from the Greene-Schwarz mechanism, and Q are the 



matter fields consisting of the quarks, the leptons and 
the Higgs. For the superpotential in the visible sector 
one may consider 

W„ = fiijQiQj + XijkQiQjQk. (60) 

Under T-duality, Q's transform as 

Qi^Qini{iciTi + dif'3i. (61) 

In general, Kjj,Hjj,iJ,jj and Xjjk are functions of the 

moduli. The constraints on r^^, are such that V is mod- 

Qi 

ular invariant. Analyses of soft SUSY breaking terms 
using modular invariance of the type above has been ex- 
tensively discussed in the literature assuming moduli sta- 
bilization. In such analyses one generically finds that CP 
is indeed violated if one assumes that the moduli are in 
general complex. 

However, minimization of the potential and stabiliza- 
tion of the dilaton VEV is a generic problem in such mod- 
els and requires additional improvements. Often this is 
accomplished by non-pcrturbativc corrections to the po- 
tential. Thus one might consider non-perturbative con- 
tributions to the superpotential so that 

Wnp = f2((7)r?(T)-^ (62) 

Here t]{T) is the Dedekind function, and we have assumed 
a single overall modulus T, and a = S + 25'^^logriiT) 
and 5'~^^ — —{?>/ Ait)5'^^ . Additionally one can assume 
non-perturbative corrections to the Kahler potential 
and treat D{z) as a function to be determined by 
non-perturbative effects. The analysis shows that for a 
wide array of parameters minima typically occur at the 
self-dual points of the modular group, i.e., T = 1 and 
T = e^'^/^. However, for some choices of the parameters 
T can take complex values away from the fixed point. 
Nonetheless CP phases arising from such points arc very 
small since in the soft parameters they come multiplied 
by the function G{T,f) = (T + + 2dlog{r]{T)/dT 
the imaginary part of which varies very rapidly as 
the real part changes. Thus large CP phases do not 
appear to arise using the moduli stabilization of the type 
above(Bailin et ai, 1997). 

The situation changes significantly if Wnp contains an 
additional factor H{T) where 




where Gn{T) and G%{T) are Eisenstein functions of mod- 
ular weight 4 and 6, m, n are positive integers and P{j) 
is a polynomial of j{T) which is an absolute modular 
invariant. Alternately H can be expressed in the form 

H{T) = {j - 1728)"*/2jn/2p^^.^ ^g4^ 

The form on H{T) is dictated by the condition that no 
singularities appear in the fundamental domain. In this 
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case to achieve dilaton stabilization with T modulus not 
only on the boundary of the fundamental domain but 
also inside the fundamental domain and thus T has a 
substantial imaginary part. In this case it is possible 
to get CP phases for the soft parameters which can lie 
in the range 10"'' - 10-i(BaiIin et at, 1997). Thus 
with the absolute modular invariant in the superpoten- 
tial large CP phases can appear in the soft breaking 
in orbifold compactifications of the type discussed above. 

In the analysis of (Faraggi and Vives, 2002) the issue of 
CP violation and FCNC in string models with anomalous 
?7(l)yi-dilaton supersymmetry breaking mechanism was 
investigated. Here scalar masses arise dominantly from 
the U{1)a contribution while the dilaton generates the 
main contribution to the gaugino masses. Further, the 
dilaton contributions to the trilinear terms and to the 
gaugino masses have the same phase. In this class of 
models the nonuniversal components of the trilinear soft 
SUSY breaking parameter are typically small and one has 
suppression of FCNC and of CP in this class of models. 



C. CP violation on D brane models 

Considerable progress has occurred over the recent 
past in the development of Type I and Type II string 
theory. Specifically D branes have provided a new and 
better understanding of Type I string theory and con- 
nection with Type IIB orientifolds. Further, the advent 
of D branes open up the possibility of a new class of 
model building (for recent reviews on D branes see (Blu- 
menhagen et al, 2005, 2006; Polchinski, 1996)). Thus 
a stack of N D branes can produce generally an SU (N) 
gauge group or a subgroup of it, and open strings with 
both ends terminating on the same stack give rise to a 
vector multiplet corresponding to the gauge group of the 
stack. Further, open strings begiiming on one end and 
ending on another transform like the bifundamental rep- 
resentations and can be chiral. Thus these are possible 
candidates for massless quarks, leptons, and Higgs fields. 
A simple possibility for model building occurs with com- 
pactification on jZi x Zi. In addition to the axion- 
dilaton field s the moduli space consists in this case of 
the Kahler (t„i) and the complex strTicture (?i„i) mod- 
uli (m=l,2,3). For the moduli fields one has the Kahler 
potential of the form 



Kq = -ln{s + s) - ^ ln{tm + tm) - 

m—l 

3 

^ln{Um + Um)- (65) 



m=l 



Consider now complex scalars c]^^' along the direction i 
with ends of the open string ending in each case on a D9- 

brane. In this case one can obtain the Kahler potential 
including the complex scalar field by the translation tm + 



t„i t„i + tm — \cm^\'^- For the case of strings with 
both ending on the same Dbi brane one can show using 
either T-duality(Ibanez et al., 1999) or by use of Born- 
Infeld action(Kors, 2006: Kors and Nath, 2004) that the 
Kahler potential is modified by making the replacement 
s + s^s + s- ici^""^"*'!^. For the case when one has 



both D9- branes and D5r, 
potential reads 



-branes the modified Kahler 



^[99+55] = _ In + ^ _ ^ 12) - 

m=l 

^ In {tm + ir, - I' - ^ E In^uACt'^^ (66) 



m=l 



To construct the Kahler potential for the case when one 
has open strings with one end on D9-branes and the 
other end on D5m— branes, or for the case when open 
strings end on two different D5 branes, one can use the 
analogy to heterotic strings with Z2 -twisted matter fields 
(Ibanez et al., 1999; Kors and Nath, 2004). Alternately 
one can use string perturbation theory (Bertolini et al., 
2006; Lust et al., 2004, 2005). The result is 



5,„1|2 



m,n,p=l 
1 ^ 



{tn + in)^'^{tp + ipY/^ 
|C^[5™5„]|2 



m.n.p—1 



{tp + tp)y^{s + sY/^ 



(67) 



Explicit formulae for the soft parameters using these re- 
sults are given in the literature. However, one needs to 
keep in mind the configurations of the type discussed 
above are the so called \BPS states, and in this case 
the spectrum of open states falls into N = 2 multiplets, 
which implies that the spectrum is not chiral. Similar 
considerations apply to open strings which start and end 
on D3 and D'j branes, and results for these can be ob- 
tained by using T dualities. 

For realistic model building one needs to work with in- 
tersecting D branes. Thus in Calabi-Yau orientifolds of 
Type IIA one has D6-branes that intersect on the com- 
pactified 6 dimensional manifold. Sometimes it is conve- 
nient to work in the T-dual picture of Type IIB strings 
where the geometrical picture of branes intersecting is re- 
placed by internal world volume gauge field backgrounds, 
called fluxes on the D9 and D5 branes. The fluxes JF™ 
where a labels the set of branes, are rational numbers, 
i.e., T'^ = / n^, in order to satisfy charge quanti- 
zation constrains. The fluxes determine the number of 
chiral families. Further, the condition that = 1 super- 
symmetry be valid is a further constraint on the moduli 
and the fluxes and may be expressed in the form (Bachas, 
1995; Berkooz et al, 1996; Kors and Nath, 2004) 



E 



~ n 



(68) 



m=l 
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In the presence of fliixes the gauge kinetic energy function 
fa is given by 

/a=fln(™)fs-l ^ 7™.pF(")FWi J . (69) 

m=l \ m,n,p=l / 

The computation of the Kahlcr metric for the case of 
an open string with both ending on some given stack a, 
Cm'^\ can be computed by dimensional reduction(Kors 
and Nath, 2004) or string perturbation theory (Lust et al., 
2004) and is given by 



I a 



[aa]|2 



mfa) 



{S + S)(tm + trn)(Um + "rn) 1 + AI"'^ ' 
n,p—l ^ ' ^ ' 

Now the technique above using the hcterotic dual or 
Born-Infeld works for ^BPS brane configurations. How- 
ever, for the bifundamental fields C^"^! that connect the 
different stacks of branes with different world volume 
gauge flux one needs an actual string perturbation cal- 
culation and here the result for the Kahler potential is 
(Lust et al, 2004) 



|C'[a*']| 



Ul=ii^m+um/'^' r(l-^i:))i/2 

(m) 



\ '''■J 



Using the above one can obtain explicit expressions for 
the soft parameters. These have been worked out in de- 
tail in several papers. One can count the number of CP 
phases that enter in the analysis. They are the phases 
arising from s,tm,Um (ni=l,2,3). These can be reduced 
with extra restrictions such as. for example, dilation dom- 
inance which would imply only one CP phase 7s. 



D. SUSY CP phases and the CKM matrix 

A natural question is if there is a connection between 
the soft SUSY CP phases and the CKM phase Sckm- 
A priori it would appear that there is no connection be- 
tween these two since they arise from two very different 
sources. Thus the Sckm arises from the Yukawa inter- 
actions (assuming there is no CP violation in the Higgs 
sector) which from the string view point originates at the 
point when the string compactifies from 10 dimensions 
to four dimensions. This is the point where we begin to 
identify various species of quarks and leptons and their 
couplings to the Higgs bosons. On the other hand soft 
SUSY phases arise from the spontaneous breaking of su- 
persymmetry and enter only in the dimension < 3 oper- 
ators. Thus it would appear that they arc disconnected. 
While this conclusion is largely true it is not entirely so. 



The reason for this is that in SUGRA models the trilin- 
ear soft term Aa/3j contains a dependence on Yukawas so 
that(Kaplunovsky and Louis, 1993; Nath et al., 1983) 



(72) 



Thus the phase of the Yukawa couplings enters in the 
phase of the trilinear coupling. However, the phase re- 
lationship between A and Y is not rigid, since even for 
the case when there is no phase in the Yukawas one can 
generate a phase of A, and conversely even for the case 
when 5c KM is maximal one may constrain A to have zero 
phase. Further, it is entirely possible that the Yukawa 
couplings are all real and Sckm arises from CP violation 
in the Higgs sector as originally conjectured (Lee, 1973, 
1974; Weinberg, 1976). A more recent analysis of this 
possibility is given in (Chen et al., 2007). 

On a more theoretical level it was initially thought that 
CP violation could occur in string theory in either of the 
two ways: spontaneously or explicitly (Strominger and 
Wittcn, 1985). However, it was conjectured later that 
CP symmetry in string theory is a gauge theory and it 
is not violated explicitly (Choi et al., 1993; Dine et al., 
1992). We do not address this issue further here. 



VIII. THE EDM OF AN ELEMENTARY DIRAC FERMION 

If the spin-1/2 particle has electric dipolc moment 
EDM df, it would interact with the electromagnetic ten- 
sor Ffj^i, through 



(73) 



(74) 



which in the non-relativistic limit reads 
£ — dfip\(T.EijjA 

where tpA is the large component of the Dirac field. The 
above Lagrangian is not renormalizable, so it does not 
exist at the tree level of a renormalizable quantum field 
theory. However, it could be induced at the loop level if 
this theory contains sources of CP violation at the tree 
level. Thus suppose we wish to determine the EDM of 
a particle with the field ipf due to the exchange of two 
other heavy fields: a spinor ipi and a scalar (pk- The 
interaction that contains CP violation is given by 

C = LikipfPLipi(t)k + RikilJjPR.'4^i(t>k + H.c. (75) 

Here C violates CP invariance iff Im{LikR*i.) ^ 0. A 
direct analysis shows that the fermion ■^j acquires a one 
loop EDM df which is given by 



'^f = le;^^"'^^*^^^^^^^*^^^^ + QfcB(^)),(76) 
where 



A{r) = 
B{r) 



1 2/nr 
2(i-r)2^^~''+T37) 
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2(1 -r): 



2rlnr ^ 
1 — r 



(77) 
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For the case of the charged lepton we find 

2 



e — char gino / ^ 



CtEM 



47rsm Ow 



^=l ' ^^<^ 
(78) 

where U and V are as defined in Appendix XVI. A 
and where T^i = {i^eU*2V*{) = \Ke\U}^2i^ lu- A direct 
inspection of T^i shows that it depends on only one 
combination, i.e., ^2 + ^/^ + ^// where the phase 9h conies 
from the Higgs sector and as discussed later is generated 
at the loop level. 



\1 



■0i 



The neutralino exchange contribution to the EDM of 
the fermion is as follows: 



^ f — neutralino / ^ 



At: sin^ 9 



where 



w 



k=l i=l 



FIG. 3 Contributions to the electric dipole moment of a lep- 
ton or of a quark from the exchange of the charginos, the 
neutralinos and the gluino. The internal dashed line in the 
loop is the scalar field (f>k, the solid line is the fermion field ipi 
and the external wiggly line is the external photon line. 



We will utilize this result in EDM analyses in the follow- 
ing discussion. 



{coXiiDf2k - KfXbiDfik) (80) 



KfXhiD*f2k) 



where ao = -\/2tan6'vK((3/ - Ts/), 5o = -V^Tsf, 
Co = V2ta.iiewQf, and in Xm b=3(4) for Tsq = -5(5). 
The following three combinations of phases appear in 
Vf^k■■ ii+0^ + Oh, ^2 + 0^ + Oh and af + e^ + en- We 
note in passing that the contribution from the neutrino 
Yukawa couplings to the lepton electric dipole moment is 
computed in (Farzan and Peskin, 2004), and the charged 
Higgs contributions to the lepton EDM in a two-Higgs 
doublet model is discussed in (Kao and Xu, 1992). 



IX. EDM OF A CHARGED LEPTON IN SUSY 

We discuss now the EDM of a charged lepton in MSSM 
using the results of the previous section. As mentioned in 
Sec. IV, in softly broken supersymmetric models as many 
as 40 additional phases can appear. However, only cer- 
tain combinations of phases appear in a given process and 
the number of such combinations depends on the process. 
We discuss now the details. 

In the computations here we use the Lagrangian of ap- 
plied N = 1 supergravity for the case of MSSM fields 
with inclusion of soft breaking(Haber and Kane, 1985; 
Nath et at; Nilles, 1984). The EDM of a charged lep- 
ton receives contributions from chargino, neutralino, and 
slepton exchanges. A discussion of the chargino and neu- 
tralino masses is given in Sec. (XVI. A) while a discussion 
of the slepton and squark masses is given in Sec. (XVI. B). 



X. EDM OF QUARKS IN SUSY 

The quarks receive contribution from the electric 
dipole operator (d^), from the chromoelectric dipole op- 
erator (c?^), and from the purely gluonic dimension six 
operator of Weinberg (d^). Thus 

^9 = < + < + (81) 
We discuss these in further detail below. 



A. The electric dipole moment operator contribution to 
EDM of quarks 

The electric dipole moment operator receives contri- 
butions from the gluino, chargino and neutralino ex- 
changes. The gluino exchange contributes to the EDM 
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of the quarks as follows 




(82) 



where qi and q2 are the mass eigenstates, and T^'^ 



g2fe-D*iA;, as 



m-g is the gluino mass, and -B(r) 



is as defined by Eq.(77). An explicit analysis gives Fj 



-FJ^ where 



Mr") = .^2 _%^2 ("^ol^gl sin(aq - 6) + 

|M|sin(0p + ^ff (83) 

which holds for both signs of — M'i^. It is easy to see 
that the combinations of phases that enter are (ag-^3) 
and ^3 + 6*^ + or alternately one can choose them to 
be aq + 6*^ + and ^3 + 6*^ + Qh- 

The chargino contribution to the EDM for the up 
quark is as follows 



B. The chromoelectric dipole moment contribution to tlie 
EDM of quarks 

For the case of the quarks one has two more opera- 
tors that contribute. These are the quark chromoelectric 
dipole moment (d*-') and the purely gluonic dimension 
six operator. For the operator we have the effective 
dimension five operator 



(88) 



where T" are the SU (3) generators. Contributions to d'^ 
of the quarks from the gluino, the chargino and from the 
neutralino exchange are given by 
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^q— neutralino 
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[<9jB(-^) + (Q„ - Oj)A(-^)](84) ^here B(r) is defined by Eq.(77) and C(r) is given by 



(91) 



Here A{r) is as defined by Eq.(77) and 



and 



-ien 



V2mw sin /? ' 



\plmw cos /3 



(85) 
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C(r) 



6(r- 1): 



:(10r-26- 



2rZnr 18/nr 



1 — r 1 — r 



)• (92) 



We note that all of the CP violating phases are contained 
in the factors Im(rg'°), Im(Fgjfe), and \vi\{r]qi].). But these 
are precisely the same factors that appear in the gluino, 
the chargino and the neutralino contributions to the elec- 
tric dipole operator. 



and explicitly 



rml(2) = \K,u\{cOS^Qdl'^WL2tUnv^ 



(87) 



The EDM here depends only on two combinations of 
phases: -h 6*^ + da and ^2+6'^ + with ^2 - o-i 
being just a linear combination of the first two. A sim- 
ilar analyses hold for the chargino contributions to the 
down quark and one gets only two phase combinations 
which are identical to the case above with aa replaced by 
(Xu- The neutralino exchange contribution to the EDM 
of quarks is given by Eq.(79). The sum of the gluino, the 
chargino and the neutralino exchanges discussed above 
gives the total contribution from the electric dipole op- 
erator to the quark EDM. 



C. The contribution of the purely gluonic operator to the 
EDM of quarks 

The purely gluonic dimension six operator which con- 
tributes to the dipole moment is (Weinberg, 1989) 



(93) 



where Ga^iv is the gluon field strength tensor. fafj~f are 
the Gell-Mann coefficients, and c'^"^'^ is the totally an- 
tisymmetric tensor with 6*^^^^ = -|-1. An analysis of dP 
including the quark-squark-gluino exchange (sec Fig. (4) 
where one of the loops contributing to this operator is 
shown) with gluino phase ^3 but with squark mass^ ma- 
trix treated real is given in (Dai et ai, 1990). Including 
the phases from At and fi in the squark mass^ matrix the 
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the following expression 




FIG. 4 The quark-squark-gluino exchange contribution to the 
purely gluonic dimension six operator. The dashed line in 
upper semi circle in the loop is the squark g, the internal 
horizontal solid line is the gluino g, the solid line on the lower 
semicircle in the loop is the quark q, while the external wiggly 
lines are the gluons. 



analysis of dP gives (Dai et ai, 1990; Ibrahim and Nath, 
1998a) 



47rm 



+m,iz1 - z!^)Im.{Tl^)H{zl z^ z,)) 



Here 



^ rrig nig 



(95) 



and H{zi, Z2, z^) is defined by 
H{zi, Z2, Z3) = — j dx 1 du j dyx{l 



where 



^1 



N1N2 



u{l — x) + Z3X{1 — x){l ~ u) — 
-2ux[ziy + 22(1 - y)], 
1 



iVa = (1 - xYil - uY +u^- -x\l - u) 

9 



|/(z|), 



(98) 



where I{z) is defined by 
1 



+ (1 - 16z - 'dz'^)lnz]. 



(99) 



The contribution of the last two operators to the EDM 
of the quarks can be computed using the naive dimen- 
sional analysis(Manohar and Georgi, 1984). This tech- 
nique can be expressed in terms of a rule using the 're- 
duced' coupling constants. Thus for example, for a cou- 
pling constant g appearing in an interaction of dimen- 
sionality (mass)'^ and containing N field operators the re- 
duced coupling is {^ivY'^ M^^^g where M is the chiral- 
symmetry breaking scale and has the value M = 1.19 
GeV. Thus the rule means that the reduced coupling of 
any term in the effective hadronic theory at energies be- 
low M is given by a product of the reduced couplings 
of the operators appearing in the effective Lagrangian at 
energies below M, that produces this term. Using this 
rule for the chromoelectric and purely gluonic dimension 
six operators one finds there contribution to the EDM of 
the quarks are given as follows 



47r ' 



(100) 



where 77^, rf^ and rf^ are renormalization group evolu- 
tion of dq, d^ and d^ from the electroweak scale to the 
hadronic scale. A discussion of how these renormalization 
group factors are computed is discussed in Sec. (XVI. C). 
Their numerical value is estimated to be rj^ « 0.61 (De- 

3.4. The alternate tech- 



(94) grassi et ai, 2005), r]^ « 77*^ 

nique to estimate the contributions of the chromoelec- 
tric operator is to use the QCD sum rules (Khriplovich 
and Zyablyuk, 1996). To obtain the neutron EDM, we 
use the non-relativistic SU (6) quark model which gives 

J — 4 J _ 1j 



D. The cancelation mechanism and other remedies for the 
CP problem in SUSY, in strings and in branes 

Thus MSSM contains new sources of CP violation 
and these phases would induce EDMs of the fermions 
in the theory. Taking the values of the parameters of 
the model at their phenomenologically favorable range 
(toi/2 ~ Too ~ lOOGeV^, tan/3 ~ 10,0^ ~ ao 1) one 
finds that the EDMs of the electron and neutron exceed 
the experimental bounds by several orders of magnitude. 
This problem is certainly a weakness of the low energy 



D = u{l - x) + Z3x{l -x){l-u)+ ux[ziy + 22(1 - ?/)]. (97)SUSY and needs to be corrected to make the theory vi- 
able. Various remedies have been suggested in the liter- 
For the case when m^, mg >> one obtains for H ature to overcome this problem. The first of these is the 
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suggestion that the first generation of sleptons and the 
first two generations of squarks are very heavy (Nath, 
1991) (see also (Kizukuri and Oshimo, 1992)). This 
means the production and study of these particles at 
LHC will be difficult if not impossible. There is another 
reason that this possibility is not attractive is that, the 
annihilation rate of the lightest supersymmetric particle 
LSP may be too low in this range of masses and as a 
result the relic density of the LSP may be larger than 
the observed dark matter density. Another suggestion is 
that the phases are small O(10~^) (Arnowitt et ai, 1991, 
1990; Braaten et ai, 1990a,b; Dai et ai, 1990; Dugan 
et ai, 1985; Ellis et ai, 1982; Franco and Mangano, 
1984; Garisto and Wells, 1997; Gunion and Wyler, 1990; 
Polchinski and Wise, 1983; Weinberg, 1989). However, 
a small phase constitute a fine tuning and there will not 
be any interesting display of CP violation in colliders. 
Moreover, electroweak baryogenesis cannot take place in 
this case(Kuzmin et ai, 1985). A third possibility first 
proposed in (Ibrahim and Nath, 1998a, b,c) is that there 
are internal cancelations among the various contributions 
to the neutron EDM and to the electron EDM, leading 
to compatibility with experiment with large phases and 
a SUSY spectrum that is still within the reach of the 
accelerators. 

This is the most interesting solution because it leaves 
room for a host of non trivial CP violating as well as 
CP conserving phenomena to be discovered at colliders 
and elsewhere. By CP violating properties we mean those 
properties that vanish in the limit of CP conservation like 
the EDMs and the neutral Higgs bosons mixing. By CP 
conserving phenomena, we mean those properties that 
exist in the absence of CP violation but they differ if 
CP violation is included like — 2. Following the work 
of (Ibrahim and Nath, 1998a, b,c) there is much further 
work on the cancelation mechanism in the literature(Abel 
et ai, 2002; Accomando et ai, 2000a, b; Barger et ai, 
2001; Bartl et ai, 1999, 2001; Brhlik et ai, 2000a, 1999a, 
2000b, 1999b, 2001; Chattopadhyay et ai, 2001; Falk and 
Olive, 1998; Falk et ai, 1999b; Ibrahim and Nath, 1998c, 
2000d; Pokorski et ai, 2000). 

As was shown above, the quark and the lepton EDMs 
in general depend on ten independent phases providing 
one with considerable freedom for the satisfaction of the 
EDM constraints. Numerical analyses show the exis- 
tence of significant regions of the parameter space where 
the cancelation mechanism holds. We describe here a 
straightforward technique for accomplishing the satisfac- 
tion of the electron EDM and the neutron EDM con- 
straints. For the case of the electron one finds that the 
chargino component of the electron is independent of 
and the electron EDM as a whole is independent of ^3. 
Thus the algorithm to discover a point of simultaneous 
cancelation for the electron EDM and for the neutron 
EDM is a straightforward one. For a given set of param- 
eters we vary ^1 until we reach the cancelation for the 
electron EDM since only one of its components (the neu- 
tralino) is affected by that parameter. Once the electric 
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FIG. 5 Two loop Barr-Zee type diagrams that contribute to 
the EDMs in supersymmetry (Chang et al., 1999). 



dipole moment constraint on the electron is satisfied we 
vary ^3 which affects only the neutron EDM keeping all 
other parameters fixed. By using this simple algorithm 
one can generate any number of simultaneous cancela- 
tions. The EDM of the atoms also provide a sensitive 
test of CP violation. An example is the EDM of Hg- 
199 for which the current limits are given by Eq.(22). 
Among the phases that enter the EDM of IIg-199 is the 
phase as- We note that enters only in to one 
loop order, and thus it can be varied to achieve a simul- 
taneous cancelation in dji and a consistency with the 
experimental limits. Illustrative examples of points in 
the parameter space where cancelations occur and all the 
EDM constraints are satisfied are given in Tables 1 and 
2 in Sec. (XVI. D). It needs to be emphasized that while 
the cancellations among the various contributions to the 
EDMs are pretty generic the suppression of the EDMs 
for the electron and for the neutron do require fine tun- 
ing. On the positive side the above, of course, leads to a 
narrowing of the parameter space of the theory. 

In theories where the Higgs mixing parameter ^ obeys 
the simple scaling behavior as the rest of the SUSY 
masses the EDMs exhibit a simple scaling behavior under 
the simultaneous scaling on mg and mi/2- In the scaling 
region the knowledge of a single point in the MSSM pa- 
rameter space where the cancelation in the EDMs occurs 
allows one to generate a trajectory in the toq — plane 
where the cancelation mechanism holds and the EDMs 
are small. Thus under the transformation toq — > Amo, 
™i/2 Atoi/2 M itself obeys the same scaling, i.e., 
/i — > A/i in the large fi region. In this case d^ exhibits the 
scaling behavior 

4 -> A"^4 (101) 

The same scaling relation holds for the electric and for 
the chromoelectric operators of the quarks 

df-^X-^df, d^->A~2^^ (102) 
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For the gluonic dimension 6 operator we find the follow- 
ing scaling 



4jG 



(103) 



Thus the scaling property of dq will be more complicated. 
However, as A gets large the contribution of will fall 
off faster than and dg' and in this case one will have 

A~^d„. Thus scaling 



A '^dq and so dr, 



the scaling dq 

property of EDMs allows one to promote a single point 
in the SUSY parameter space where cancelation occurs 
to a trajectory in the parameter space. With the scal- 
ing property one can arrange the cancellation mechanism 
to work for the EDMs over a much larger region of the 
parameter space (Ibrahim and Nath, 2000d) than would 
otherwise be possible (Pospelov and Ritz, 2005). The 
scaling phenomenon also has implications for the satis- 
faction of the EDM constraints in string and D-brane 
models (Ibrahim and Nath, 2000d). As stated already in 
general only certain phase combinations appear in the 
analysis of a given physical quantity. Some examples of 
such combinations are given in Table 3 in Sec. (XVI. E). 
For other solutions to the SUSY CP problem see (Abel 
et al., 2001; Babu et ai, 2000b; Dimopoulos and Thomas, 
1996; Nir and Rattazzi, 1996). 



E. Two loop contribution to EDMs 

Two loop contributions to the EDMs can be quite sig- 
nificant. Thus the analysis of (Barr and Zee, 1990; Gu- 
nion and Wyler, 1990) showed that significant contribu- 
tions to the EDM of the electron and of the neutron can 
result if the Higgs-boson exchange mediates CP viola- 
tion. A more recent analysis in the same spirit is given 
by (Chang et al, 1999) for the MSSM case. Here the CP 
phases arising from the Higgs boson couplings to the stop 
and the sbottom enter and these are not stringently con- 
strained by data. Thus CP phases in the third generation 
could be quite substantial consistent with the EDM con- 
straints. We discuss now the two loop analysis in further 
detail. We assume that the large CP phases arise only in 
the third generation trilinear soft parameters A-^^u, and 
the relevant two loop interactions arise via the CP-odd 
Higgs a{x) (see Fig. (5)) whose interactions with fermions 
and sfermions are given by 
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RfiaHf + v£,fa{-flh + /2V2), (104) 



where g is related to the W boson mass by = gvjl, 
Rf = cot/3(tan/3) for tI = |(-^). The diagrams of 
Fig. (5) give the following contribution to the EDM of a 
fermion at the electroweak scale 
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FIG. 6 Estimate of the size of two loop contribution to the 
EDMs in supersymmetry with phases only in the third gen- 
eration (Chang et al., 1999). 



where Xia = {mqjma)'^ (i=l,2), (q=t,b) are defined 



by 



sin 2/3 ' 
mt sin 2^t-fm(/ie'*') 



6 = 



sin^ /? 



(106) 



and where 5q = arg{Aq + Rq^i*). The function F{x) is 
given by the loop integral 

Similarly the contribution to CEDM at the electroweak 
scale is given by 

" q=t,b 

A numerical analysis of the EDM is given in Fig. (6) and 
indicates that one can satisfy the EDM constraints in 
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certain ranges of the parameter space. However, it re- 
mains to be seen how one can naturaUy suppress phases 
in the first two generations while aUowing them only in 
the third generation. The reader is also directed to sev- 
eral other works on two loop analyses of EDMs: (Chang 
et al, 1990, 1991; Degrassi et al, 2005; Feng et ai, 2006, 
2005; Pilaftsis, 2002). Specifically, a complete account of 
all dominant 2-loop Barr-Zee type graphs in the CP vi- 
olating MSSM is given in in (Pilaftsis, 2002). The anal- 
yses of EDMs given in this section were based on the 
assumption of R parity conservation. For analyses of 
EDMs without R parity see (Faessler et al, 2006; Hall 
and Suzuki, 1984; Keum and Kong, 2001a, b). 
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XI. CP EFFECTS AND SUSY PHENOMENA 



As noted earlier with the cancelation mechanism the 
phases can be large, and thus their effects could be visible 
in many supersymmctric phcnomena(Aoki et oi., 1999; 
Asatrian and Asatrian, 1999; Baek and Ko, 1999; Barr 
and Khalil, 2000; Choi and Drees, 1998; Choi et al., 
2000b; Choi and Lee, 2000; Choi et a/., 2000c; Dedes 
and Moretti, 2000a,b; Goto et al., 1999; Huang and Liao, 
2000a,b, 2002; Kneur and Moultaka, 2000; Kribs, 2000; 
Ma et ai, 1999; Mrcnna et ai, 2000; Okada et al, 2000). 
Below we discuss several of these phenomena and refer 
to the literature above for others. 
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A. SUSY phases and — 2 

The effects of CP violating phases on the supersym- 
mctric electroweak contributions to g'^ — 2 have been 
investigated(Ibrahim et al., 2001; Ibrahim and Nath, 
2000a, c). The parameter = ^"^"^ is induced by loop 
corrections to the muon vertex with the photon field. In 
MSSM the muon interacts with other fermions ■i/'j and 
scalars <f)k through 

C = LikJj.PLi'i(j>k + RikJiPRi'i(l)k + H.c. (109) 

where ipi stands for the neutralino (chargino) and 0k 
stands for the smuon (scalar neutrino). The one loop 
contribution to is given by 



(110) 



Here comes from the neutralino exchange contribution 
and comes from the chargino exchange contribution 
so that 
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FIG. 7 Chargino and neutralino exchanges contributing to 
the muon g-2 which generate dependence oi — 2 on phases. 



and 
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3:-,{\Lik? + \Rik?)h{'^,'^)- (112) 
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In the supersymmctric limit the soft breaking terms van- 
ish and should vanish as well(Barbieri and Giudice, 
1993; Ferrara and Remiddi, 1974). A careful limit of 
Eqs.(lll) and (112) shows that in the supersymmctric 
limit the sum of the W exchange contribution, in the 
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tz (rad.) 

FIG. 8 An exhibition of the dependence of on a SUSY 

CP phase. The curves correspond to the four cases below 
(Ibrahim and Nath, 2000c) : (1) mo=70, mi/2=99, tan/3=3 
, |^o|=5.6, 6=-l, 6=0.62; 9^^ 2.35, aAo = -4; (2) mo=80, 
mi/2=99, ton/3=5 , |^o|=5.5, ,^i=-0.8, 6=0.95; 6^= 1.98, 
aAo= 0.4; (3) mo=75, mi/2=132, ton/3=4 , |Ao|=6.6, 6=-l, 
6=2.74; 1.2, = -1.5; (4) m(,=70, mi/2=99, tan/3=6 
, |Ao|=3.2, $1=0.63, 6=0-47, 6l^= 2.7, uaq = --4 where all 
masses are in GeV units and all phases are in rad. 



standard model part, and of the chargino exchange con- 
tributions, in the supersymmetric counterpart, cancel. 
Thus 







(114) 



Similarly one can show that the Z boson exchange and 
the contribution of the massive modes of the neutralino 
sector in the supersymmetric limit cancel. 



when loop corrections to the effective potential arc in- 
cluded. To calculate such mixings we use the one loop 
effective potential as given by Eq.(28). We assume that 
the SU{2) Higgs doublets Hi^2 have non- vanishing vac- 
uum expectation v\ and so that we can write 



fu \ ^ ( vi + (f)i + ixpi 



-v/2 \V2 + 4>2+ ill>2 



(116) 



For the present case with the inclusion of CP violat- 
ing effects, the variations with respect to the fields 
<l>i,4>2, i'1,'4'2 give the following 

^ ,dAV 2.9l+gl,2 2x, 2. a Q 

(-53— )o = + — -{vi — V2)+m^ tanpcos^ij, 



vi d(j)\ 



1 ,aAF, 
— )o = m 



(117) 



V2 d(j)2 



2 92 9l / 2 2\ , 2 ,r, n 

— ^ — K - ^^2) + m^cotp COS Oh, 

(118) 



1 ,dAV. 

Vi Otp2 



TOg sm Oh 



1 ,dAV. 

V2 Olpi 



(119) 



where the subscript means that the quantities are 
evaluated at the point (pi = (t>2 = 'fpi = 4'2 = 0- As 
noted in (Demir, 1999) only one of the two equations in 
Eq.(119) is independent. 



(massive) = 



One can show that the massless part of the neutralino 
spectrum in the supersymmetric limit gives the value of 
— aeTO/27r. Thus it gives the same magnitude but is op- 
posite in sign to the famous photon exchange result. 

The CP dependence of arises from the effect of the 
phases on the sparticlc masses, and on their effects on La- 
and Rik and significant variations can arise in as the 
phases are varied. An illustration of this phenomenon is 
given in Fig. (8). Because of the significant dependence 
of on the phases it is possible to constrain the CP where 
phases using the current data on a^. This is done in 
(Ibrahim et ai, 2001). Further details on the analysis of 
this section are given in Sec. (XVI. F). 



(115) One can have sizable contributions to the potential 

corrections from top-stop, bottom-sbottom (Choi et ai, 
2000a; Demir et ai, 2000b; Ibrahim and Nath, 2001a; Pi- 
laftsis and Wagner, 1999), W — — sector (Ibrahim 
and Nath, 2001a) and from the x° - Z - h° - H° sector 
(Ham et ai, 2003; Ibrahim and Nath, 2002). The mass- 
squared matrix of the neutral Higgs bosons is defined by 



B. SUSY CP phases and CP even -CP odd mixing in the 
neutral Higgs boson sector 

Another important effect of CP violating phases is 
their role in determining the spectrum and CP properties 
of the neutral Higgs fields arising due to mixings of the 
CP even-CP odd Higgs(Pilaftsis, 1998a,b; Pilaftsis and 
Wagner, 1999). 

Such mixings between CP even and CP odd Higgs 
bosons cannot occur at the tree level, but are possible 



(a = 1 — 4) are defined by 

{*a} = {01,02,V'1,V'2}, 



(120) 



(121) 



and the subscript means that we set = 4)2 = ipi — 
il>2 = 0. The dominant contributions come from the stop, 
sbottom and chargino. With the inclusion of the stop, 
the sbottom, and of the chargino contributions one finds 
that 6h is determined by 



2/-'>tir-M--H i-Ji\---ti' ---ts/ ' 2' 

,2 2 \ 92 I. .11- I „; j: /_2 2 



|Ab|sin7b/i(mj^,m^J - Yg;^|/^||TO2| sin 72/1(771^^, m^J, 

(122) 
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where 



7, Phi, 



3hl 



167r2' 167r2' 
7t = o:At + 0^1, lb = oiA^ + 9^, 72 = 6 + 9^i, (123) 



and /i (a;, t/) is defined by 



h{x,y) = -2 + ;og^ + ^— ^Zoff 



xy v + x y 
y - x a; 



(124) 



The inclusion of the stau and the neutralino sectors in 
the analysis would contribute extra terms to Eq.(122) 
that are dependent on the phase 7r = + ^/x and 
7i = Ci + ^/n- The tree and loop contributions to M^^ are 
given by 



(125) 



where M^^^^ are the contributions at the tree level and 
AM^^ are the loop contributions where 



327r2 'd^a d'^b 



and where e=2.718. Computation of the 4x4 Higgs 
mass^ matrix in the basis of Eq.(121) gives 

f Afii+Aii -M12 + A12 A13S/3 AisCp \ 
-M12 + A12 M22 + A22 A23S/3 A23C/3 

A13S/3 A23S/3 M33SI M33S/3C;3 

y Ai3C^ A23C^ M33S0C/3 M33CI y 

(127) 

where Mn = Mfc^ + M^sj, M12 = (M| + M^sisc^, 
M22 = M|s| + M^f:|. Cff.sp = cos /3, sin /3, and M33 = 
M\ + A33. and where (c^jS^) = (cos /3, sin /3). Here the 
explicit Q dependence has been absorbed in m\ which is 
given by 



m\ = (sin (3 cos /?) ^ (— TO3 cos ^ + 
J^t I cos jtfi {m\ , m| ) 

+ 7;Ph,\Ai,\\fi\ cos76/i(m? m? ) 



+ 2^'»x l-4r I ImI cos7t-/i (m?^ , m?J 
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^^|to2||/x| cos72/i(TO^+,m^+)) + A^. 



(128) 
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FIG. 9 An exhibition of the phenomenon of CP even -CP 
odd Higgs mixing via the SUSY CP phases. The CP even 
component 0i of Hi (upper curves) and the CP odd compo- 
nent ViD of Hi (lower curves) including the stop, sbottom, 
stau, chargino and neutraUno sector contributions as a func- 
tion of 9^. The common parameters are ruA = 300, Q = 320, 
mo = 100, mi — 500, ^2 = .5, ao = .3, |^o| ~ 1. For curves 
with diamonds tan/3 = 15, ^1 = 1.5, for circles tan/3 = 8, 
^1 = 1.5, and for triangles tan/3 = 8, = 0.5 where all 
masses are in GeV and all angles are in radians. 



2)0, (126) exchange and 



167r2 



i=i 



Ml. Ml 

^(M^)-i) 



+Ml{gl{\m,\ 
+<??(|m2p 

-gl\mi\^\n\^\rh2\ cos 72 - gl\m2\^\fj,\^\rhi\ cos72)(129) 



'- |/xp)|/i||m2|cos72 
|^P)|^||mi|cos7i) 



where A^ = ^1 — ^2- The first term in the second brace 
on the right hand side of Eqs.(128) is the tree term, while 
the second, the third and the fourth terms come from the 
stop, sbottom, stau and chargino exchange contributions. 
The remaining contributions in Eq.(128) arise from the 



neutralino sector. For Ay one has 



A,:. = A, 



ijt 



(130) 



Here A^ is the contribution arising from the neutralino {Mho , M^o ) — {Mz , 0) 



where A^jf is the contribution from the stop exchange 
in the loops, A^j^ is the contribution from the sbottom 
exchange in the loops, Aijf is the contribution from the 
stau loop, Aij^+ is the contribution from the chargino 
sector and A^j^o is the contribution from the neutralino 
sector. For illustration A^J^ are listed in Sec. (XVI. G). 

We note that the phases come to play a role here 
through the squark, slepton, chargino and neutralino 
eigen values of their mass matrices. We note that 
in the supersymmetric limit M^a = (0,0, Mz, Mz), 
M + =Mh+ = M, 



w and Mq- = 
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iTiq. With this in mind one can sec that aU the radiative 
corrections to the potential vanish in the supersymmetric 
limit. By introducing a new basis <pi,<p2,'4>iD,'4>2D where 



ipiD = sin/3V'i + cos/3^/'2 
V'2£) = - cos/3V'i + sin/3V'2, 



(131) 



one finds that the field tp2D decouples from the other 
three fields and is a massless state (a Goldstone field). 
The Higgs mass^ matrix -M^.^^^ of the remaining three 
fields is given by 



Mil + All 

-Mi2 + Ai2 
Ai3 



-Mi2 + Ai2 

M22 + A22 
A23 



Ai3 

A23 

Ml + A33 



(132) 



We note that the basis fields {(j)i,(j)2,ipiD} of the above 
matrix are the real parts of the neutral Higgs fields and 
a linear combination of their imaginary parts ipi . Thus 
these states are pure CP states where ^1,2 are CP even 
(scalars) and Vid is CP odd (a pseudoscalar). What 
we are interested here is the mixing between the CP 
even and the CP odd Higgs states in the cigcn vectors 
of the above matrix and this mixing is governed by the 
off-diagonal elements A12 and A23. These are found to 
be linear combination of sin7(, sin7f,, sin 71, sin 72, and 
sin7T- where these phases are defined as in Eq.(123). In 
the limit of vanishing CP phases the matrix elements A12 
and A23 vanish and thus the Higgs mass^ matrix factors 
into a 2 X 2 CP even Higgs matrix times a CP odd el- 
ement. The effect of phases on CP even-CP odd Higgs 
boson mixings have been studied by (Choi et al., 2000a; 
Demir et al, 2000b; Ibrahim and Nath, 2001a; Pilaft- 
sis and Wagner, 1999) and found to be significant. It is 
shown that if a mixing effect among the CP even and the 
CP odd Higgs bosons is observed experimentally, then 
it is only the cancelation mechanism of EDMs that can 
survive(Ibrahim, 2001a). A more accurate determination 
of the VEV of the Higgs fields would require use of two 
loop effective potential. An improved accuracy and scale 
dependence should be obtained with the full two-loop ef- 
fective potential (Martin, 2003). 



C. Effect of SUSY CP phases on the b quark mass 

The running b quark mass is another object in MSSM 
where CP phases could have an impact, nib can be writ- 
ten in the form 

rribiMz) = hbiMz)^ cos P{1 + Ab) (133) 
v2 

where hb{Mz) is the Yukawa coupling for the b quark at 
the scale Mz and Ab is the loop correction to rub- The 
SUSY QCD and electroweak corrections are large in the 
large tan/3 region(Carena et al., 1994; Hall et al., 1994; 
Pierce et al, 1997). At the tree level the b quark couples 
to the neutral component of Hi Higgs boson while the 



coupling to the H2 Higgs boson is absent. Loop correc- 
tions produce a shift in the couplings and generate a 
non- vanishing effective coupling with i/j ■ Thus the effec- 
tive Lagrangian would be written as (Babu et al., 1999; 
Carena and Haber, 2003) 

-jCeff = {hb + 5hb)bRbLH^ + AhbbRbLH^* + H.c{134) 

where the star on H2 is necessary in order to have a 
gauge invariant Lagrangian. The quantities 6hb and Ahb 
receive SUSY QCD and SUSY electroweak contributions. 
The QCD contribution arises from the corrections where 
gluinos and sbottoms are running in the loops. In the 
elcctrowak contributions, the sbottoms (stops) and the 
neutralinos (charginos) are running in the loops. The 
basic integral that appears in the expressions of 5hb and 
Ahb involving heavy scalars Si, S2 and a heavy fermion 
/is 



rrif- 



■7m 



(27r)4 (A;2 - mj){k^ - m|J(A;2 _ to|J 



(135) 



In the approximation of the zero external momentum this 
integral could be written in the closed form 



(47r)- 



(136) 



where the function f{rn?, mf,rrij) is given by 

/ TM^ Tfl^ 771^ \ 

X [ m^m^ ^^^'^ m^m- In + m- m- In ) (137) 



for the case i ^ j and 



2 2 2\ 
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(to^ In ^ -I- (to^ - uii)) 



(to? — TO^)2 



TO 



for the case i = j. In the SUSY QCD the heavy fermion is 
the gluino and the heavy scalars are the sbottoms. In the 
chargino contribution, the chargino is the heavy fermion 
and the heavy scalars are the stops. In the ncutralino 
part, the neutralino is the heavy fermion and the heavy 
scalars are the sbottoms. 

The couplings Shb and Ahb are generally complex due 
to CP phases in the soft SUSY breaking terms. Elec- 
troweak symmetry is broken spontaneously by giving ex- 
pectation values to Hi and H2. Thus one finds for the 
mass term 



-Crn = MbbRbL+H.C, 



(138) 



where 



hbV cos B,^ Shb Ahb /,„„n 
Mb= ^ ' {1 + ^ + --^ tan 13) (139) 



V2 



hb hb 
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Here M^, is complex. By rotating the b quark field 



b = e'/^'^'-^'b', tan 



Xb 



ImMb 
ReMh 



one gets 



= mbb'nb'^ + H.c, 



(140) 



(141) 



where 



1 



Cl = -^{Re{hb + 5hb)Rii + {-Im{hb + 6hb) sin/3 
v2 

+Im{Ahb) cos/3}i?23 + Re{Ahb)R22) 

cl = — \={-Irn{hb + 6hb)R2i + {-ReQib + 5hb) sin/3 
V2 

+Re{Ahb) cos/3}i?23 - Im.{Ahb)R22) 



where rub is real and positive and b' is the physical field. 

hbV cos P - , 2 -ON i 

S^ = Re{^-^)+Re{^)t&nf3 
hb hb 

5j = Im{^) + Im{^) tan p (142) 
hb hb 



Thus one finds for the mass correction 



Re— — tan p + Ke- — 
hb hb 



(143) 



The SUSY CP violating phases in the SUSY QCD correc- 
tions are ^3, a a,, and 0^. These come from the vertices of 
bbg and bbH. In the chargino part one finds the phases ^2, 
a At and 0^. In the case of ncutralino wc have ^2, Cij 01 
and 9^. The corrections of the b quark mass are found 
to be very sensitively dependent on 0^, ^3 and a^o as the 
values of these phases affect both the sign and the magni- 
tude of the correction. Thus the correction can vary from 
zero to as much as 30% in some regions of the parame- 
ter space and can also change its sign depending on the 
value of these phases. The effect of ^2 is less important 
and ^1 is found to be the least important phase(Ibrahim 
and Nath, 2003c). Similar results hold for the r lepton 
mass and for the top quark mass. For the r lepton the 
numerical size of the correction is as much as 5% and for 
the top quark is typically less than a percent. 



D. SUSY CP phases and the decays ft — » 66, ft. — > rr 

As was mentioned above, the spectrum of the neutral 
Higgs sector and its CP properties are sensitive to the 
CP violating phases through radiative corrections. The 
couplings of the quarks with the Higgs are also found to 
be dependent of these phases. Thus one can deduce the 
corrected effective interaction of the b quark with the 
lightest Higgs boson H2 as 



-Cint = b{Cl +i^zCP)bH2 



where 



cosxf, -smxh 
sinxb cosxb 



Cl 



(144) 



(145) 



The matrix R is the diagonalizing matrix of the Higgs 
mass^ matrix 

RMl^gg, R^ = diag{m%^ , , ml^ ) (146) 

where we use the convention that in the limit of vanish- 
ing CP phases, one has Hi ^ H, H2 ^ h and A. 
These elements Rij and the corrections Shb and Ahb are 
found to be sensitive functions of the CP violating phases 
and their values are all determined by SUSY radiative 
corrections of MSSM potential. The quantity i?b/r de- 
fined as 



R 



b/7 



BRQi bb) 
BR{h f r) 



(147) 



is found to be an important tool to discover supersym- 
metry. In Standard Model, it is given by 



m. 



(148) 



where (1 + w) is the QCD enhancement factor (Gorishnii 
et al, 1990). 



1 -I- w = 1 -I- 5.67— -h 29.14 



(149) 



By identifying nih with , the lightest Higgs boson in 
MSSM, we find a shift in iJ^/r value due to supersym- 
metric effect including the eflFects due to CP phases as 
follows 



ARb/r = 



R> 



■b/T 



r/SM 



R> 



■b/T 



(150) 



The quantity Rb/r in MSSM depends on the CP phase 
via and . Thus if a neutral Higgs is discovered and 
Rb/r measured and found to be different from what one 
expects in the Standard Model, then it would point to 
a non-standard Higgs boson such as from MSSM (Babu 
and Kolda, 1999). The analysis of (Ibrahim and Nath, 
2003b) that the supersymmetric effects with CP phases 
can change the branching ratios by as much as 100 % for 
the lightest Higgs boson decay into bb and fr. Similar 
results are reported for the other heavier Higgs bosons. 
Thus the deviation from the Standard Model result for 
Rb/r depends on the CP phases and it can be used as 
a possible signature for supersymmetry and CP effects. 
Similar analyses can also be given for the decay of the 
heavy Higgs, e.g., for tt,bb and to x~^X~ (Eberl 

et al., 2004; Ibrahim, 2007) if allowed kinematically. 
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E. SUSY CP phases and charged Higgs decays H — > th. 

In the neutral Higgs sector, the ratio = BR{hP — > 
bb)/BR{h° tt) is found to be sensitive to the super- 
symmetric loop corrections and to the CP phases. In 
an analogous fashion we may define the ratio = 
BR{H- bt)/BR{H- rvr) and it is also affected by 
SUSY loop corrections, and is sensitive to CP phases. 
Thus the tree level couplings of the third generation 
quarks to the Higgs bosons 

-C = CijhbRHiQ^^ - e^.httaW^Qi + H.c. (151) 

receive SUSY QCD and the SUSY electroweak loop cor- 
rections which produce shifts in couplings similar to the 
case for the neutral Higgs bosons. Thus the general ef- 
fective interaction may be written as 

-£e// = ei,{hb + ShlfbRHlQi + AhlbnW^*Ql 
-Cijih + 5hi)iRHiQ{ + AhitRHl*Ql + H.c. (152) 

We note that in the approximation 

Sh] = 5h), Ah)- = Ah) (153) 

one finds that the above Lagrangian preserves weak 
isospin. This is the approximation that is often used 
in the literature(Carena and Haber, 2003). However, 
in general, the above approximation will not hold and 
there will be violations of weak isospin. In the neutral 
Higgs interaction with the quarks and leptons of third 
generation, we examined 5h\^, ^^6r' ^^"t Ah\. In 
the charged Higgs interaction with these particles we 
should similarly examine 5h1 ^, Ah\ ^, 5h\ and A/i|. The 
latter corrections have SUSY QCD contributions when 
gluinos, stops and sbottoms arc running in the loops 
and SUSY electroweak contributions when neutralinos 
and/or charginos, stops and/or sbottoms are running in 
the loops. The CP violating phases that enter and 
Ahlf are ^3, a At, cha,, and 0^. The phases that appear 
in and A/i^f are ^1,6) "At, aA^, 01 and The 
phases that enter the corrections 8h\ and A/i^ arc the 
same as in (5/i^ and A/i^. One can measure the size of 
the violation of weak isospin by defining 

n = {\Ahl\^ + \5hl\^)H\Ahl\'' + \5hW')-'^ (154) 




(«) (b) 

FIG. 10 Loop diagrams with CP dependent vertices that con- 
tribute to charged Higgs decays into charginos and neutrali- 
nos. 

where 

Bli = -^{hb + ShDe-'^"* sin/3 ^Ahle''""* cos/? 
~{ht + 5hl*)e'''''* cos/3 iA/ife*^"* sin/3 
Bl = ~{ht + Shl*)e'^''* cos/3 ]^Ahfe'^'' sin/3 
-{■^{hb+ 5hl)e-'^'" sin/3 - ^Ahle-'^"' cos/3 
= + 5/i2)e-*Xx/2 sinP+^Ahle-'^-/^ cos/3 

(156) 

and where Bp^ = -Bl^, and Oh = (Xb + Xt)/2. The 
same holds for Xr with b replaced by t. For tan X/ a 
similar expression holds with b replaced by t. The loop 
corrections to the charged Higgs couplings can be quite 
significant. Also the loop corrections can generate signif- 
icant violations of the weak isospin in this sector. 



Similar ratios could be defined for the top and tau, rj 
and Tt- The deviation of these quantities from unity is 
an indication of the violation of weak isospin in the Higgs 

couplings. It is found that such deviations from unity F. SUSY CP phases and charged Higgs decays X*X° 
can be as much as 50% or more depending on the region 

of the parameter space one is in. It is also seen that ^j^^ ^^^^^ ^± ^ ^±^0 sensitive to CP violation 

these measures are sensitive functions of CP violating ^^^^^ ^^^^ ^^^^ ^^^^1 Inclusion of the loop cor- 

phascs (Ibrahim and Nath, 2004). The mteractions of ^^^^-^^^ ^^^^^^^ ^^^^^^^ ^q^^^^ CP phases, 

the charged Higgs are thus governed by the Lagrangian r^^^ ^^^^ ^^^^^ lagrangian for ff±x^X° is 

-C = b{Bt, + Bl^^)tH- + f{Bt, + Bl,^^)vH- + H.c. 

(155) C = ^j^H^X^^jPLxt + i',iHlx\PRxt + H.c, (157) 
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where ^ij and ^- are given by 



G. Effect of CP phases on neutralino dark matter 



^11 



9 
9_ 

%/2 



tan 9wXijV*2 



(158) 



The phases that enter the couplings ^j, and are ^1, ^2 

and 9^. The loop corrections produce shifts in the cou- 
plings and the effective Lagrangian with loop corrected 
couplings is given by 

Hi'ji + Seji)H!x°jPRxt + +AC,,H^*x°jPRxt + H.c. 

(159) 

The phases that enter the corrections A^jj, 5£,ij are ^1, 
^2, ciAti CKAi and 6^. This dependence arises from the 
shifts in the vertices of the charginos with top and sbot- 
toms, charginos with bottoms and stops, neutralino with 
bottom and sbottoms, neutralino with tops and stops, 
W bosons with charginos and neutralinos, Z bosons with 
charginos and neutralinos, charged Higgs with neutrali- 
nos and charginos and charged Higgs with stops and sbot- 
toms. All these vertices enter in the loop corrections. 
Thus Ceff may be written in terms of the mass eigen- 
states as follows 



where 



4 = i^i^'ji + ^£.'oi) sin/3 + 2^^^i cos/3 
+ + COS/3 + ^A^ji sin/3, 

afi = + sin/3 + ^A^^ cos/3 

-^fei + ^^ji) cos/3 - sin/3. (161) 

From the above Lagrangian one can write down the decay 
rate of the charged Higgs into charginos and neutralinos. 



x(0.5(|af/ + |af/)(M|,_ - m^o - m^+) 

-0.5(|4|2 _ \af,\^){2m^^m^o)) (162) 

The charged Higgs decays are found to be more sensitive 
to the phases that enter both at the tree level as well as at 

the loop level such as 6^ (Ibrahim et al, 2004) relative to 
the phases such as aA which enter only at the loop level. 



If the lightest neutralino is the LSP then with R parity 
invariance it is a possible candidate for cold dark matter, 
and in this case the relic density (Chattopadhyay et al., 
1999; Falk and Olive, 1998; Falk et al, 1995; Gomez et al., 
2005) as well as the rates in experiments to detect neu- 
tralinos will be affected by the presence of CP phases 
(Chattopadhyay et al, 1999; Falk et al, 1999a, 2000). 
We give a brief discussion of neutralino dark matter and 
highlight the effects of CP on neutralino dark matter 
analyses. A quantity of interest in experimental mea- 
surements is Vidmhl where Vldm = Pdm/Pc, where pdm 
is the dark matter density, and pc is the critical matter 
density needed to close the universe where 

Pc = SH^/SttGn ~ 1.88 xhlx IQ-'^^gm/cm^. (163) 

Here Hq is the Hubble constant, and ho is its value in 
units of lOOkm/sec.Mpc, and Gjv is the Newtonian con- 
stant. The current limit from WMAP3 on cold dark mat- 
ter is (Spergel et al., 2006) 



i'cdm"' 0.iU'd:0_g Qggg. 



(164) 



In the Big Bang scenario the neutralinos will be pro- 
duced at the time of the Big Bang and will be in thermal 
equilibrium with the background till the time of freeze 
out when they will go out of equilibrium. The procedure 
for computation of the density of the relic neutralinos is 
well known using the Boltzman equations. In general the 
analysis will involve co-annihilations and one will have 
processes of the type 



//, WW, ZZ, WH, 



(165) 



Additionally co-annihilations with staus, charginos, and 
other sparticle species can also contribute. Thus the relic 
density of neutralinos n = rii is governed by the Boltz- 
man equation (Gondolo and Gelmini, 1991; Griest and 
Seckel, 1991; Lee and Weinberg, 1977) 



dn 



-3Hn- 



<7ijv){ninj - nl'^nY) , (166) 



Here cry is the cross-section for annihilation of particle 
types and n^'^ the number density of Xi in thermal 
equilibrium. Under the approximation m/n = nl'^/rf'^ 
one has the well known result 

where (Toff = X^i j ^ijlilji and 7^ are the Boltzman sup- 
pression factors 7j = Explicitly one finds that the 
freeze-out temperature is given by 



Xf = In 



45 



Sn^NfGN 



(168) 
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where Nj is the number of degrees of freedom at freeze- 
out and Gat is Newton's constant. The relic abundance 
of neutrahnos at current temperatures is then given by 



1.07 X lO^GcV" 



dx 



.(169) 



Here Xf = rrii/Tf^ Tf is the freeze-out temperature, 
Mpi = 1.2 X 10^^ GeV, and < av >\s the thermal average 
of av so that 



< (TV >- 



Jo 



dvv^ 



,u /4a; 



(170) 



The diagrams that contribute to < cr?; > include the s 
channel h, and A'^ poles, and the t and u channel 
squark and slepton exchanges. The Higgs boson, and 
sparticle masses are affected by the CP phases of the soft 
parameters. Further, the vertices are also affected. In- 
clusion of the loop corrections to the vertices further en- 
hances the dependence on phases (Gomez et ai, 2004b). 
Specifically the Yukawa couplings of bottom quark and 
neutral Higgs bosons are found to be sensitive to ^3 if one 
includes SUSY QCD corrections in the analysis. A de- 
tailed analysis to study the sensitivity of dark matter to 
the b quark mass and to the neutral Higgs boson mixings 
is given in (Gomez et ai, 2004b). It is found that the relic 
density is very sensitive to the mass of the b quark for 
large tan (3 and consequently also to the GP phases since 
the b quark mass is sensitive to the phases. In Fig. (11) 
we give an exhibition of the relic density and its sensitiv- 
ity to phases. In the analysis presented in Fig. (11), the 
relic density was satisfied due to the annihilation through 
resonant Higgs poles, and one observes the sensitivity of 
the relic density to CP violating phases. The analysis of 
the relic density with inclusion of Yukawa unification con- 
straint with inclusion of CP phases is given in (Gomez 
et ai, 2005). An analysis of relic density in the pres- 
ence of CP phases is also given in (Argyrou et ai, 2004; 
Belanger et ai, 2006a; Falk and Olive, 1998; Nihei and 
Sasagawa, 2004). 

Typical dark matter experiments involve scattering of 
neutralinos of the Milky Way that reside in our vicinity 
with target nuclei. The basic lagrangian that governs 
such scattering is the neutralino-quark scattering with 
neutralino and quarks in the initial and final states. The 
relative velocity of the LSP hitting the target is small, 
and so, one can approximate the effective interaction 
governing the neutralino-quark scattering by an effective 
four-fermi interaction 

^eff = X7m75X97^(-4Pl + BPR)q + Cxx^qm 
+-DX75X'7^g9759 + ExnsXmqqq + FxX'^q^^lbq (171) 

The deduction of Eq.(171) requires Fierz rearrangement 
which is discussed in Sec. (XVI. H) and further details 
are given in Sec. (XVI. I). The first two terms A^B in 
Eq.(171) are spin-dependent interaction and arise from 
the Z boson and the sfermion exchanges. The effect of 
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FIG. 11 An exhibition of the satisfaction of the relic den- 
sity constrains with large phases from the analysis of (Gomez 
et al., 2005). The curve labeled (i) is for the case mo = 1040, 
\Ao\ = 0, tan/3 = 40 , O^u = 2.9, ua = 0,Ci = 1.0, 
^2 = 0.15, ^3 = 0.5, while the curve labeled (ii) corresponds to 
mo = 1080, 1^0 1 = 0, tan/? = 40, 9^ = 0.6, ua = 0,Ci = 0.5, 
^2 = —0.6, ^3 = 1.6. For case (i) EDM constraints are sat- 
isfied when mi/2 = 1250 and for case (ii) they are satisfied 
when mi/2 = 1100. All masses are in units of GeV and all 
angles in radians. 



CP violating phases enter via the neutralino eigen vec- 
tor components and the matrix Dq that diagonalizes the 
squark mass matrix. Then the phases that play a role 
here are 0^, ^1, ^2 and a^, . The C term represents the 
scalar interaction which gives rise to coherent scattering. 
It receives contributions from the sfermion exchange, and 
from the exchange of the neutral Higgs Hi mass eigen- 
states. The term D is non vanishing in the limit when CP 
phases vanish. However, this term is mostly ignored in 
the literature as its contribution is suppressed because of 
the small velocity of the relic neutralinos. In fact the 
contributions oi D,E and F are expected to be rela- 
tively small and could be ignored. A significant body 
of work exists on the analysis of detection rates in the 
absence of CP phases (Arnowitt and Nath, 1996; Nath 
and Arnowitt, 1995, 1997), but much less so with inclu- 
sion of GP phases. Inclusion of the CP phases shows a 
very significant effect of CP phases on the detection rates 
(Chattopadhyay et ai, 1999; Falk et al, 2000). The CP 
effects can be significant even with inclusion of the EDM 
constraints (Gomez et ai; Nihei and Sasagawa, 2004). 



H. Effect of CP phases on proton stability 

GP violating phases can affect the nuclcon stability in 
supersymmetric grand unified models with baryon and 
lepton number violating dimension five operators ((Sakai 
and Yanagida, 1982; Weinberg, 1982). For a recent re- 
view see (Nath and Perez, 2007)). Thus in a wide class 
of unified models including grand unified models, string 
and brane models baryon and lepton number violation 
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arises via dimension LLLL and RRRR chiral operators 
of the form 



and Pp is defined by PpUl 



l^hL = ■J^£abc(PfiV)i.j{f2)kl{uLbidLcj) 
{elk{VUL)al - l^kdLal) + ■■ + H.C . 



(172) 



^5R = -J^eabc{V''r)ij{PVf'^)kl{e'iiiURajURckdRbl + 

.. + H.C. (173) 

Here L^l in the LLLL and L^n is the RRRR lepton and 
baryon number violating dimension 5 operators, V is the 
CKM matrix and fi arc related to quark masses, and 
Pi appearing in Eqs(172) and (173) are the generational 
phases given by Pi = (e*'''* ) with the constraint J2i 7i = 
(i=l,2,3). 

Using the above one generates the baryon and the lep- 
ton number violating dimension six operators by dressing 
the dimension five operators by the chargino, the gluino 
and the neutralino exchanges. The dressing loops con- 
tain the CP phases both via the sparticle spectrum as 
well as via the vertices. This can be explicitly seen by 
elimination the sfermion fields above via the relations 



UiL = 2 j [A^iL^i + Af 

UiR = 2j [A^^R^i + Af (174) 

where Lui — SLj/Sul^, Rui — SLi/Suji^. Here L/ is 
the sum of fermion-sfermion-gluino, fermion-sfermion- 
chargino and fermion-sfermion-neutralino and A's are 
the propagators. A detailed analysis of the specific mode 
p DK^ which is typically the dominant mode in su- 
persymmetric decay modes of the proton is then given by 
the following with the inclusion of CP phases 

\A.,K\''AliA^sf 



1(1 + "^-f + ^^^ ^ii + yl^ + (e-^3;y^ + y^)^^^ + 

oVTlB 

§y?5i,) + l^Dii + yf - ie-^i^y, - y,)Si, 

+§y2%3)\' (175) 



where 



A.,K = {sm2(3M^)-^alP2mcmfV^V2iV22 
{J^{c;di;W)+J^{c; ei;W) 



(176) 



In the above A]:^{Ag) are the long (short) suppression 
factors, D,F, are the effective Lagrangian parameters. 



where U2 is the proton wavefunction 



eabceap < 0Kl''4l''^Jl\p > 
Theoretical 

determinations of /3p lie in the range 0.003 — 0.03 GeV^. 
Perhaps the more reliable estimate is from lattice 
gauge calculations which gives(Tsutsui et al., 2004) 
\pp\ = 0.0096(09)(l|o) GeV^. 

CP violating phases of the soft SUSY breaking sector 
enter in the proton decay amplitude. The CP phases 
enter the dressings in two ways, via the mass matrices of 
the charginos, the neutralinos and the sfermions, and via 
the interaction vertices. Taking account of this additional 
complexity, the analysis for computing the proton decay 
amplitudes follows the usual procedure. This effect is 
exhibited by considering Rt 



R-r 



t{p - 



K+) 



To{p^9 + K+) 



(177) 



where t{p — » P + K'^) is the proton lifetime with CP vio- 
lating phases and To{p P + K~^) is the hfetime without 
CP phases. This ratio is largely model independent. All 
the model dependent features would be contained mostly 
in the front factors which cancel out in the ratio. Since 
the dressing loop integrals enter in the proton decay life- 
time in GUTs which contain the baryon and the lepton 
number violating dimension five operators, the phenom- 
ena of CP violating effects on the proton life time should 
hold for a wide range of models of GUTs. The baryon 
and lepton number violating operators must be dressed 
by the chargino, the gluino and the neutralino exchanges 
to generate effective baryon and lepton number violat- 
ing dimension six operators at low energy. These dress- 
ing loops have vertices of quark-squark-chargino, quark- 
squark-neutralino and quark-squark-gluino. Prom this 
structure one can read the phases that might enter the 
analysis. The chargino one has the phases of 0^, a^, and 
^2- The neutralino vertex has beside the above set, an 
extra phase . The gluino vertex has the set of , , 
^3. Following the standard procedure (Nath et ai, 1985; 
Weinberg, 1982) one can obtain the effective dimension 
six operators for the baryon and lepton violating interac- 
tion arising from dressing of the dimension five operators. 
By doing so and estimating Rt, one finds that this ratio 
is a sensitive function of CP phases (Ibrahim and Nath, 
2000b) . Modifications of the proton lifetime by as much 
as a factor of 2 due to the effects of the CP violating 
phases can occur. It is found also that the CP phase 
effects could increase or decrease the proton decay rates 
and that the size of their effect depend highly on the 
region of the parameter space one is in. 



I. SUSY CP phases and the decay B° 



The branching ratio of the rare process 
is another area where CP violating phase effects arise. 

It is known that the standard model value is rather 
small while in supersymmetric models it can get three 
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FIG. 12 The counter term diagram which produces the lead- 
ing term in amplitude proportional to tan'^ (3 in the branching 
ratio B° l+T 

orders of magnitude larger for large tan/3 (Arnowitt 
et at, 2002; Babu and Kolda, 2000; Baek et al, 2003; 
Bobeth et ai, 2001; Buras et ai, 2002; Chankowski and 
Slawianowska, 2001; Choudhury and Gaur, 1999; Dedes 
et ai, 2002; Huang et ai, 2001; Isidori and Retico, 2001; 
Mizukoshi et ai, 2002; Xiong and Yang, 2002). 

Detecting such large values of would be a positive 
test for SUSY even before any sparticles are found. This 
decay is governed by the effective Hamiltonian 

X {CsOs + CpOp + C'sO's + C'pO'p + CioOio)q (178) 

where C's are the coefficients of the Wilson operators 
O's defined by 

Os mb{d'aPRba)ll, Op ^ mb{d'aPRba)h5l, 
O's = md'id'aPLba)ll, O'p = nid'{d'aPLba)h5l, 

Olo = id'o.l''PLbc.)hf.l5l, (179) 

and Q is the scale where the coefficients are evaluated. 
The branching ratio is a function of the coefficients Cs,p 
and C'g p. In the counter term diagram (see Fig. (12) 
which contributes to this ratio one can find vertices of 
bbHi, sx~t and fl^Hi. The first two vertices are sensi- 
tive functions of the CP violating phases as was explained 
in the different applications above. The phases that play 
a major role here are 0^, ^2 and a^, • Gluino and neu- 
tralino exchange diagrams also contribute which brings a 
dependence on additional phases and ^3. Inclusion of 
these (Ibrahim and Nath, 2003a) shows that the branch- 
ing ratio can vary in some parts of the parameter space 
by up to 1-2 orders of magnitude due to the effect of 
CP phases. A demonstration of the strong effect of the 
phases on B decay branching ratio is given in Fig. (13). 
An analysis of this process using the so called resummed 
effective lagrangian approach for Higgs mediated interac- 
tions in the CP violating MSSM is given in (Dedes and 
Pilaftsis, 2003). 
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FIG. 13 An exhibition of the strong dependence on ua of 
the ratio of the branching ratios B{Bs — > ^'^ fi~)/B{B^ 
/i^/i~)o, where B{Bg jj.'^fi~)o is the branching ratio when 
all phases are set to zero (Ibrahim and Nath, 2003a). The 
curves in ascending order are for values of |ylo| of 1,2,3,4,5. 
The other parameters are mo — 200 GeV, mi/2 ~ 200 GeV, 
tan/3 = 50, Ci = 6 = 7r/4, ^3 = 0, and 61^ = 2. 

J. CP effects on squark decays 

The interactions of qq[Xj and qqi^j do have CP vi- 
olating phases at the tree level. These interactions are 
important for squark decays into fermions and such de- 
cays are expected to show up in the Large Hadron Col- 
lider when squarks become visible. The Lagrangian that 
governs the squark decays is given by 

C = gt{RujPR. + LujPL)x+jbi + gb{RujPR + 
LujPl)x''jU + gtiKujPp + MujPl)x°jU 

+gb{Ku,PR + MujPL)x°jh + H.c. (180) 

where 



v2miy sin /3(cos/3) 

and where 

Lhij = KtV*2Dbii 

Rbij — —{UjlDtli — KbUj2Di,2i) 

Kbij = -V2[l3bjDbu + a*bjDb2i] 

Mb^J = ~V2[abjDbu - JbJDb2^] (182) 

The corresponding quantities with subscript t can be ob- 
tained by the substitution b ^ t, U < > V. The cou- 
plings R and L are functions of the phases 9^, ^2 and 
aAg- The set of phases that enter the couplings K and 
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as in B ^ J/^K^ to 0(A2) (Grossman and Worah, 1997). 
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FIG. 14 A sample of one loop diagrams with CP phase de- 
pendent vertices that contribute to the decay of the stops. 



M is the same above set with an extra phase . The loop 
corrections produce shifts in the couplings as follows 



^eff = 9t{RbijPR + LujPL)X^ji>i 

+gb{RtijPR + LujPL)x''jti + gi{KujPR + MujPl) 

X°jii + gb{KbijPR + MbijPL)x°jbi + H.c. (183) 

where Ruj = iibij + Ai?(,ij where ARhij is the loop cor- 
rection and other tilde are similarly defined. The loops 
that enter the analysis of A's have gluinos, charginos, 
neutralinos, neutral Higgs, charged Higgs, squarks, W 
and Z boson exchanges. The masses of sparticlcs as well 
as the vertices where they enter are sensitive to the CP 
phases. The analysis using the loop corrected lagrangian 
enhances the CP dependence of the masses and the ver- 
tices that already appear at the tree level. Recent anal- 
yses of stop and sbottom decays can be found in (Bartl 
et al, 2003, 2004c; Ibrahim and Nath, 2005) 



K. S — » (pK and CP asymmetries 

Like B — > Xg + ^, the decay B — > (pKs has no tree 
level contribution and proceeds only via loop corrections. 
Thus the process presents a good testing ground for new 
physics since new physics also enters at the loop level. 
An interesting phenomenon concerns the fact that in the 
SM, the CP asymmetry predicted for BcpKs is the same 



The current value of the B — > J/^Kg experimentally 



is 



Sj/^K, = 0.734 ± 0.055 



(184) 



which is in excellent agreement with SM prediction of 
sin 2/3 = 0.7151q;o45. Although currently the experimen- 
tal value for S^Ks (Aubert et al., 2004) 



S^K = 0.50 ± 0.25{stat)toM{syst) 



(185) 



is consistent within la of the SM prediction, its value 
has significantly in the past showing a 2.7a deviation 
from the SM prediction, which triggered much theo- 
retical activity to explain the large deviation (Agashe 
and Carone, 2003; Arnowitt et ai, 2003; Back, 2003; 
Chakraverty et al., 2003; Cheng et al., 2004; Chiang 
and Rosner, 2003; Ciuchini and Silvestrini, 2002; Datta, 
2002; Dutta et ai, 2003; Hillcr, 2002; Kane et ai, 2003; 
Khahl and Kou, 2003; Kundu and Mitra, 2003). 

Although the discrepancy has largely disappeared it is 
still instructive to review briefly the possible processes 
that could make a large contribution to the B 4>Ks 
process. It should be noted that the branching ratio 
Bi?(BO ^ (j)Ks) = (8.0 ± 1.3) X 10"*^ is quite consis- 
tent with the SM result. 

The time dependent asymmetries in B ^ (p^S sje 
defined so that 



A4,K{t) 



r(B(i) ^ (l>Ks) - T{Bit) ^ <j)Ks 



T{B{t) ^ cjyKs) + T{B{t) ^ <j)Ks 
= —C^K cos{AmBt) -h S^K sin(AmBt) (186) 



where S^Ks and C^k g are given by 



4>Ks 



1 + \^4'Ks? o 

2 ) i^ifiKs 



2ImX 



1 + \^(l>Ks. 

where X^Ks is defined by 



1 + \X4,Ks 



A{BO ^ <j>Ks) 



(187) 



(188) 



where /? is as defined in the SM, and 6fi iw any pos- 
sible new physics contribution. Much of the work 
in trying to produce large effects within supersym- 
mctric models has focusscd on generating corrections 
from flavor mixing in the quark sector using the mass 
insertion method(Gabbiani et al., 1996; Hall et al., 1986). 

Thus, for example, an LL type mass insertion in the 
down quark sector will have the form 



d V. 
LL)^J 



iLvi 



(189) 



Here {MJ)ll is the LL down squark mass matrix, Vj^ 
is the rotation matrix that diagonalizes the down squark 
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mass matrix, and rh is the average squark mass. Sim- 
ilarly one defines the mass insertions {S'^fi)ij, {^LRjij 
and {5'jiL)ij. Among the supersymmetric contributions 
considered are the gluino-mediated h sqq with q = 
u, d, s, c, b and Higgs mediated b — * sss. Typically it is 
found that the LL and RR insertions give too small an ef- 
fect but chirality flipping LR and RL insertions can gen- 
erate sizable corrections to B(f)Ks. Thus, for example, 
I('^l_r)23| — 10~^ can significantly affect B (pKs while 
the constraints on B ^ Xgj and AM^ are obeyed. The 
analysis in B — > ^Ks in supergravity grand unification 
with inclusion of CP phases was carried out by (Arnowitt 
et al., 2003) and it was concluded that significant correc- 
tions to the asymmetries can arise with inclusion in the 
trilinear soft parameter A with mixings in the second and 
third generations either in the up sector or in the down 
sector. A similar analysis of asymmetries in B ij' K 
have also been carried out by (Gabrielli et ai, 2005). 

L. T and CP odd operators and their observability at 
colliders 

In the previous sections we have discussed the effects of 
CP violation on several phenomena. The list of CP odd 
or T odd (assuming CPT invariance) is rather large (For a 
sample, see, (Bernreuther and Nachtmann, 1991; De Ru- 
jula et ai, 1991; Kane et ai, 1992; Valencia, 1994).). We 
discuss briefly now the possibilities for the observation of 
CP in collider experiments. First we note that CP phases 
affect decays and scattering cross sections in two differ- 
ent ways. Thus in addition to generating a CP violating 
contribution to the amplitudes, they also a£Fc;c;t the CP 
even part of the amplitudes which can affect the over all 
magnitude of decay widths and scattering cross sections. 
However, definite tests of CP violation can arise only via 
the observation of T odd or CP odd parts. As an exam- 
ple of the size of the effects induced by CP odd operators 
in supersymmetry on cross sections consider the process 
e+e^ tt. Here an analysis in MSSM including loop 
effects with CP phases gives (Christova and Fabbrichesi, 
1993) 

where k (p) are the center of mass momentum of one of 
the initial (final) particles and J is the unit polarization 
vector of one the produced t quarks perpendicular to the 
production plane, c depends on the details of the sparti- 
cle spectrum and can vary significantly depending on the 
sparticle spectrum. The choice c ~ .1 give the correction 
of the T-odd observable to be of size (10^^^) which is 
typically of the same size as the radiative corrections from 
the Standard Model. More generally in e^e~ colliders in 
the process e+e^ X with momenta pi,p2,pa, product 
of the type (^^ x ^j).S,k where is either a momentum or 
a polarization will give a T-odd observable. For example 



one will have T-odd operators of the type (Gavela et ai, 
1989) 

Ti = {pi X P2)-Se-, 

%=p.{S,- xS,+ ) (191) 

More generally with several particles {i = 1, .., n, n > 4) 
one can form T odd operator such as 

ea/SjSPtP^plpt. (192) 

An example of such an operator is the squark decay 
t ^ t + + Xi which can also lead to an observable 
signal at the LHC (Langacker et al, 2007). A study 
of the effects of CP-violating phases of the MSSM on 
leptonic high-energy observables is given in (Choi et al., 
2004a). An efficient way to observe CP violation is via 
use of polarized beams in e+e^ colliders which is of 
interest in view of the proposed International Linear 
collider. A discussion on tests of supersymmetry at 
linear colliders can be found in (Baer et al, 2004) and a 
detailed discussion of tests of CP asymmetries is given 
in (Moortgat-Pick et al., 2005). A number of works 
related to the effects of CP on the Higgs and sparticle 
phenomena discussed in this section arc (Accomando 
et al., 2006; Akeroyd and Arhrib, 2001; Alan et ai, 2007; 
Bartl et al, 2006, 2004a,b; Boz, 2002; Cheung et al., 
2006; Choi et al, 2004b; Ghosh et ai, 2005; Heincmeycr 
et ai, 2004; HoUik et al., 1999, 1998). An interesting 
issue concerns the possibility of expressing CP odd 
quantities in terms of basis independent quantities for 
the supersymmetric case similar to the Jarlskog invariant 
for the case of the Standard Model. Recent works in this 
direction can be found in (Dreiner et al., 2007; Lebedev, 
2003). 

Finally we note that the computation of SUSY phe- 
nomena with CP phases is more difficult than computa- 
tions without CP phases. In Sec. (XVI. J) we give a brief 
discussion of the tools necessary for the computation of 
SUSY phenomena with CP phases. 

XII. FLAVOR AND CP PHASES 

CP violation can influence flavor physics (for recent 

reviews sec (Bigi, 2007; Fleischer, 2006; Schopper, 
2006)) and thus such effects could be used as probes 
of the SUSY CP violation effects. This can happen 
in several ways. This could happen in CP violation 
effects in K and B physics, or if EDMs of leptons are 
measured and turn out to be in violation of scaling, and 
in possible future sparticle decays which may contain 
flavor dependent CP violating effects. Let us consider 
first CP violation in K and B physics. Essentially all 
of the phenomena seen here can be explained in terms 
of the CP violation with a Standard Model origin, i.e., 
arising from the phase Sckm- This means that unless 
some deviations from the Standard Model predictions 
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are seen, the supersymmetric CP violation must be 
small. On the other hand if significant deviations occur 
from the Standard Model predictions then one would 
need in addition to the large CP phases a new flavor 
structure. An example of this is flavor changing terms 
arising from the off diagonal component in the LR 
mass matrix {6ij)LR{d) — {'m\j^{d))ij /m^ (Demir et al., 
2000a,b; Dine et al, 2001, 1993; Khalil and Kobayashi, 
1999; Masiero and Murayama, 1999). 

Further, if one adopts the view point that the entire 

CP phenomena in the K and B system arise from the 
supersymmetric CP phases (Brhlik et al., 2000a; Frere 
and Helen Gavela, 1983) then one will need a new flavor 
structure. But such an assumption appears to be drastic 
since Yukawa couplings arising from string compactiflca- 
tion will typically be complex. However, there are other 
ways in which CP violation can act as strong probes of 
flavor physics and vice versa. For instance, SUSY CP 
effects would be relevant in flavor changing neutral cur- 
rent processes such as 6 ^ s + 7 and in /z ^ e + 7. Also 
if the EDM of the electron and the muon are eventually 
determined and a scaling violation is found, then such 
effects give us a connection between CP violation and 
flavor. Similarly the connection between and CP and fla- 
vor can be obtained from collider data in the decays of 
sparticles. In the following we discuss two speciflc phe- 
nomena where CP and flavor affects can be significant. 
Issue of flavor and CP violation is discussed in many pa- 
pers (Ayazi and Farzan, 2007; Chang et al., 2003; Demir 
and Farzan, 2005; Demir et al, 2000b; Masiero and Mu- 
rayama, 1999). Additional papers that discuss these is- 
sues are(Ellis et al, 2006; Farzan, 2007; Gronau, 2007; 
Pospelov et al., 2006a, b). CP and flavor violation in 
SO(IO) is discussed in (Babu et al., 2005; Chen and Ma- 
hanthappa, 2005; Dutta et al., 2005; Harvey et al., 1980; 
Nath and Syed, 2001). 

A. dpi vs de and possible scaling violations 

The EDM of the muon and the electron are essentially 
scaled by their masses, so that 

d^/de ^ mn/rrie (193) 

The current experimental limits on the muon EDM are 
much less stringent than on the electron EDM, and thus 
it is reasonable to ask if the EDM of the muon could 
be much larger than the EDM of the electron. If so the 
improved experiments on the muon EDM may be able 
to detect it. Thus we explore the conditions under which 
signiflcant violations of scaling may occur. Now we recall 
from our discussion of the EDM of the electron, that 
large EDM for the electron generated by the chargino 
exchange may be canceled by the contributions from the 
neutralino exchange. Thus one possibility in generating 
a large muon EDM is to upset this cancelation for the 
muon case. This appears possible by inclusion of flavor 



dependent nonuniversalities in the soft parameters. To 
make this idea concrete we consider the chargino and 
neutralino exchange contributions to a lepton EDM are 

2 w? . 

di = 9- > m^+ZmfKoKi )A( — ^) 

47rsin2e^m? ^ ^ ^ m? > 

+w^EE^-(^-)wQi^(]t^)(^94) 

fe=l t=l If. If. 

where A, B and k; are deflned earlier, and where r]\^. is 
given by 

rilk = [-y/2{tan0wiQi - ni)Xu + niX2i}D*i,^ - 
-KiX:iiD*i2k]{V2tBxi6wQiXiiDi2k - K/XsjD/uXlQS) 

Here X diagonalizes the neutralino matrix M-^q , and Di 
diagonalizes the slepton (mass)^ matrix. The chargino 
exchange contribution depends on the single phase com- 
bination ^2 + 0^1, while the neutralino exchange contri- 
bution depends additionally on the phase combinations 
+ Ci' ^iid ^iJ- + Non- universalities can be intro- 
duced in two ways: via sneutrino masses which enter in 
the chargino exchange, and via slepton masses that enter 
in the neutralino exchange diagram. One efficient way 
to introduce nonuniversalities in the slepton sector is via 
the trilinear coupling parameter Ai which can be cho- 
sen to be flavor dependent at the GUT scale. In this 
case the cancelation in the electron EDM sector would 
not imply the same exact cancelation in the muon sector 
and significant violations of the scaling relation can be 
obtained. 

Since the violations of scaling arise from the neutralino 

sector we discuss this in further detail. Here the leading 
dependence of the lepton mass arises from n'^, while sub- 
leading dependence arises from the outside smuon mass 
factors in Eq.(194), Thus to understand the scaling phe- 
nomenon and its breakdown we focus on n\^. which can 
be expanded as follows using Eq.(195). 

'Hik = CLoCoXfiD*^^Di2k + boCoXiiX'2iD*i^Di2k 
—KiaoXiiX:u\Dijk\ - KiboX2iX3t\Di2k\ — 

KiCoXuX3i\Di2k\'' + KfXlDiikD^2k (196) 

where ag, bo and cq are independent of the lepton mass. 
The first two terms on the right hand side of Eq.(196) 
are linear in lepton mass through the relation 

Im{D*iiDi2i) = -Im{D*i2Di22) = 

Till 

_]^2 (mo I A; I sin a/ + |^| sin 6*^ tan/3).) (197) 
11 12 

The third, fourth and fifth terms on the right hand side 
of Eq.(196) have a leading linear dependence on the lep- 
ton mass through the parameter k/ and have additional 
weaker dependence on the lepton mass through the diag- 
onalizing matrix elements Dij. The last term in Eq.(196) 
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FIG. 15 An exhibition of the strong flavor dependence via A 
nonuniversahties in enhancing the muon EDM relative to the 
electron EDM in the cancelation region from the analysis of 
(Ibrahim and Nath, 2001b). Plotted are the electron EDM 
(solid line with squares), of the neutron EDM d„ (dashed line 
with plus signs), and of the muon EDM as a function of 
l^ol for the case when tan/3 = 20, mo = 200, mi — 246, — 
.28, C2 = -.51, = -.11, Of, = .4 and a^, = 1.02 where all 
masses are in GeV. The curve with dashed line with triangles 
pointed down is a plot of the muon EDM which have all 
the same parameters as for de and dn except that qa^ = 0.0 
and the curve with dashed line with triangles pointed up is a 
plot of the muon EDM which have all the same parameters 
as for de and dn except that \A^\ = 6.0, and oa^ = —2.0. 



B. SUSY CP phases and the FCNC process B Xsj 

There are other effects of the CP violating phases on 
the phenomenological constraint arises from the measure- 
ment of the rare decay B Xs"f. This decay only oc- 
curs at the one loop level in the standard model (Al- 
tomari et ai, 1988; Casalbuoni et ai, 1993; Colangelo 
et ai, 1993; Deshpande et at, 1987; Dominguez et al, 
1988; Falk et ai, 1994; Grinstein et ai, 1988). The su- 
persymmetric radiative corrections might be of the same 
order of magnitude as the standard model contribution 
(Baer et ai, 1998; Barbieri and Giudice, 1993; Barger 
et ai, 1993; Bertolini et ai, 1991; Bertolini and Vis- 
sani, 1995; Diaz, 1993; Garisto and Ng, 1993; Goto and 
Okada, 1995; Hewett, 1993; Lopez et ai, 1993; Nath and 
Arnowitt, 1994). It is recently been recognized that su- 
persymmetric contributions can receive significant contri- 
butions from the next-to-leading order corrections (NLO) 
which are enhanced by large tan/3. These are typically 
parameterized by e's. In addition to the e's there are 
other two loop (NLO) corrections which, however, are 
small and can be absorbed in a redefinition of the SUSY 
parameters(Carena et ai, 2001; Degrassi et ai, 2000). 
Currently the branching ratio oi B Xg^f is fairly ac- 
curately known cxperimently(Abe et ai, 2001; Aubert 
et ai, 2002a,b, 2005; Barate et ai, 1998; Chen et ai, 
2001) and imposes very significant constraints on model 
building. The current experimental value is 



BR{B Xsj) = (355 ± 2AZ% ± 3) x 10"^ (198) 



is cubic in the lepton mass. However, in most of the pa- 
rameter space considered, the first term in Eq.(196) is the 
dominant one and controls the scaling behavior. Thus for 
the case when all the soft SUSY breaking parameters in- 
cluding A arc universal (i.e., Ai — A) in Eq. (197)), one 
finds that scaling results, i.e., df^/de ~ m^/me- However, 
for the nonuniversal case since the contribution from the 
A parameter is flavor dependent we have a breakdown 
of scaling here. An analysis is given in Fig. (15). This 
breakdown can be seen by comparing d^ for the nonuni- 
versal cases (dashed line with triangles pointed down and 
dashed line with triangles pointed up) with de (solid line 
with squares) in Fig. (15). With the inclusion of nonuni- 
versahties the d^ can be several orders of magnitude 
larger than de . Specifically values of c?^ could be as large 
as (lO"^"* — 10~^'^)ecm and within reach of the proposed 
experiments which extend the search for the muon EDM 
to the range 10~^^ ecm. An enhanced EDM for the muon 
relative to the electron EDM in excess of what scaling 
law allows can be generated with large neutrino mixings 
arising from the See-Saw mechanism (Babu et al, 2000a). 
Another analysis where lepton flavor violations are used 
to generate an enhancement of the muon EDM is given 
by (Feng et ai, 2001). 



as given by the Heavy Flavor Averaging group(Barberio 
et al, 2006). The Standard Model result with QCD 
corrections(Chetyrkin et ai, 1997) including NLO gives 
(Gambino and Misiak, 2001) BR{B X,-f) = (3.73 ± 
.30) X lO^'*. A similar robust prediction for supersym- 
metric models is needed. To analyze the NLO corrections 
for the supersymmetric case (for recent analyses see, one 
has to examine the effective Lagrangian describing the 
interaction of quarks with the charged Higgs fields 
and the charged Goldstones which we display below 
(see, e.g., (Belanger et al, 2002; Demir and Olive, 2002; 
Gomez et al, 2005, 2006)) 



^cff 



V2A% 1 + ett cot ^ 

T/ 1 + efc(")tan/3_ 



e* tan/3 
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where 



(ii) 



as 



O2 = {CL7^bL){SL7nCL 
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FIG. 16 A sample of diagrams with CP dependent vertices 
that contribute to the NLO corrections to the epsilons in 6 — > 
s + 7 decay. There are a total of 20 such diagrams. 



where 



A/i? Jh\ 
ht ht 

^ +cot/3- ^ 



h* 
"■6 



h* 



A/i? Jhl 
.__L+cot/3— ^ 
ht ht 



K 



tan/3— 
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and where e^b and Ctt is given by 



Ahl Jhl 
£66 = -r-^+cot/3— ^, 
hb hb 

^h\ , Jh'i 
ett = ^ + tan/3— ^. 
ht hi 



(201) 



Using the above Lagrangian along with the interaction 
of quarks and W bosons one can write down the contri- 
butions to Wilson coefficients C7 and Cs in the effective 
Hamiltonian that governs the decay b ^ (for further 
details see (Belanger et al, 2002, 2006b; Kagan and Neu- 
bert, 1998, 1999)) 



AC 

ff = --^VZVtbJ2a{Q)0i{Q) (202) 

^ i=l 



Crfs {Qw) 



^ l + e;(s)*tan/3 ^(2) 



3tan2/3 1 + e^tan/? 



where xt and yt are defined by 



Xt = 



yt = 



w) 



Mi 



7,8 (yt) 
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and and are given by 
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(207) 



The C7 and Cg terms receive dominant exchange con- 
tributions from the W, charged Higgs and the charginos. 
The gluino and ncutralino exchange terms can also con- 
tribute. The gluino exchange contributions have also 
been computed (Everett et ai, 2002). However, it turn 
out that in the minimal flavor violation (MFV) sce- 
nario, the contributions from the gluino and neutralino 
exchanges are indeed relatively small. The analyses of 
b ^ s -1-7 in beyond the MFA scenario where generational 
mixings are taken into account have been carried out in 
the work of (Foster et ai, 2005a,b; Hahn et ai, 2005). 
The most complete analyses of B ^ in SUSY with 
the inclusion of NLO' effects is given in (Buras et al., 
2003; Dcgrassi et ai, 2006; Gomez et ai, 2006). Specif- 
ically in the analysis of (Gomez et ai, 2005, 2006) it is 
shown that the e's as well as the decay B Xg^ are 
sensitive to the GP phases. 



XIII. CP PHASES IN u PHYSICS AND LEPTOGENESIS 

Recent experiments discussed later in this section show 
that neutrinos are not massless. In general neutrinos 
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could have cither a Dirac mass, a Majorana mass or per- 
haps a mixture of the two. For a neutrino to have a 
Dirac mass there must be a corresponding right handed 
neutrino to give a mass term of the type mni^Lt^R + H.c. 
On the other hand, one can generate a Majorana mass 
term from just the left handed neutrinos, i.e., a mass 
term of the form vJ^C~^mLi'L + H.c, where C is the 
charge conjugation matrix. For the case of three neu- 
trino species the Majorana mass matrix is in general a 
symmetric mass matrix of the from (Mohapatra et al., 
2004, 2005; Nunokawa et al., 2007) 

(TOee mefj. TOer \ 
"T'MT (208) 
rriTT J 

The Majorana neutrino mass matrix can be diagonalized 
by an orthogonal transformation so that 

V'^M^V = (209) 

where V can be written asV = UK and where the matrix 
U is similar to the CKM matrix and K is a diagonal 
matrix with two independent Majorana phases. For U 
one can use the parametrization 

(C1C3 C3S1 536"*'' \ 

-S1C2 - CiS2S3e'^ C1C2 - SiS2S3e"' S2C3 I 
-S1S2 - c-iC2Sze^^ C1C2 - siC2S3e'* C2C3 / 

(210) 

where ci = cos ^12, C2 = cos 023, C3 = cos 6*13 and simi- 
larly for si, S2 and S3, where 6ij and 8 are constrained so 
that < 0ij < 7r/2 and < (5 < 27r. The matrix K is 
diagonal and can be taken to be 



K=\q e* (211) 
\0 e''^^ / 

Thus we have three diagonal masses, three mixing angles 
and three phases which together exhaust the full 9 pa- 
rameter set of the Majorana neutrino mass matrix. The 
Majorana CP phases do not enter in the neutrino os- 
cillations, and only the Dirac phase 5 does. Thus the 
oscillation probability from flavor Va to is given by 

P{ya - yf)) = 5o,f)-^Y. UciU0jU*jU*pi 

i>j 

xsin^{AmlL/4:E^) (212) 

where Amfj = \mf — mjl. From the solar neutrino 
and the atmospheric neutrino data (Abdurashitov et al., 
1999; Ahmad et al., 2002a,b; Altmann et al, 2000; Am- 
brosio et al., 2001; Fukuda et al., 2000; Hampel et al., 
1999) one finds that the neutrino mass^ differences are 
given by 

Am^^, = (5.4 - 9.5) x lO"'"^ eV"^ , 
Amlt^ = (1.4 - 3.7) X 10-3 eV\ (213) 



A fit to the solar and atmospheric data using the three 
neutrino-generations gives constraints only on the neu- 
trino mass differences and on the mixing angles. One 
has 

AmL = l|m2|^-|min, 
Am2,„ = ||m3p-|m2n, 
sin^ 6*12 = (0.23 - 0.39), sin^ 6*23 = (0.31 - 0.72), 

sin^ 6»i3 < 0.054. (214) 

An interesting aspect of Eq.(214) is that the mixing an- 
gles 6*12 and 6*23 are large with ^23 being close to maximal 
while ^13 is small. This feature was rather unexpected 
and quite in contrast to the case of the quarks where the 
mixings are small. An important point to note is that 
the neutrino oscillation experiments do not give us any 
information on the absolute value of the neutrino masses. 
Other experiments arc necessary to provide information 
on the absolute values such as from cosmology and from 
the neutrinoless double beta decay. Thus from cosmol- 
ogy one has the following upper bound on each species of 
neutrino masses (Hannestad, 2003, 2004; Hannestad and 
Raffelt, 2004; Spergel et al, 2006) 

^ |m^.J < (0.7 - l)ey (215) 

i 

Similarly the neutrinoless double beta decay gives the fol- 
lowing upper bound on the effective neutrino mass |mee| 
(Bilenky, 2004; Klapdor-Kleingrothaus et al, 2001) 

|mee| < (0.2-0.5)6^ (216) 

where (Mohapatra et al, 2005) 

\7n!'J\ = |C0S^ 6*13(177111 COS^ 6*12 -h 

|m2|sin26»i2e^''^i) +sin26ii3|m3|e2*'^^| (217) 

Several scenarios for the neutrino mass patterns have 
been discussed in order to explain the data. One 
possibility considered is that the third generation mass 
is much larger than the neutrino masses for the first two. 

Among these are the following: (i) 1 7)7,^3! >> ItRj^^^jI, 

(ii)l I >> \m^3V (iii) I ~ 

, 1^2,1^31 ll'^i^J ~ I'^i^jll- Neutrino 

oscillations are sensitive to 5 but not to the Majorana 
phases(Barger et al, 2002b, 1980). As is clear from 
Eq.(217), Majorana phases do enter in the neutrinoless 
double beta decay, but an actual determination of CP 
violation in Qv(}f3 appears difficult (Barger et al, 2002a). 

We discuss now briefly the possible determination of 5 
in the next generation of neutrino experiments such as 
N07^A(Ayrcs et al, 2002, 2004) and T2KK (Hagiwara 
et al, 2006). We begin by noting that under the con- 
dition that CPT is conserved, the conservation of CP 
would require P(va Vfj) — P{va Vjj) = 0. In the 
presence of CP violation this difference is non- vanishing. 
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Thus specifically one has(Barger et al, 2007; Nunokawa 
et ai, 2007) 



-IGJsinf 
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) sin( 
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■), (218) 



where E is the neutrino beam energy, L is the oscillation 
length and J is the Jarlskog invariant for the neutrino 
mass matrix similar to the one for the quark mass matrix 



J Sl2Cl2S23C23Sl3C?3 sin ^. 



(219) 



We note that J depends on ^13 and 5 both of which 
are currently unknown and thus one has only an upper 
limit for J so that J < 0.04. Thus the observation of 
a CP violation via Eq.(218) depends on other factors. 
For example, J vanishes if 613 vanishes and thus the 
effect of CP violation via Eq.(218) would be unobserv- 
able. Similarly, if there was a degeneracy in the neutrino 



masses, for example if 



I , then again the ob- 



servation of CP violation via Eq.(218) would be diffi- 
cult. However, aside from these extreme situations the 
process Eq.(218) holds the strong possibility that long 
baseline experiments should allow one to observe CP vi- 
olation due to S in the neutrino sector. Two experiments 
are ideally suited for this observation. One of these is 
NOi^A (Ayres et ai, 2002, 2004) which wiU be 25kton 
liquid scintillator detector placed 810 km away from the 
NuMI neutrino beam in Fermilab (see Sec. XIV). The con- 
figuration will allow runs in the neutrino as well as in 
anti-neutrino mode. The second possibility is the T2KK 
detector (Hagiwara et ai, 2006) which is discussed in 
Sec.XIV. 



A. CP violation and leptogenesis 

As already mentioned in Sec.l, achieving baryon asym- 
metry in the universe requires three conditions: viola- 
tion of baryon number, violation of C and of CP, and 
departure from thermal equilibrium. Quantitative anal- 
yses show that the Standard Model falls short of fulfilling 
these conditions. Specifically, the amount of CP violation 
is found not sufficient. Thus in the framework of the elec- 
troweak baryogenesis the effective CP suppression factor 
that enters is fcp where (Farrar and Shaposhnikov, 1993; 
Shaposhnikov, 1986) 



fcp = T, 



1\l 2 2\i 2 2\ 



I 2 2\i 2 2\i 2 

(m^ - toJ(to6 - ■ma){m^ 



TOd)si2S23S3i sin(5 (220) 



where = sin( 



\j and 9ij are the three mixing angles. 



and 8 is the CKM phase, and Tq is the temperature of 
the electroweak phase transition (EWPT). The EWPT 
is supposed to occur at values Tc ^ 100 GeV, which 
leads to 8cp ^ (10^^^ - 10"^°). A rough estimate of 
baryon asymmetry in EWPT is B ~ 10~*/cp and the 
Standard Model in this case leads to B ~ (lO'^^- IQ-^S) 




FIG. 17 Generation of lepton number asymmetry via decay 
of the right handed neutrino (N) by interference between the 
tree, the vertex, and the self energy loop diagrams. 



which is far too small compared to the desired value of 
B - 10^1°. Additionally there are stringent constraints 
on the Higgs mass which are already in violation of 
the current limits. Analysis of baryogenesis in MSSM 
relieves some of the tension both because there are new 
sources of CP violation, and also because the Higgs 
mass limits are significantly larger, e.g., rrih < 120 GeV. 
However, the analysis requires a significant fine tuning 
of parameters. 

An attractive alternative to conventional baryogenesis 
(for reviews see (Cohen et ai, 1993; Riotto and Trod- 
den, 1999)) is baryogenesis via leptogenesis ((Fukugita 
and Yanagida, 1986). For recent reviews see (BuchmuUer 
et ai, 2005; Chen, 2007; Nardi et ai, 2006; Mr, 2007a)). 
The essential idea here is that if one can generate enough 
lepton asymmetry, then it can be converted into baryon 
asymmetry via sphlcron interactions which violate B + L 
but preserve B — L. Leptogenesis is a natural conse- 
quence of the See-Saw mechanism (Gell-Mann and Slan- 
sky, 1980; Glashow, 1979; Minkowski, 1977; Mohapatra 
and Senjanovic, 1980; Yanagida, 1979) which is a popular 
mechanism for the generation of small neutrino masses 
(see also (Schechter and VaUe, 1980, 1982; Valle, 2006) 
for early work on the See-Saw phenomenology). To gen- 
erate a See-Saw one needs heavy Majorana neutrinos and 
one can characterize the Lagrangian for the Majoranas 

by 



(221) 



where Ni are the Majorana fields, and A are in general 
complex and thus the A terms violate CP. Further, Ln 
violates lepton number and B — L. Thus the Lagrangian 
of the above type has the general characteristics that 
might lead to the generation of baryon asymmetry via 
leptogenesis. The CP violation occurs in the decay of the 
Majoranas because of the overlap of the tree and the loop. 
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One can define a CP asymmetry parameter so that 



ei 



(222) 



For the case of just two Majorana neutrinos the analysis 
of ei gives 



ei = Ci 



Ml Im,{XX'<) 



12 



M^' (AAt) 



11 



(223) 



where C{z) = Ci{z) + C2iz) where (Covi et al, 1996) 



Ci{z) = {Sw)-^^[1 - (1 + z)ln{^^) 



C2iz) = i8n)-'^ (224) 



For the case of two singlets and Mi < M2 one has 



- 3 Ml Jm(AAt)' 
S^^Ih' (AAt)ii 



(225) 



Next consider the case when the initial temperature Tj 

is larger than the mass of the lightest singlet neutralino 
Ni . In this case neglecting the effect of the decays of the 
heavier neutralinos, one can write the Boltzman equa- 
tions that govern the number densities nj^^ and ub-l 
so that (BuchmuUer et al, 2002, 2005; BuchmuUer and 
Plumacher, 2001) 



dn 



dx 



-{D + S){nN,-n''l) 



dn 



B-L 



dx 



eiD{nN, - n^\) - Wub-l (226) 



Here x = Mi/T, and W ~ Tw/iHx) is the washout 
term. The processes contributing to the Boltzman equa- 
tions are the decays, inverse decays, scattering processes 
with AL = 1 and processes with AL = 2, where 
D = Vd/{Hx) includes decays and inverse decays and 
5 = Ts/{Hx) includes AL = 1 scattering. Two pa- 
rameters that enter prominently in the analysis are the 
effective mass mi which is defined by 



mi = 



(AAt)iil <0> |2 



Ml 



(227) 



and the equilibrium neutrino mass m* defined by 



m* 



3^5 



M, 



PI 



(228) 



where (/* is the total number of degrees of freedom ((/* = 
106.75 for SM). Numerically m* ~ IQ-^eV. 

The ratio ifii/m* controls whether or not A^i decays 
are out of equilibrium. When mi < m* (the weak 



washout region), Ni decay is slower than the Hubble ex- 
pansion and leptogenesis can occur efficiently. For the 
case mi > m* (the strong washout region) the back re- 
actions that tend to washout are fast and leptogenesis is 
rather slow. However, even for mi/m* >> 1, a sufficient 
amount of lepton asymmetry can be generated. The so- 
lution to hb-l can be obtained in the form 



nB-L{x) = Wg_j^exp{— I dx'W{x')) — -j-eiK 



(^)(229) 



where n;g_^ = nB-hix = 00), and where k{x) is given 
by 



k{x) 



D diiNi 



D + S dx 



f ea;p(- / dx'W(a;"))(230) 
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The [B — L) asymmetry is converted into baryon asym- 
metry by spheleron processes so that 



^sph 
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where agph is the spheleron conversion factor {usph = 
28/79), and f is a dilution factor / = n^^^^/n* which de- 
pends on the photon production from the beginning of 
leptogenesis till the point of recombination, and numeri- 
cally f= 2387/86. One then has 



r]B — 10 eiKf 



(232) 



Now an upper limit on ei can be obtained assuming that 
iVi decay dominates the asymmetry as assumed above 
with a hierarchical pattern of heavy neutrino masses, 
and assuming that the decay of A^i occurs for T > 10^^ 
GeV. In this case one can deduce, under the assumption 
M1/M2 « 1, the result (Davidson and Ibarra, 2002) 



3 Mi(m3-m2) 



167r 
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With |m3 — m2| < \/ Am 
bound on Mi so that 



.05 eV, one finds a lower 



Ml > 2 X 10^ GeV 



(234) 



This result implies a lower bound on the reheating tem- 
perature, and this bound appears to be in confiict with 
the upper bound on the reheating temperature to con- 
trol the gravitino overproduction for the supersymmet- 
ric case. Consequently several variants of leptogenesis 
have been studied such as resonant leptogenesis (Pilaft- 
sis, 1997; Pilaftsis and Underwood, 2004, 2005), soft lep- 
togenesis (Boubekeur et al, 2004; Grossman et al, 2003, 
2004), and non- thermal leptogenesis (Fujii et al., 2002). 
The type of CP violation that occurs in leptogenesis in- 
volves neutrinos which are Standard Model singlets, and 
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hence have no direct gauge interactions with the nor- 
mal particles, and in addition are very heavy. Thus a 
direct observation of CP violation that enters leptoge- 
nesis would be essentially impossible in laboratory ex- 
periments. However, in unified models CP phases could 
be inter-relatcci across different sectors and thus indirect 
constraints on such phases could arise in such models. 



B. Observability of IVIajorana phases 

In the previous subsection we found that the lepto- 
genesis does depend crucially on the Majorana phases 
(for a review of Majorana particles and their phases, see 
e.g.,(Kayser, 1984, 1985)). However, these phases arise 
from heavy Majoranas and are not the same as the Ma- 
jorana phases that arise in light neutrino mass sector. 
It was noted in our discussion of the neutrino masses 
that Majorana phases do not enter in neutrino oscilla- 
tions which depend only on the Dirac phase. The Ma- 
jorana phases do enter in the neutrinoless double beta 
decay. However, they do so only in a CP even fashion 
and further their observation in the Qv(3[3 appears diffi- 
cult. The question one might ask is in what processes 
the Majorana phases can enter in a manifestly CP odd 
fashion. It is known that one such process is neutrino- 
anti-neutrino {y v) oscillations (Bernabeu and Pas- 
cual, 1983; dc Gouvea et ai, 2003; Schcchtcr and Valle, 
1981). The analysis of (de Gouvea et ai, 2003) sets out 
some simple criteria for their appearance in scattering 
phenomena. Thus consider the amplitude for the process 
X where 



A, 



(235) 



where we have pulled out a common phase factor e*^^ 
so Ai has no phase dependent factor multiplying it, 6 is 
a CP even phase and is a CP odd phase. Then the 
mirror process X has the following amplitude 



= e'^^ {Ai + A2e'^''-'^>) 



(236) 



where Ai^2 are assumed not to contain any CP violating 
effects and are the same in processes X and X. The 
difference ATcp = ~ I^^P is then given by 

ATcp = 4Ai^2 sin((5) sin((?i) (237) 



The above simple analysis points to the following three 

conditions necessary for CP odd effects to arise in the 
process X vs its mirror process X. These are (a) the 
existence of two distinct contribution to the amplitude, 
(b) the two contributions must have a non-vanishing 
relative CP odd phase, and (c) they must also have a 
non- vanishing relative CP even phase. The analysis of 
(de Gouvea et ai, 2003) considers application to the pro- 
cess 



itw- 
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(238) 



for which one has the amplitude 



(239) 



where E is the energy of the intermediate state which 
propagates a microscopic distance L, U is the mixing ma- 
trix, and S depends on the initial and the final states and 
on kinematical factors. For the CP-conjugate process 
1-W+ l+W- one has 



Ax = JZiXiUaiUfSiY^e-'^S 



(240) 



where the combination {XiUaiUpi) is free of the phase- 
convention (Bilenky et ai, 1984; Kayser, 1984; Nieves 
and Pal, 1987, 2003). Limiting the analysis to the case 
of two generations we write 
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cos t/ sm 
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e*"^ 
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Under the approximation Ai = l = A2, \S\ = \S\, a = e 
and (3 = 1-1 this leads to 



ATcp = \Aj,f - \Axf = '^\Sfsm'2e 



4^2 
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(242) 



The above example satisfies all the criterion set forth 

earlier for a CP odd effect to appear. CP odd effects can 
also appear in lepton number violating meson processes 
such as — > TT^n^fi^. Thus, for example, if we write 

Ak+ = e^^-+ {A^K + A2Ke*('-+^-'), 
Ak- = e'«^- {AiK + A2Ke'^^''-'t"'^) (243) 



one will have AF^p 



AF^p DC 4AikA2k sm{6K) sm{<pK) 



\Ak+? given by 



(244) 



AL = 2 contributions do arise with R parity violation in 
supersymmetry and contribute to AF^p. However, the 
effect turns out to be extremely small. Some possible cos- 
mological effects of CP violation in neutrino oscillations 
are considered in (Khlopov and Petcov, 1981). 



XIV. FUTURE PROSPECTS 

A. Improved EDM experiments 

There are good prospects of improving the EDM 
bounds significantly. Thus future experiments may im- 
prove the sensitivity of EDM experiments by an order of 
magnitude or more (Dzuba et al, 2002; Kawall et al., 
2004; Kozlov and Derevianko, 2006) and in some cases 



by a significantly larger factor (Lamoreaux, 2001; Se- 
mertzidis, 2004; Scmcrtzidis et ai, 2004). A recent re- 
view on the current experimental situation and future 
prospects regarding the electron electric dipole moment 
can be found in (Commins and DeMille, book chapter). 
Regarding the neutron EDM a sensitivity at the level of 
1.7 X 10~2^ ecm could be achieved (Balashov et ai, 2007) 
and even a sensitivity of 10~^® ecm is possible (Harris, 
2007). Regarding the EDM of ^^^Hg improved mea- 
surements are in progress and a factor of 3-4 improve- 
ment over the next year or so is possible. Beyond that 
there are various projects aimed at improving the limit 
with diamagnetic atoms, using Xe-129, radioactive Ra 
or Rn. However, all of them are still in the develop- 
ment phase, so when one may expect better limits from 
these experiments is unclear. Regarding the muon EDM 
one proposed experiment (Semertzidis, 2004; Semertzidis 
et al., 2004) feasible at JPARC (Japan Proton Accelera- 
tor Reseach Complex) could extend the sensitivity to as 
much as 10^^^ ecm. However, it appears that the earliest 
muons may become available at JPARC is 2016. How- 
ever, recently another proposal for muon EDM has been 
made where the existing muon beam /lEI at PSI could be 
used. It is claimed that the muon EDM with a sensitivity 
of better than rf^ ^ 5 x 10~^^ ecm within one year of data 
taking is feasible (Adelmann and Kirch, 2006). Currently 
there is also an exploration underway regarding the possi- 
ble determination of the deutron EDM using techniques 
similar to the ones used for muon EDM with the goal 
of reaching a sensitivity of 10~^®ecm(Semertzidis, 2007; 
Semertzidis et al, 2004). 



B. B physics at the LHCb 

LHCb is one of the four detectors at the LHC, the other 
three being ATLAS, CMS, and ELLIS. Of the these AT- 
LAS and CMS are the main particle physics detectors 
dedicated to the search for new physics such as super- 
symmetry or extra dimensions. While the ATLAS and 
CMS can also study B physics their capabilities in this 
respect are rather limited. On the other hand LHCb 
is an experiment which is specifically dedicated to the 
study of B physics. Thus the B mesons produced in col- 
lisions at the LHC are likely to lie in angles close to the 
beam directions and a detector ideal for the study of 
B physics should be able to detect such particles. This 
is precisely what the LHCb is designed to do. Specif- 
ically the detection of charged particle will be accom- 
plished by its Ring-Imaging Cherenkov (RICH) detector. 
The precise identification of the interaction region uti- 
lizes a vertex locator or VELO which can be used for 
B-tagging, and more generally for the separation of pri- 
mary and secondary vertices. The number of B mesons 
produced at LHCb will be enormous. Even a luminosity 
of 10^^cm~^s^^ will lead to a number of bb events at the 
rate of 0(10^^^^^) per year. Thus the LHCb will have an 
unprecedented opportunity to study B physics in great 
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C. Super Belle proposal 

The B factories are an ideal instrument for the study 
of elements of the CKM matrix including the CP phase 
5c KM ■ The analyses provided by the B factories at SL AC 
(BaBar) and at KEK (Belle) have given a wealth of data 
and have improved the measurements of the CKM ele- 
ments. Specifically they have been able to measure time 
dependent CP asymmetries with good precision. Further 
improvements in the measurements of these elements will 
come only with significantly greater luminosity. The Su- 
per Belle proposal aims at achieving that by an upgrade 
of the KEKB collider to a luminosity of 10^^~^^cm~^s~^. 
Such an improvement will also require an upgrade of the 
vertex detector for the Super Belle and specific proposals 
are under study (Kawasaki et al., 2006). 



D. Superbeams, v physics, and CP 

The answer to the question whether or not CP phases 
appear in neutrino physics is of crucial relevance to our 
understanding of fundamental interactions. The observa- 
tion of such phases in the light neutrino sector is possi- 
ble using long baseline experiments and intense beams 
(Diwan et al., 2006; Marciano and Parsa, 2006; Mar- 
ciano, 2001) and its observation will give greater cre- 
dence to the hypothesis that such phases also appear 
in the heavy neutrino sector which enter in leptogene- 
sis. Thus the AIP-2004 study recommends "as a high 
priority, a comprehensive U.S. program to complete our 
understanding of neutrino mixing, to determine the char- 
acter of the neutrino mass spectrum, and to search for 
CP violation among neutrinos" (Freedman and Kayser, 
2004). Such high priority efforts could include improved 
0i//3/3 experiments, and super beams to study neutrino 
oscillations and detect CP phases. Specifically the study 
recommends "a proton driver in the megawatt class or 
above and neutrino superbeam with an appropriate very 
large detector ca,pa,ble of observing CP violation and, mea- 
suring the neutrino mass-squared differences and mixing 
parameters with high precision" . One such proposal is an 
upgraded Fcrmilab Proton Driver (FPD). Such an up- 
grade will improve the study of Ug oscillations by 
a significant factor (Geer, 2006). Thus the current Fcr- 
milab NuMI proton beam has 10^'^ protons at 120 GcV 
(a beam power of .2 megawatts). A secondary beam of 
charged pions is generated from the proton beam, and 
the pions then decay producing a beam of tertiary as 
they propagate along a long corridor to the target 735 km 
downstream. With .2 megawatt of proton beam power 



See, e.g., http;//www-pnp. physics. ox.ac.uk/ Ihcb/ 
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one generates only 10~^ interactions in a Ikt detector 
at the far end. Thus an upgrade of the proton beam to 
deliver several megawatt of proton beam power coupled 
with an upgrade of the detector to 10 kt will significantly 
enhance the sensitivity of the detector to observe possi- 
ble CP effects. A similar idea being discussed is T2KK 
where the far detector is put on the east coast of Korea 
along the Tokai to Kamioka (T2K) neutrino beam line 
(Hagiwara et al, 2006). 



XV. CONCLUSIONS 

We have attempted here to give a broad overview of 
CP violation and the effect of CP phases arising from 
physics beyond the Standard Model. We know that 
CP violation beyond what is allowed in the Standard 
Model must exist in order that one generate the desired 
amount of baryon asymmetry in the universe. We have 
examined the origin of such CP violation in some of 
the leading candidates for physics beyond the Standard 
Model. These include models based on extra dimensions, 
supersymmetric models with soft breaking, and string 
models. Specifically supersymmetric models and string 
models generate a plethora of new CP phases and one 
problem one encounters is that such phases lead to 
EDMs for the electron and for the neutron in excess 
of the current limits. One way to limit to these is to 
fine tune the phases to be small which, however, is 
not satisfactory from the point of generation of baryon 
asymmetry. What one needs is a mechanism which 
allows at least some of the phases to be large while 
suppressing their contribution to the EDMs. One 
possibility is suppression of the EDMs by having a heavy 
sparticle spectrum. However, this possibility puts the 
sparticle masses at least for the first two generations in 
the several TeV range and thus outside the reach of the 
LHC. An alternative possibility of controlling the EDMs 
is the cancelation mechanism which allows for large 
phases consistent with the stringent limits on the EDMs 
from experiment. If the cancelation mechanism is valid, 
then the effect of CP phases will show up at colliders 
in a variety of supersymmetric phenomena. We have 
discussed some of these phenomena in this report. One 
important such phenomenon is CP even -CP odd Higgs 
mixing which would lead to discernible signals at Hadron 
colliders and at a future International Linear Collider 
(ILC). Effects of CP could also be visible in i?° M^A*~) 
in Higgs decays /i° 66, rf and in sparticle decays. 
Dark matter analyses are also affected, specifically the 
detection cross section for neutralino-nucleon scattering. 

The future proposed experiments will investigate CP 

phenomena with vastly increased data. Chief among 
these is the LHCb experiment which is dedicated to the 
study of the B mesons. The proposed Super Belle will 
further add to these efforts. These will pin down the 
CKM matrix elements to a much greater precision than 



BaBar and Belle, and may shed light on the possibility 
whether or not new sources of CP violation arc visible. 
However, if the sparticles are indeed observed, as one 
expects they will be, then a study of their branching ra- 
tios is likely to put significant limits on CP phases from 
sparticle decays. 
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XVI. APPENDICES 

A. Chargino and neutralino mass matrices with phases 

Here we give some details on the diagonalization of 
the chargino and neutralino mass matrices which arc in 
general complex. These appear in the analysis of Sees. 
(IX) and (X). We consider the chargino mass matrix 
first. Here we have 

M^^( rh, V2mw sin I3\ 
\V2mwcosp fi ) 

The chargino matrix Mc is not hcrmitian, is not sym- 
metric and is not real since /i and ra^ are complex. For 
simplicity we analyze its diagonalization for real m2 and 
complex \i. Generalization for complex rhi and yt, is 
straightforward. Mq can be diagonalized by using the 
following biunitary transformation 

U'^McV'^ = Md (246) 

Here V and V arc hcrmitian matrices and Md is a di- 
agonal matrix which, however, is not yet real. V and V 
satisfy the relation 

V[mIMc)V-^ = diag(|m^+|2, |m^+|2) 

= [/'*(McM^)(;7'*)-i (247) 

We may parameterize U' so that 

cos^ sin%e*\ 
^ - \^_sin|-e-<^i cos I- ) 
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where 

tan^i = 2V2mw{'rht ~ l/^l^ ~ 2m^ cos2f3)~^ 
(m2Cos^/?+ |/u|^sin2/3+ |/i|m2sin2/3cos6'^)^ (249) 

and 

tanc^i = |/u| sinOfi sin/3 (m2 cos/3 + |/u| cos 6^ sin/3)~^ 

(250) 

Similarly we parameterize V so that 



V 



cos ^ sin ^6-^"^^ 



■ sm ^e" 



cos 



^2 



(251) 



where 

tan ^2 = 2\/2mw{rh2 — \nf + 2m^ cos 2/3)""^ 
(m2sin^/3+ |/i|^cos^/3 + |/i|m2sin2/3cos^^)^ (252) 

and 

tan 02 = — l/i] sm9fi cos/3 (m2 sin/3 + cos6^ cos/3)~^ 

(253) 

We wish to choose the phases of U' and V so that the 
elements of Md will be positive. Thus we define U = 
H X U' where 



H={e 



(254) 



where 71,72 are the phases of the diagonal elements of 
Mn in Eq.(246). With the above choice of phases one 
has 



U*McV ^ = diag{m^+,m^+) 
Our choice of the signs and the roots is such that 

|2 , o™2 



(255) 



^^{rh +)(m +) 



= l[rhl + \l^{' + 2mU{+)i-) 

^[{rhj - + 4m^ cos2 2/3 + 4to^ 

(m^ + + 2m2|/i| cos 6*^ sin 2/3)] ^ 



where the sign chosen is such that rh+ < m + if 

ml < \iJ.f + 2m^ cos 2/3. 
For the neutralino mass matrix M^o one has 



(256) 
(257) 



-Mzswcf! Mzswsfi \ 



mi 


-Mzswcp Mzcwcp 







(258) 



In the above sw = sin 9w, S/3 = sin/3 where 9w is the 
weak angle, and = cos /3, and S/s = sin /3. The matrix 
M-^o is a complex non hcrmitian and symmetric matrix, 
which can be diagonalized by a unitary transformation 
such that 



(259) 



B. Squark and slepton mass^ matrices with pliases 

In this appendix we give details on the diagonaliza- 
tion of the squark and slepton mass matrices that appear 
in Sees. (IX) and (X). We consider the squark (mass)^ 
matrix 



Mi 



-'"g21 -''^922 



(260) 



For the up squark case one has 



^lii =Ml+ml + M|(- - Q„s^) cos 2/3 



Mi^2 = my,iAl- n cot 13) 
^^21 = ™«(^« - cot/3) 
^l22 =ml+ + MlQuslv cos 2/3 (261) 

Thus the squark mass^ matrix is hermitian and can be 
diagonalized by the unitary transformation 



Z?tM|D„ = diag(M?i,M?2) 
where one parameterizes D„ so that 



cos ■ 



COS ■ 



(262) 



(263) 



Here M?2i = |-^|2ik"'^" ^^'^ choose the range of 6u 

so that ^ < 6ij, < f where tan6i„ = . The 

eigenvalues and M?2 can be determined directly 

from Eq.(260) so that 

K'(i)(2) = ^(^»n+^«'22)(+)(-) 



[(M?ii-M?22)'+4|M?2ini 



(264) 



The (+) in Eq.(264) corresponds to the case so that for 
*^lii > -^122 one has M|j > and vice versa. For 
our choice of the 9u range one has 



tan^„ = 



2m„|A„mo - [I* cot/3| 



and 



sm 0„ = 



/tf 2 _ /If 2 

mo|A„| sinct^ + sin 6*^7?,; 
|mo A„ - /X* cot /3| 



(265) 



(266) 



where i?„ = cot /3. The analysis for the down squark case 
proceeds in a similar fashion with the following changes 

^lii = ^1 + '^'^^ -Mzi\+ QdSw) COS 2/3 
Ml^=mdiA*i- ^ltl^np) 
^hi = "^diAd - /i* tan/3) 
MI22 =m% + + MlQdsl^ cos 2/3 

(267) 
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The other changes are the modification of expressions for 
04 and (pd- They read 



2md\Aumo — ji* tan/3| 

dll 422 



and 



sin <j)d = 



mo\Ad\ sin ad + sin0^Rd 



(268) 



(269) 



\mo Ad — tan l3\ 
where Rd = tan (3. Finally for the case of the slectrons 

Mill =Ml+ml~Ml{^- sir) cos 2/3 

M|i2 = We (A* - /itan/3) 
M|i2 = me{Ae - jJL* tan/3) 
^6-22 ='fn% + ml- M|s^ cos 2/3 (270) 

Expressions for 9^ and 0e are identical to the case of the 
down quark with the replacement of d by e. 



C. RG evolution of electric dipole, color dipole and purely 
gluonic operators 

In this Appendix we discuss the renormatization group 
(RG) evolution of the EDMs discussed in Sec.(X). As 
discussed in the text, there are three competing operators 
that contribute to the EDM of the neutron. These are 



6 



(271) 



The one loop RG evolution of the electric dipole and 

of the color dipole operators can be easily obtained us- 
ing their anomalous dimensions since these operators are 
cigcnstatcs under the rcnormalization group. Evolving 
these operators from a high scale Q = Mz to a low scale 
^ one finds 



o,(/x) = r-^*/'^o,(Q) 



(272) 



where 



gsiQV 

IE = 8/3 



7c = 



(29- 

: (33 



2Nf)/S, 



(273) 



where Nf is the number of light quarks at the scale /i. 
Regarding the purely gluonic dimension six operator it 
obeys the following renormalization group equation(Boyd 
et al, 1990; Braaten et al, 1990a,b; Dai et al., 1990; 
Weinberg, 1989) 

^^§J;0G = ^hG0G-6^m,{^,)0,c), (274) 



where jg = ^3 — 2A^/. The gauge coupling and the 
running quark mass satisfy the RG equations 



and 



where 7„ 



M^5.(m) = -/3^^5.(m), 



M^™g(M) = 7m^^™g(M)- 



(275) 



(276) 



-8. The above operators contribute to the 
CP violating Lagrangian multiplied by coefficients which 
must cancel their /i dependence. This allows one to ob- 
tain for the coefficients the following relations 

rf(B,C,G)(^) ~ rT(«,c.G)//3rf(i5,C,G)(Q)^ (277) 

where Q is the high scale. In implementing the RG evo- 
lution one uses the matching conditions due to crossing 
the heavy thresholds for q=b,c. Thus, for example, 

. ^ I rt — 1 nrrt > i i /t — i iv-r\ i ^ ^ ^ 



dP(m-) = 1 (278) 

OTT m„ 



Using this technique one can evolve the EDMs from the 
electroweak scale Q = Mz down to the hadronic scale 
fj.. A more up-to-date discussion of the RG evolution of 
operators including the mixings between the electric and 
the chromoelectric operators is given in (Degrassi et al., 
2005). 



D. Satisfaction of the EDM constraints in the cancelation 
mechanism 

Here we give some examples of the parameter points 

where the cancelation mechanism discussed in Sec.(X.D) 
works to produce de, dn and dng consistent with the 
current limits. Table 1 gives three sets of points (a), 
(b) and (c) for which the corresponding EDMs de, dn 
and Chq are listed in Table 2 where Chq is related to 
the dPd,d'^u,dPs by C„g = \dPd - dPu - 0.012rfC?,|. 
Using the experimental constraints on dug one obtains 
the following constraint on Cug 



Cug < 3.0 X lO^^^cm. 



(279) 



The values of Cng listed in Table 2 are consistent with 
the above experimental constraint. 

Tal)lc 1. Thiw parameter sets willi .In in uuits of mo. 



case 


mo, mi, |Ao| 


aA, 6: ^2, ?3 


(a) 


200,200,4 


1,. 5, .659, .633 


(b) 


370,370,4 


2, .6, .653, .672 


(c) 


320, 320, 3 


.8, .4, .668, .6 



Table 2. Electron, neutron and Hg EDMs. 



case 


de (ecm unit) 


dn (ecm unit) 


(cm unit) 


(a) 


1.45 X lO"'-''^ 


9.2 X 10^^^ 


7.2 X 10-^^ 


(b) 


-1.14 X 10-^' 


-7.9 X 10"^' 


2.87 X 10"^'=' 


(c) 


-3.5 X 10"^' 


7.1 X 10"^' 


2.9 X 10"^" 
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E. Combination of CP phases in SUSY processes 



where 



The various phenomena discussed in Sees. (IX) and (X) 
involve several specific combinations of CP phases. Below 
we exhibit these combinations. 

Table 3: Examples of CP phases in SUSY phenomena 



SUSY Quantity 


Combinations of CP phases 


p ^ UiK^ 


6,2,3 + Ol, aAt^t + 0\ 


s + 7 


"At.,,,,, + ^1, 6,2.3 + Oi 


mixing and spectrum 


"At,b,, 6,2 +0l 




"At, 5 + 6*1, 6,2 + di 


<7m-2 


6,2+6*1, a A, +6*1 




"A, +6*1, 6,2,3 + Ql 


Dark matter 


"A, +6'l, 6+^*1 




6 +6'i, aA,,t +6*1, 6+6*1 


de (dij,) 


6,2 + 6*1, +6'i(aA„ +6*1 ) 


dn 


6,2,3 + 0\, aA^i + ^l,"Adi + Q\ 



In the Table (3) Q\ is defined so that Oi = 6^ + Or and 
the rest of phases are defined as in Eqs.(31) and (33). 



yj2Mw cos (3 



(285) 



Next we discuss the neutralino exchange contributions to 
ttfj^. These are given by 



(286) 



where 



11 _ m^aEM ^ ^ 1 , k ("^l ^ 



27r sin"" 9w ^ M^o 

j—i k=i 



-^3 Xi Xi 



and 



(287) 



4 2 

EE- 

j=l fe=l '^X'i 



a.. . o . > > .i2iT72-'T72-J' 



47r sin^ ew^.t^. '^^ Ml, 



Here 77^^ is defined by 



(288) 



F. Details oi — 2 analysis in SUSY with CP Phases 

Here we give further details of the analysis of dis- 
cussed in Sec. (XI. A) but limiting ourselves to the case 
when the muon mass can be neglected relative to other 
masses. The chargino exchange contribution is given by 



'V2 

(\/2tan6'vF-'^ii-D2fc + n^X^jDik) 



(289) 



and X^- is defined by 



af + 



(280) 



where for a^} and a^,^ we consider now the limit where 
I^{a,l3) and /4(a,/3) that appear in Eq.(113) have their 
first arguments set to zero. In this case one has 

^3(0, a;) = -\Fi{x), h{Q,x) = -\Fi{x) (281) 
where 



2M^ cos2 /? 

2 
1 



+ - tan^ 6w\Xi3\{\DikV + ^D^kV) 

|2 1 n _ |2 I j.„„ /)___ In _ |2 1 



+ -|X2,fl-Difer + tangly |-Diferi?e(Xi,X2*,) 
+ ^?^^i?e(X3,Xr,iPi,I?2\) 



Mi^A cos /3 



MvF cos /3 



Re{XijXip^kD*2k)- 



(290) 



i^3(x) 



(x-iy 



r(3a;2 - 4a; + 1 - 2xHnx) 



Fi{x) = 



{x-iy 



-{2x^ + 3a;2 - 6a; + 1 - 6a;^Zna;).(282) 



In the limit when the muon mass is neglected relative to 
other masses and the first argument in the double integral 
is taken to be zero one finds a simplification of the form 
factors so that 



In the limit considered above one has the following ex- 
plicit expressions for the chargino contributions 



/i(0,a;) = ii^i(a;),72(0,a;) = ^^^2(2;), (291) 



,21 



and 



47rsin^^VK ^ 



where 



{x-l) 



— g(l — a; + 2a;/na;), 



(292) 



22 'm-n'^EM 



1 

»=1 Xi xj 

(284) 



and 



F2{x) = TT4 (-a;^ + 6a;^ - 3a; - 2 - 6a;Zna;). (293) 

[x — i) 
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G. Stop exchange contributions to Higgs mass^ matrix. 



H. Fierz rearrangement relations involving Majoranas 



For completeness we give here an analysis of the one 
loop contributions from the stop sector with inclusion of 
CP violating effects in the analysis of CP even-CP odd 
Higgs mixings discussed in Sec.(XI.B). The contribution 
to the one loop effective potential from the stop and top 
exchanges is given by 



Ml 



a=l,2 



-12mt{log'-^ '-)) (294) 



Using the above potential our analysis for A^^^ gives 



2| |2(| ^t|cos7t - \iJ,\cotl3)'^ 

(mj -m? )2 J^vH,^-t, 



/2(rn? ,m|) 
(295) 



A22t = -2/?fc,m? 



(m? — m? )2 

^ tl 12' 



- \fi\cot p COS jt? 
? 

2 2 

/2(m|, m|) + 2l3hMln{^^) 



\At\[\At\ - \iJL\cotl3 cos-it] 



(m? — rn? 

^ ti ' 



m~ 

ln{^) (296) 



2 MU M COS 7t - |M|cot/3] 

/ 2 2 \ V 2 / 

— m~ ) mf 



^i2i = -'^Phtmt 

2 r / 2 2 

|/i||At|[|At| cos7f - |/;i|cof/3][|>lt| - \ iJ,\cot (3 cos jt 



(rrij - m? )2 



(297) 



^i3t = -2/?fc.m: 



2 iM pl^tl sin7t[|A^|cot/3 ~ cos7t] 
sin/3(m? — m? )2 

f2«,ml) (298) 



A - 9/^ rr,2|„iu |2 sin7i(|At| -|^|cot/3cos7i) 



sin/3(m|^ — m|^)^ 

sin7i 



/2(m^ ,mn + 2/3;,, . Ii^;^^)(299) 

^ tl 12' t2 



and 
A 



33* 



m^|Mp|Apsin^7t 2 2 
■ 2P?it 2 ^/ 9 5rT7 /2 y^ii ' "^^2 ^ ^"^^^^ 



* sin^ /3(to? — m? )2 " 



In the analysis above the D terms of the squark (mass)^ 
matrices are ignored to obtain approximate independence 
of the renormalization scale Q similar to the analysis of 
(Carena et al, 2000; Demir, 1999). 



Fierz rearrangements are known to be very useful 
when manipulating interactions involving four fermions. 
Specifically such Fierz rearrangements are needed in the 
analysis of Sec. (XI. G). Here we give these relations for 
the case when two of the fermions are Majoranas (such as 
neutralinos) and the other two are quarks. Thus any four 
Fermi interaction with two Majoranas and two quarks 
can be written as involving the following combinations 



X7''75Xg7/i75g, XlhXm, XXQlaq- (301) 

For convenience define the 16 gamma matrices as follows 

= {I,70,i7',n%5,y75,75,^<70\a'^} : i,i = 1 - 3 

(302) 

with the normalization 

tr(T^r^) = 4(5-^^ (303) 



The Fierz rearrangement formula with the above defini- 
tions and normalizations is 



(tz-ir^n2)(«-3r^«4) = ^F^^(^z-ir^w4)(w3r^w2) 

(304) 



C,D 



where Uj are Dirac or Majorana spinors and 

1 



(305) 



and where the -|-ve sign is for commuting u spinors and 
the -ve sign is for the anticommuting u fields. In our 
case we have to use the -ve sign since we are dealing with 
quantum Majorana and Dirac fields in the Lagrangian. 
We give below the Fierz rearrangement for four combi- 
nations that appear commonly in neutralino-quark scat- 
tering. These are 

1_ _ 1_ _ 1_ „ _ 
XQQX = --^XXQQ - ^X75Xg'759 + ^X7'^ 75X97/* 75 9 

_ _ 1- u - 1- - 1- - 

X75OTX = ^X7'^75X97/*9 - ^XX975Q' - ^X75X99 

1_ „ _ 1_ _ 1_ _ 

XQQ15X = -^X7 75X97/^9 - ^XX9759 - ^X75X99 

_ _ 1_ _ 1_ _ 1- u - 

xihqqibX = --^xxqq - -^xisxqjsq - -^xi'^i^xqiixi^q 

(306) 

The metric used above is ?7,j;y = (1,— 1,— 1.— 1), and since 
x's are Majoranas we have used the properties X7mX = 
and X(J^,^,x = 0. 



I. Effective four-Fermi interaction for dark matter 
detection with inclusion of CP phases 

In this appendix we give a derivation of the four fermi 
neutralino-quark effective Lagrangian with CP violating 
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phases given in Sec. (XI. G). Wc begin by discussing the 
squark exchange contribution. From the fundamental su- 
pergravity Lagrangian of quark-squark-neutraUno inter- 
actions 

-£ = q[CqLPL+CqRPRWl+€C'qLPL+C^RPR]Xq2+H.C. 

(307) 

the effective lagrangian for g — x scattering via the ex- 
change of squarks is given by (Chattopadhyay et al., 
1999; Falk et al., 1999a). 



B = 



4M^ 



[\X30f - \X4of]egSin'^ew + 



\C,l\' 



A{Ml - Ml) 4(M| - Ml) 



(313) 



The terms C, D, E and F receive contributions from 
sfermions and from neutral Higgs and can be calculated 
using similar techniques. 



91 X 

q[CgLPL + CgRPulx 



1 



+ + C^rPlMC^lPl + C^rPr]x 

(308) 

where 

CqL = V2{aqoDgll - "fqoDg2l), 
CqR = V2{0qoDqii — SqQDq2l), 
^'qL = V2{aqoDqi2 - 7qoDq22), 

C'qR = V2{l3qoDql2 ' SqoDq22), (309) 

and where a, (3, 7, and 6 are given by 



a 



u{d)j 



9^u(d)X4(3)j 

2mw sin /3(cos/3) ' 



Pu{d)3 — (iQu{d)]X[* + 



cosOw 

X (T3M(<i) - Qu{d) sin^ Ow), 

-.n y gQu(d) sin' Ow ^, 

lu{d)j - eQu{d)jXij X2j, 

^ -9mu{d)X^s)j 
~ 2miysin/3(cos/J)" 

Here g is the SU{2)l gauge coupling and 



X't^ = Xij cost 

Xij sin 6w + X2j cos 6w 



i-ij — yi- ij COS aw + X2j sin 9w , 



(310) 



(311) 



The effect of the CP violating phases enter via the neu- 
tralino eigenvector components Xij and via the matrix 
Dqij that diagonalizes the squark mass' matrix. 

Using the Fierz rearrangement one can obtain now the 
coefficients A, B, C, D, E and F that appear in Eq.(171) 
in a straightforward fashion (Chattopadhyay et al, 1999; 
Falk et al., 1999a). The first two terms {A, B) are spin- 
dependent interactions and arise from the Z boson and 
the sfermion exchanges. For these one has 



A = 



^[1X30 r - \XMTm - eqSmHw] 



\CgR? 



\CqR? 



A{Ml-Ml) mi -Ml) 



(312) 



J. Computational tools for SUSY phenomena with CP 
phases 

The numerical analysis of supersymmetric phenomena 
with CP phases is significantly more difficult than for 
the case when the phases arc absent. First most nu- 
merical integration codes for the renormalization group 
evolution, sparticle spectra and for the analysis of sparti- 
clc decays and cross sections are not equipped to handle 
phases. Second any physically meaningful set of param- 
eters which include phases must necessarily satisfy the 
stringent EDM constraints which also require care. A 
significant progress has been in this direction by the so 
called CPsuperH(Lee et al., 2004), which is a Fortran 
code that calculates the mass spectrum and decay widths 
of the neutral and charged Higgs bosons in MSSM with 
CP phases. Obviously there is significant room for fur- 
ther progress in this area. 
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